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Analysis of Symmetrical Waveguide Junctions 
David M. Kerns 


Previous theory of consequences of symmetry in waveguide junctions has made limited 
use of mathematical tools available in the theory of group representations and has been 
limited to the consideration of nondissipative junctions. In this paper group-theoretical 
technique is utilized more fully (in much the same way as in the analysis of the vibration of 
symmetrical molecules) in the formulation of an improved and general technique for the 
analysis of symmetrical waveguide junctions subject merely to the requirement of linearity. 

A waveguide junction, for the purpose of this paper, is a linear electromagne tic system 
possessing ideal w aveguide leads, and is considered to be subject to excitation solely through 
the effects of nonattenuated modes in the waveguide leads. Under the conditions of the 
problem, an arbitrary electric (or magnetic) field in a waveguide junction is expressible 
linearly in terms of a finite number of linearly independent electric (or magnetic) basis fields. 
From any given ordered pair of electric (or magnetic) basis fields one can in principle calculate 
a complex number—an element of the admittance (or impedance) matrix characterizing the 
junction (relative to the choice of basis fields). The geometric concept of rotation and 
reflection of fields (and structures) is discussed in terms of a rotation-reflection operator, and 
the symmetry of a junction is characterized by a group of rotation-reflection operations under 
which the structure is invariant. A general procedure is given for the construction of a basis 
in which the basis fields transform according to irreducible representations of the symmetry 
group involved. Such basis fields are said to be of particular symmetry species and from the 
special properties of such fields follow the physical results, of which perhaps the most 
conspicuous is the vanishing of the matrix element between two fields of distinct symmetry 
species, 


I. Introduction A fA 

el 
This paper is concerned with symmetry properties | - 7 P il uy a 

in “waveguide junctions.” A waveguide junction, | f | v 
for the purpose of this paper, is a linear electro- 
magnetic system possessing ideal waveguide leads 
and is considered to be subject to excitation solely 
through the effects of nonattenuated waveguide 


modes in the waveguide leads. The electromagnetic 
boundary-value problem presented by a waveguide | 
junction is, in general, impracticably difficult to 
solve. Nevertheless, important information con- 
cerning the characteristics of a waveguide junction 

in its primary function as a device for transferring | 


power from one waveguide-mode to another is 
derivable with relatively little labor from general 
properties, such as reciprocity, losslessness, and, in 
particular, symmetry. Many waveguide junctions 
used in microwave practice do in fact possess useful 
and interesting properties in virtue of symmetry. A 
few simple examples of such junctions are shown in 
figure 1. 

The literature on the present subject is not exten- 
sive. The book, Principles of microwave circuits, 
contains a valuable and fairly comprehensive treat- 
ment applying to nondissipative waveguide junctions. 
A report by Slater® is mainly concerned with the 
analysis of nondissipative T-junctions having essen- 
Ui ally a single symme try element. A paper by 

. Nemes q 4 . : 
——. pm rete Bley lr © > Bao sais general theory of the consequences of samatry a 
ati> euaidien of a Gietatenahelien eatuedh waveguide junctions of a general class: in the interior 
as , 6 * | of a waveguide junction media that may be non- 

A dissertation submitted to the faculty of the Graduate School of Arts and homoge neous and anisotropic are pe ‘rmitted; dissi- 
Sciences of the Catholic University of America in partial fulfillment of the ° -r } } 
requirements for the degree of Doctor of Philosophy pation, by reason of finite conductivity or radiation 

?C. G. Montgomery, R. H. Dicke, and E. M. Purcell, Principles of microwave to infinit Vv (or both), is pe rmitted: fulfillme nt of the 


circuits, chapter 12 (McGraw-Hill Co., New York, N. Y., 1948 


J. C. Slater, Technical Report No. 37, Electronics Research Lab., Mass. Inst reciprocity condition is not re quired, Linear be- 
lech. (Cambridge, Mass., 1947 ) 


‘'M. Chodorow, E, L. Ginzton, and J. F. Kane, Proc. IRE 32, 634 (1949 | havior and freedom from internal sources are as- 
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Ficgure 1. Simple junctions possessing symmetry. 


There is to be found also the occasional and limited 
use of symmetry arguments in discussions primarily 
concerned with other matters. 

The analytical technique employed in the book 
referred to in footnote 2 is partly formulated in 
general terms and partly indicated by the consid- 
eration of a series of examples. The technique used 
involves the restriction to nondissipative junctions 
as an explicit condition. The discussions contained 
in the papers referred to in footnotes 3 and 4 are of a 
more or less specific nature and, in the form given, 
are likewise subject to the restriction of no dissipation. 

The object of the present paper is to develop a 








(the former being essential, the latter non- 
A broad outline of the discussion follows. 


sumed 
essential). 

The needed physical and mathematical formula- 
tion of the electromagnetic problem is given in sec- 
tion Il. Under the conditions of the problem, an 
arbitrary electric (or magnetic) field in a waveguide 
junction is expressible linearly in terms of a finite 
number of linearly independent electric (or mag- 
netic) basis fields. From any given ordered pair of 
basis fields of the same kind (electric or magnetic) 
one can in principle calculate a complex number 
an element of the admittance (or impedance) 
matrix of the junction—which is a measure of the 
field of opposite kind associated with the second (or 
first) of the given pair of fields. The characteriza- 
tion of a junction by means of matrices is relative to 
the choice of basis fields; the basis fields first chosen 
might be described as “‘simple with respect to excita- 
tion.”” Formulas for change of basis are given. 

The geometric concept of rotation and reflection 
of fields (and structures), taken up in section III, is 
discussed in terms of a rotation-reflection operator 
applicable to tensor point-functions.” The sym- 
metry of a waveguide junction is characterized by a 
group of rotation-reflection operations under which 
the structure is invariant. 

A method of symmetry analysis and the results 
obtained for the class of problems considered are 
presented in general terms in section IV (with some 
further results in appendix 2). The basis fields set 
up in section I] do not necessarily exhibit particu- 
larly simple transformation properties under opera- 
tions of the symmetry group; however, as is shown, 
it is possible linear combinations of the 
original basis fields to form new basis fields that do 
exhibit special transformation properties. Such 
fields are said to be of particular symmetry species. 
(If, for example, the symmetry group consists of 
only two operations, the identity and reflection in a 
plane, say, the two possible species are the familiar 
“even” and “odd.” The general definition of sym- 
metry species is provided by the theory of group 
representations, which theory indeed provides the 
natural mathematical the analysis.) An 
important property of fields of the new basis is the 
vanishing of the matrix element between two fields 
of distinct symmetry species. 

As might be expected, the method used here is in 
some respects very similar to methods used in the 
analysis of the vibration of symmetrical molecules. 
However, because dissipation (as well as failure of 
reciprocity) is permitted in the waveguide problem, 
the main part of the analysis here is formulated 
without reference to the question of eigenvalues of 
the matrices of a junction. The eigenvalue problem 
is discussed briefly in the latter part of section IV. 

Three 


to select 


tools for 


illustrative examples are considered in 
section V. In one of these examples the theoretical 
results previously obtained (footnote 4) for the 
waveguide Wheatstone bridge are presented in a 
more general context. 


Tensor is used in the general sense, a scalar being a tensor of rank 0, etc. 


II. Electromagnetic Formulation 


A waveguide junction may be described brief 
as an electromagnetic system comprising an arbitray 
number, n, of ideal waveguide “leads”, whi 
individually may be of arbitrary cross section, and 
“coupling region’’ from which the waveguide lea: 
emerge. Various aspects of the theory of wavegui: 
junctions have been considered in recent years by 


number of authors, and a considerable body « 
systematic theory ° centering on the use « 


impedance, admittance and scattering matric 
has been built up. Nevertheless, for the purposs 
of this paper a formulation, which can be brief, bu 
which is in some respects new and more complet: 
is needed. Part of the formulation will depend, « 
course, on certain rather well-known general results 
of the theory of waveguides (see footnotes 2 and 7 
The domain of the electromagnetic field in 
waveguide junction will be denoted by V, the com 
plete boundary of V will be denoted by S, and the 
inward normal unit vector on S will be denoted by k. 
The surface S and also, in part, the boundary cor 
ditions to be imposed may be described (in two 
typical cases) as follows. If the domain is of infinite 














I IGURE 2 Slotted waveguide 


Illustrating Sa,Sua (n=2); Se is not shown 


extent (fig. 2), Vis bounded internally by n closed 
surfaces S,,+S,. (m=1, 2, n), where S,, is a 
transverse surface (the terminal surface) in the m™ 
waveguide and S,, S,,.” together enclose the termina- 
tion of the m™ waveguide. Although V in this case 
is externally unbounded, it is convenient to employ 
a large spherical bounding surface S.. (of radius 7 
appropriate limiting processes being implied. The 
complete boundary of V is then S=S..+S,+S8,+- 
S,+8,. On all parts of S except the terminal 
surfaces, the field is to satisfy homogeneous boundary- 
conditions: on Ss, lim (rE) is bounded, and E satis- 
fies the outward-radiation condition; on S,, the 
tangential component E, of E vanishes (EF denoting 
the electric field). If the field is confined to a finite 
domain by a perfectly conducting metal surface (fig 
3), then S=S)+S,+ - +--+ +8,, where S) coincides 
with the metal surface. In this case the homogene- 
ous boundary condition is simply E,=0 on S,. (if 
metal walls are considered finitely conducting, but 
are sufficiently thick, Sp may be taken on the outer 
surface where E, hence E,, is substantially zero 
The whole of the space and structure within V 


can be regarded as a linear, source-free medium, 
See footnote 2 
7J.C. Slater, Rev. Mod. Phys. 18, 441 (1946 
*'D. M. Kerns, J. Research NBS 42, 515 (1949 (The references cited 


believed to be the ones most useful in connection with the present paper 
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which is in general nonhomogeneous and aniso- | single frequency w, are caused to appear on the 
tropic. The conductivity and the electric and mag- | terminal surfaces. If, as is assumed, the terminal 
ef! netic inductive capacities of the medium are to be | surface S,, is sufficiently far from any discontinuity 














“al given by the real point-functions o=o(r), e=e«(r), | (or any other departure from the conditions defining 
vi u=u(r), respectively. Anisotropy is taken into | a waveguide lead), the contributions of attenuated 
id account by considering ¢, €, uw to be, in general, | waveguide modes to the field on S,, will be negligible. 
aT tensors of rank two. Thus only nonattenuated modes need be considered 
ic The interior of the m™” waveguide lead is a cylin- | in describing the field on a terminal surface.” It is 
Dy drical, source-free domain +, (of finite length) in | understood that the field of a single nonattenuated 
( which o=0, wu, «=scalar constants. +, is bounded | mode may and in general will involve both incident 
by a cylindrical surface (also of finite length) on | and emergent progressive components. Waveguide 
ee which the conduc tivity becomes infinite. The ter- | and frequency being given, the number »y,, of non- 
Os minal surface S,, is a transverse section of 7,,; S, | attenuated modes supported in the m™ waveguide 
bu is necessarily finite or zero; it is naturally assumed 
let SYMMETRY that »,, 51. In general, TeM, TM, and TE modes 
- PLANE ~ will occur among the »,, nonattenuated modes in 
sults -$3 waveguide m, but it will not be necessary to distin- 
17 | guish the several types of modes in the notation. 
n ‘ | = Of essential interest will be suitable expressions for 
om hal \ the tangential ( pny ype dew ongeerr E, Hi, of the 
the * most general E,H on §,, consistent with the above 
vk. ~ ot conditions. Let the index a (o=1,2, - +, Vm) 
CO! } , } ‘ identify the nonattenuated modes supported ¢ it. fre- 
two $ ; | quency w in waveguide m. From waveguide theory 
nits / | we know that £, on S,, can be expressed in the form 

















iis 7 
] t Evr,,): p> F AY (F,, ON Sm), (2) 
S., pel 








Figure 3. Hybrid T-junction. 


IHustrating Sar S ‘ 


is to be a closed, connected but not necessarily 
simply connected) plane surface bounded by a curve 


where r,, denotes r on S,,, the v,, are scalar coeffi- 
cients, and the ¢. are derivable from eigenfunc- 
tions of certain two-dimensional boundary-value* 
problems formulated for S,, and its perimeter C,,. 
The vector ef,, like E,, lies in the plane of S,,; 
e°,, is considered to be defined only for r on S,,. 
No coordinate-dependence is indicated for the v 


risa 


osed C,, (which may consist of one or more distinct | since S,, is considered to be in a fixed position in its 
is a parts As the figures indicate, S,, is to be located | waveguide lead. Similarly, H, on S,, can be ex- 
m™ within its waveguide lead at some distance from any | panded in the form 
Lina- discontinuity. (A portion of a waveguide lead may 
case be included bodily within V, as in fig. 2; in any case Hm 
. : = v. ( > 0 ) ( ar (3) 
et the interior of a waveguide lead from the inner H(r,,); —'" uh’, (Pm); Pm OD Sy 0 
- . . A 
portion of the junction out to the terminal surface 
The is an integral part of the domain V.) here the 7,, are scalar coefficients and the A®, are 
Harmonic time-dependence, at angular frequency | defined by 
sina] w, is assumed; we shall deal with the complex electric 
lary- and magnetic field (amplitudes) E=E(r), H=H(r) REAP m) = (FO mu) Km X Culm +) 
atis- itting the t dependent factor exp(jot). With- ' 
at —_ , ne he 7 dependen ane = Je . where k,, denotes k On Sm, Nmu IS the Wwave-admit- 
the in V, then, £,H satisfy Maxwell’s equations in the ; : 
ons wd tance of mode w in waveguide m, and ¢)=1 ohm 
ee , ; (see the following paragraph). Thee®, and the A®, 
“f VXE jou, ) , may be assumed to be real and to satisfy the ortho 
iw { : ° . . 
. ; 7 gonality and normalization relation 
rides 7 XH=(jwe+o)-E. j 
rene- * y " ‘ . . (1,A=p? ‘ 
(If Tt . » A x hd, ‘k,dS Our é , \v 
The rationalized meter-kilogram-second system of 2 J Sm 10, Ay) 
ruil 
units is assumed. ; ~ ' . 
puter To de ee a +e. : Y Equations ft and 5 each represent combinations of 
Under the specified conditions holding on S and : betes 
D ,  £ ae — what is necessary with what is convenient. In par- 
in V the electromagnetic field within V will be ae : _— — 
im | wee . ve ticular, eq 5 is automatically satisfied if modes yu,A 
determined by the boundary conditions on the term- ; a 
ium, * eda "eb a field, | are not mutually degenerate or if one is a TF and 
inal surfaces. The ultimate sources of the field, t] se cee 
rd » other a 7; 4 
which are to be found within the waveguide ter- — ' 
- minations (pictured in fig. 2, implied in other figures), Attentuated modes are excluded mainly because these higher-1 
a . actions are usually avoided in practice he finiteness of the 1 
are of interest only insofar as certain fields, of a | tons ie ufually at | putas a0 tar as the enhenenent ona , 
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Clearly, v,, and 7, are respectively linear meas- 
ures, relative to the standard terminal fields €&, and 


h°,, of the contributions of mode » to E, and H, on 
the terminal surface S,. The variables a,,, dn, 
defined by 

2A.np=Vag t fetes, 20 p= Uap — Sobmp (6) 


will be employed to a very limited extent in this 
paper. a,, and b,,, are respectively linear measures 
(relative to e&.,) of the electric fields of the incident 
and emergent progressive components of mode u 
at S,. The second of eq 6 shows that if 5,.,.=0, 
then the corresponding value of Ump/ime (called the 
characteristic impedance of mode yu in waveguide m) 
is equal to fp. So far as the present paper is con- 


cerned, { may be considered to be primarily a 
dimensional constant inserted to bring about an 
attractive dimensional scheme. Indeed, from eq 4 


and 5 one finds (considering 6,, and k,, to be dimen- 


sionless) [e2,)]—[A%,]—meters', and then from eq 
2 and 3, [v,,]—volts, [?,,,]—amperes. 
It is appropriate to review the electromagnetic 


situation in the junction as a whole. For convenience 
let v and i denote the column matrices whose ele- 
ments are respec ‘tively the Nv,,’s and the N i,,’s, 
N=y,+%+4 +v,. v uniquely determines and 
is uniquely detemntned by E, on all the terminal 
surfaces; similarly, i determines and is determined by 
H, on all the terminal surfaces. Now the speci- 
fication of either E, or H, on all terminal surfaces 
(together with the homogeneous boundary condition 
holding elsewhere on S) is just sufficient to determine 


E and H throughout the domain V. Thus, if 2 is 
given, E,H, and i are determined; if given, 


determined.” The existence of a 
relation connecting v and i 
there are (under the conditions 
linearly independent 


HE and v are 
homogeneous linear 
implied. Moreover, 
of the problem) exactly N 


Is 


electric fields possible in V; similarly, there are 
exactly N linearly independent magnetic fields 
possible in V. This fundamental property is 


expressed analytically in the following paragraph. 

We define the electric basis-field e,,, as the electric 
field '' in V corresponding to the special boundary 
condition implied by ?,,=—1,®,=0 (Amy). The 
connection, as well as the difference, between Cn 
and e8, is to be noted: e,,, is defined throughout V;, 
its tange ntial component (@,,,.), reduces to e&, on S,, 
to zero on other terminal sur i ‘es. From the bound- 
ary conditions it follows that the N e,,,’s are linearly 
independent. The Ein V corresponding to arbitrary 
v may be written 


E(r) =>) Cnu(¥) 0 m.=er, (7) 
hy 
where e is the row matrix e=(e,, - - = 
Cnr) Similarly, the magnetic basis- field hy 1S de- 


fined as the magnetic field in V corresponding to the 


10 It may be assumed, with no appreciable loss of generality, that the electro- 
magneti« 
unique, and that E=H=0 (throughout V 
in footnote 2, p. 134; also paper in footnote 8, p 


field in V corresponding to arbitrarily prescribed oe or i exists and is 
corresponds to e=0 or i=0. 


nor 
535 


See paper 
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The 


boundary condition i,j=1, ixn=0 (kx). 
in V corresponding to arbitrary 7 is then 


Hir) DS Ahal(rin hi, (es 
Dn 


where A is the row matrix h=(h,, -- - Ay - 
h,,,). The possible E’s and the possible Bs « al 
elements of linear vector spaces,” of dimension A 
E has the coordinates v relative to the basis e; H ha 
the coordinates i relative to the basis h. The pal 
ticular basis fields introduced in eq 7 and 8 are thy 
simplest ones to start out with; e, h,or both togethe: 
will accordingly be called a primitive basis. 

The usefulness of the following definition will be 
come apparent. The bracket |E', H’) of any = - tri 


field E' in V and any magnetic field H’ is 
defined by 
(e,a\—) >> | &'xH)-kadS, (9 


m=! J Om 


where, as always in this paper, the superposed bar 
denotes the complex conjugate. It is to be empha- 
sized that E' and H? are by no means necessarily 
associated electric and magnetic components of the 
same elocteetionnetia field. The most important 
algebraic property of the bracket is exemplified by 
a\E', H"| 


\(aE' + bE’), H°| +- bE? H°), 


(E', ld, 


where a, 6, c,d are any constants. (There will be 
no need to define or use brackets of the type [H, E}.) 

The orthonormalization of the standard terminal 
fields (as expressed in eq 5) has as an immediate con- 
sequence a corresponding property of the basis fields: 


(E', (cH? +d] —[E', HZ 


[Emy, Rar] Smy, try (10) 
where bap, v=! for mu=/), Bmp, td ) for mul). 
The coordinates ire ‘lative to the given ‘tele of arbi- 


trary E and Hin V may be defined by 
[E, hy] Dirs [@my, Hf] 
for if E=ev and H=hi, then, with the aid of eq 10, 


one finds that eq 11 do in fact yield v, and im,. 
Let H(e,) denote the magnetic field associated with 
€,, and let £(A,,,) denote the clectric field asso- 
ciated with A,,, (so that VX e,=— jou - A(e,) and 
VX hn,=(jwet+o)-E(h,,)). Replacing FE and H in 
eq 11 by £(A,,,) and H(e,), respectively, we write 


[Z (h,,,), Rol Zn muy 


thereby defining the N-dimensional square matrices 
Z and Y. (The first index- pair attached to the 


Possibly one might prefer to say that €m, is a field of electric type inasmuch 
as the units of em, are those of e°., and not those of 
see, @ F. D 


(11) 


bmuy 


(Emu, £1 (€,)] 7 (12) 


? For the mathematical postulates defining suc h Spaces, 
Murnaghan, Theory of group representations, p. 11 (Johns Hop kins P A ss, Balti 
more, Md, 1938) 
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matrix element labels the row in both cases.) Now 
if E and @ are associated (E=£(H), H=(B)), it is 


clear that 


E=E (D Mwir) = DDE (Nes) ine, 
Kx 


kx 


H FD) Cex ex)= >) Fh (Cx) Vex, 
ke kx 


since £ and ® are linear functions of their argu- 
ments. By inserting these forms for E and H into 
eq 11 one finds 


— ° > ss. on , 
Vir > Zur rtess Umu p a4 Pape oe (13) 
ke kx 


Thus the coordinates of associated E and H are 
related by v= Zi or, equivalently, by i=Yv. Z and 
Y are, respectively, the impedance and the ad- 
mittance matrices (relative to the primitive basis 
characteristic of the waveguide junction. 

The bracket [E, Hi) of (associated) E and H is 
(directly from eq 9) the integral of the inward 
normal component of the complex Poynting’s vector 
extended over the aggregate of the terminal surfaces. 
Thus the total (complex) power influx W across the 
terminal surface is 

W =[E, B= >> tmp l€ms, Ral taxr=0* v. (14) 


Mp, lr 


Here and subsequently the star is used to denote 
the Hermitian conjugate ( transposed complex- 
conjugate) of a matrix. The additional expressions 
W=i*Zi=v*Y*v follow immediately from eq 13. 

Consider the introduction of a new electric-field 
basis e’ =(e;{ ex.) related to the primitive basis 
e=(e, €,,,) by means of the linear trans- 
formation e’=ea, where a is a unitary matrix 
(a*=a™')." The postulate that eq 7, 14, 8 shall have 
invariant meaning determines the transformations of 
v, 7, and A relative to that of e. Namely, from 


E=ev=e'v’, (7) 
W a*vo=(2’)* 0’, (14) 
H—hi=-h’i’, (3 


and the given e’ =ea one finds easily 


e’=€a, o’=a~'g, 
(15) 
h’ ha, "‘y a a. 
And from eq 12 
Z'=a™"'Za, Y’=a"JYa. (16) 


*s 


(It should be noted that the formulas for A’, 7’, 
Z’, and Y’ are written for unitary a and do not hold 
unless @ is unitary.) 





rransformations with an arbitrary nonsingular a could be considered but 
will not be needed in this paper 


With the above provision for change of basis the 
electromagnetic formulation, as far as needed here, 
is essentially complete. It may be remarked that 
Zand Y can be regarded as metric tensors of the 
vector spaces of eq 7 and 8, respectively (the metric 
is not in general Hermitian, to be sure).“ The 
scheme acquires additional meaning when it is 
recognized that for each bracket there is an expres- 
sion involving volume integrals extended throughout 
the domain V (appendix, 1). 

Although the discussion will relate primarily to 
Z, Y, and the corresponding basis fields, the results 
to be obtained for Z and Y will hold also for the 
scattering matrix S, which is defined as follows. Let 


2a=v+fot, 2b=v—fyi, (17) 
be the matrix form of eq 6. Then the scattering 
matrix furnishes the relation 


b=WSa, (18a) 


and a simple calculation shows that 


S=(Z—¢,(Z+¢,)", (18b) 
where {) is to be interpreted as a multiple of the 
N-dimensional unit matrix. Equations 17 and 18 
will have invariant meaning under a change of 


basis provided 
a’=a™'a, b’=a~'b, S’=a "Sa. (19) 


The matrix ¢), which must transform like Z (from 
eq 17), is invariant under a unitary change of basis. 
So far as the subsequent symmetry analysis is 
concerned, Z, Y, or S is an arbitrary nonsingular 
matrix subject only to the consequences of structural 
symmetry of the waveguide junction. (It may be 
understood that Z, Y, and S are such that the real 
part of W, Re(W), can not be negative, but this 
condition is not used in the analysis.) The following 
special conditions (of electrical origin) are of interest. 
I. The nondissipative condition (Re(W)= } 
0): Z, Y are skew-Hermitian; S is uni- 
tary: Z*=—Z, Y* Y,S*=S- 


Il. The nonreactive condition (Im( W)=0): 
Z,Y,S are Hermitian: e.g., Z*=Z. 


r (20) 


Ill. The reciprocity condition: Z,Y,S are 
symmetric: e.g., Zmurc—Zar, ma, OF, iN 





the matrix notation to be used, Z=Z i 


These conditions may be incorporated, at will, after 
the main results have been obtained: of course, the 
more usual cases are III (alone) and the combination 
of III] and I. The matrix conditions in I and II are 

4 For the mathematical postulates k ading to a positive-definite Hermitian 


metric-tensor, see, e. g., Murnaghan, p. 17. To make connection with Mur- 
naghan’s notation, one may define 


(B'|E°) =(B', @(E°) (or (|) =(& (1), WH). 
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invariant under transformations of the type needed, 
but the property Z=Z is not (unless a is real)." 
In general Z’ is not simply equal to Z’ but rather, 
from eq 16, 

Z'=aZi &(aZ’ am"\& &a)Z'(&a)~'. (21) 
Obviously Y’ and S’ are subject to the same condi- 
tion when Y and S are symmetric. 


Ill. The Rotation Operator; Symmetry Groups 


The geometric concept of spatial rotation and 
reflection of a waveguide junction or of an electro- 
magnetic field is represented and in fact analytically 
defined by certain transformations of the tensor 
point-functions involved (e. g., o(F), @n,(r), Hr), 
ete These transformations are basic for what 
follows and will be set down presently. In the 
present context tensor components will be distin- 
guished by the use of the letters z, y, 2 as indices, 
and a fixed, orthogonal three-dimensional basis is 
to be understood The three unit vectors of such a 
basis plus an arbitrary origin O define a rectangular 
Cartesian coordinate-system Oxyz (r,—2). 

For the present purpose the intrinsic geometric 
properties of a particular rotation are conveniently 
characterized by a tensor R (independent of r) 
whose components ?,, are such that R=(,,) is a 
real orthogonal matrix. The determinant of R may 
be +1 or 1; a rotation (as the term ts used here for 
convenience) may be termed proper or improper 
according as the determinant of ? is +1 or —1. 

As a preliminary to the complete expression of the 
desired transformation, we consider the transforma- 
tion F(r)—>F’(r) of a tensor F (of rank 0,1, or 2) 
furnished by 


F’(r)=Fi(r (22a) 
F'(r)= >> RFP) (22b) 
y 
or 
FY (r)= SO RR FP (0) (22c) 


This may be called “local’’ rotation, inasmuch as the 
components of F’ at the point r are related directly 
to those of F at the same point. It will be conven- 
ient to use the notation r’=Rr for eq 22b when 
F(r)=r. 

It may be noted that in electromagnetic theory it 
is customary (perhaps invariably so) to consider 
electric charge a scalar, thereby determining E as a 
vector and H as a pseudovector. Since a pseudo- 
vector is a (antisymmetric) tensor of rank two, eq 22c 
applies; but in terms of the usual pseudovectorial 
components, say //, of H, eq 22¢ becomes 

Hi(r)=det(R)S> RH Mr), (23) 
. 


~ 
manifested in Z=Z7 (in the primitive basi 


“% The fact that reciprocity is is 
by hypothesis) real 


depends upon the fact that the standard terminal fields are 


where det (22) is the determinant of ?. It is well t 
note also that in the present framework EX 4H is 
vector and EX H-k is a scalar. 

By combining a local rotation and a suitab| 
related functional transformation, which may b 
called “parallel transport’’, we obtain the transforma 
tion to be termed rotation of a tensor field (or, briefly 
rotation). In the sense of this term, the rotatio 
corresponding to R applied to the tensor point 
function F produces a new tensor point-function, t 
be denoted by PF, whose components at the poin 
r’=Rr are equal respectively to those of the local] 
rotated F at the point r. Thus, for all r in th 
domain V of definition of F, 


Pak (r’)= Fr) (24a 
[PeF(r’)|-= Rp Fr (24b 
[P,F r’)|, D> RR, F (0). (2 4¢ 


izy — 


P,F is defined in the rotated domain conveniently 
denoted by PgV, PeV being such that r’ is in Pe! 
ifrisin V. It should be observed that the meaning 
of a rotation depends upon the location of the origin 
O of Oxyz (at which point F undergoes only local 
rotation). By what amounts to a change in the 
designation of the independent variable, one may 
write, say for eq 24b, 


[P,F(r)|.= Ry F(R“), (25 
v 


which holds for r in PgV. In the later applications 
of eq 24 it will in fact generally be convenient to 
have r as the argument of PpF. 

The analytical construction of eq 24 belies the 
simplicity of the underlying geometric picture: for 
proper rotations, at least, the transformation defines 
what might be described as a rigid motion (in which 
one point is fixed) of a tensor field; and for scalars 
and vectors, at least, the transformation is easily 
visualized. (Simple illustrations are given in fig. 4 
and in connection with example A; section V.) 

We impart operational meaning to the symbol Pz 
by saying that it stands for the operation (local 
rotation plus parallel transport) by which the fune 
tion P,F is produced from the function F. It is 
convenient to describe the properties of the trans- 
formations 24 by describing the properties of Pe. 

The properties of Pe to be listed now are all 
rather evident from the geometric picture corre- 
sponding to eq 24 and to some extent from eq 24 
themselves. No proofs will be needed here. (a) Px 
is linear: 

P,(cF +dG) 


cPpF +dP 2G, 


* Equation 24a essentially reproduces Wigner’s definition; eq 24b,c give the 

generalization needed here k Vigner, Gruppentheorie und ihre Anwendung 

suf die Quantenmechanik der Atomspektren, p. 113 (Vieweg and Sohn, Braun- 
Ls) 


schweig, 1931 
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where c, d are scalar constants and F, G are tensors. 
(b) Pe is distributive with respect to tensor products | 


| The equality of the first and last members, holding 
identically in r’’ for an arbitrary vector function’F, 


in that | shows in the vector case that PsP, as defined is 


P,(AFG)=(P2F)(P 2G), (26) 


| 


where on both sides the same type of tensor product 











is implied (e. g., w»H,EXH). (c) If the product 
i 
' 
1 PRE (r')= E(r) 
i] 
‘ 
‘ 
‘ 
' 
' 
i 
y ; 
, 
P r’=Rr 
/ 
‘ r 
a 
0 x 
FIGURE 4 Proper rotation of a vector. 


FG isa scalar (ce. g., F=E-0, G=E, FG=E-c- E) then 


FGdV (PF) (PeG)d V, (27) 
Jv J PeV 
The same property 


(d) The 


and in this sense P, is unitary. 
obviously holds also for surface integrals. 
identities 

V.[Pef(r)| PIV, f(r), 
V,-[PrE(r)] 


VX [PrE(r)] 


P,|V,-E(r)], (28) 
P,[V,X E(r)), 


in which the subscript r denotes differentiations 
with respect to r, lead to the statement that Pe, 
commutes with the V-operations. 

The product PsP, (of Pp followed by Ps) is de- 
fined by (PsP,2)F = Ps(P,F) for arbitrary F, where 
Ps(PpF) denotes the function produced from F by 
first applying P, to F and then Ps to PpF. The 
product Of two P-operations is itself a P-operation 
and in fact PsPp= Psp, where the matrix represent- 
ing SR is the matrix product SR. It is worth while 
to verify this multiplication law, say in the vector 
ease. Letting r’’=Sr’=SRr and thrice using eq 
24b, one may write 


(Ps[PeF P| }2= 2 SalPeF OV), 
y 


20 Sry R,. P(r) 
y.2 


Ds R),.F Ar): =|PspF (v’’)),. 


equal to Psp. Since the collection {R) of all three- 
dimensional real orthogonal matrices, upon matrix 
multiplication, constitutes a group, it is clear that 
the collection {P,} of all rotations of a given tensor 
function, upon multiplication as defined, also con- 
stitutes a group; {R} and {P»} are abstractly iden- 
tical and may be identified as the three-dimensional 
rotation-reflection group. The groups of interest in 
the present discussion are subgroups " of this group, 
as will become clear from the content of the next 
few paragraphs. 

A tensor function F will be said to be invariant 
with respect to the rotation P, (in which 0 is a fixed 
point) if 

P,F (r)=Fir), (29) 


for all r in the domain V of definition of F. (For 
this equation to hold it is necessary that PpV=V.) 
When eq 29 holds, Pg will be said to be a covering 
operation of F. 

For any choice of O, the corresponding set of 
covering operations P9,Pr,Ps,--- of a given 
function constitutes a group. For, associative mul- 
tiplication is defined within the set: (PePs)F 
P,(P.F)=P,F =F; moreover, the set contains the 
identity operation P, (obviously) and the inverse of 
each of its operations: Ps"'F==P,"'P;sF=F. 1t may 
happen, of course, that a group of covering operations 
actually consists only of the single element P,; this 
trivial case is obviously of no interest in what follows. 

It may be recalled that the structure of a wave- 
guide junction is described by three functions, o(r), 
u(r), e(r), which are in general tensors of rank two, 
defined in a region V. If o,y,e are invariant with 
respect to Pp, we say that the structure is invariant 
with respect to Pg and that P, is a covering operation 
of the structure. For any choice of O, the set of 
covering operations of a given structure will evidently 
be a group. If 0 is chosen so as to make the group 
of covering operations as large as possible, then this 
group may be said to characterize the symmetry of 
the structure and be called the symmetry group of 
the structure." The symbol {P,}, used above for 
the three-dimensional rotation-reflection group, will 
be used also for symmetry groups. 


IV. Symmetry Analysis 


Consider a waveguide junction having the sym- 
metry group {P»}, and let BE, H satisfy 
VX E=—joyu-H 
>in V, 


E,=FonS, (30) 
VX H= (jue | ¢)-E) 





’ For a discussion of these groups, see, e. g., J. Rosenthal and G. M. Murphy, 
Rev. Mod. Phys. 8, 317 (1936 

* Ify,eare scalar constantsand «=O in V, then the necessary condition Pg V= V 
contained in eq 29 becomes also sufficient and the symmetry group will be de- 
termined by V. In this way the present formulation includes cases in which 
symmetry is determined by geometrical figure alone 
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where the vector function F, defined on the complete 
b@undary S of V, denotes in a unified notation any 
prescribed set of electric fields on the terminal 
surfaces and the homogeneous boundary condition 
holding on all other parts of S. Let the whole 
system—structure and field—be subjected, in imag- 
ination, at least, to the (proper or improper) rotation 
in space corresponding to Pg. One obtains an ana- 
lytical description of the resulting situation by apply- 
ing Pp to V, to S, and to both sides of each of the 
three equations above. 


VX(PrRE)=—jo(Pru)PeH 
in PeV, 
VX (Pe) =(joPre + Peo PE) 
P,AE;) P,¥ on P,S. (31) 


These equations hold formally whether or not the 
structure is invariant with respect to P,; indeed, if 
P», is a proper rotation, the configuration of the sys- 
tem relative to itself is in no way altered. However, 
the structure is by hypothesis invariant, so that eq 
31 is equivalent to 


VX (P rE) jou-P,pH 
inV, (PrE),—P,F on S 
((P,H (jwe T o)-PE) 


(32) 


The differential equations here are the same as those 
in eq 30; that is to say, the macroscopic Maxwell’s 
equations of the problem are invariant with respect 
to Py. The field PpE, PgH appearing in eq 32, 
although closely related to the original E, W, is 
nevertheless a new field relative to the fixed axes 
Oxyz; and since the rotated structure is indistinguish- 
able from the nonrotated one, PrE, PpH may be 
counted as a new field relative to the structure. One 
may, of course, consider that only the field, and not 
the structure, is rotated.” The essential observations 
to be made here are: (a) If BE, H is a possible field 
in V, so also is P,E, PeH (so that the vector spaces 
of eq 7, 8 are closed with respect to symmetry opera- 
tions), (b) PeE, PpH is essentially determined by 
the transformed terminal fields prescribed in eq : 
since the stated boundary conditions holding on ali 
parts of S other than the terminal surfaces are not 
affected by Pp ; 

It follows from (a) above that the set of electric 
fields coinciding initially with the primitive basis e 
spans a representation of the symmetry group {P,}. 
Namely, since é,,, is a possible electric field, PrEmy 
is also a possible electric field, and so Pra, is @x- 
pressible as a linear combination of primitive electric 
fields: 


Ppem AP) > ealr) D(R)a mu. (33) 


tr 


(This equation is not to be regarded as defining a 
change of basis; the coefficient D(R), ,, is to be 
* This point of view is helpful in fixing the significance of the waveguide index 


m(m= 1,2 n), which is to be understood to identify a waveguide lead in a 
fixed spatial lo ation 
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regarded as a tensor component defined relative | 
the basis e of the N-dimensional vector space intre 
duced in section II.) We take D(R)y », as tl 
element in the (JA)™ row and the (mu)™ column « 
an NXN matrix D(R). These matrices constitut 
a representation of {P,}. For, calculation of 


Po Prem.) = >> (Pe evn) (Rin, mu 
Th 


D> Cx D( Qe nD(R)» mu 


th, Ke 


>) Cxel D(Q) DCR) ex, mus 
k 


a 
and comparison of this result with 


(Pe Penn D>) Cee (QR) xx muy 
Ke 


shows that the matrix associated with the product 
PoP» is in fact the product, 
the matrices individually associated with Pg and P, 
To determine further essential properties of 
D(R) we utilize observation (b) above and consider 
the transformation of the terminal-surface boundary 
conditions that define e,,,. Upon transformation 
by Per, the standard terminal-field e®,, defined on 
terminal surface S,, in waveguide m, goes over into 
the field P,pe’,, defined on (and tangential to) the 
terminal surface in an equivalent waveguide whose 
index may be denoted by R(m). Pre, must accord- 
ingly be expressible linearly in terms of the standard 
terminal fields in waveguide R(m). Thus, letting 
R(m)=k to simplify the typography, 
PO (Pe) = >* Gel?) DCR) ex, mus (r, on §,). 
<= 


(34) 


Since the tangential components of Pre, are zero 
on terminal surfaces other than S;, eq 34 determines 
Prem, Hence D(R)», »=0, 1#AR(m); also, since 
the standard terminal-fields are real, D(R) must be 
real. The nonvanishing elements of D(R) can be 
presented in a set of n gee es D(R). 
(D(R ey. m iS square and of dimension »,,). The arrange- 
ment of these submatrices within D\ R) i is in accord- 
ance with the scheme of permutation, m—R(m), 
of terminal fields among equivalent waveguides 
(In example A, section V, the permutation |-—2, 
2—>1, 3->3, 4-4 for n=4 occurs.) Thinking pn- 
marily of the general case v,, >1, but not excluding 
Ym==1, we determine the nature of the submatrices 
in the following manner. For the moment let two 
modes (in the same or in equivalent waveguides) 
be termed equivalent if they have equal wave- 
admittances, so that equivalent modes are not only 
“degenerate’’ but also of the same kind (TE, TM, 
or TEM). It is clear that the right-hand side of 
eq 34 can involve only modes equivalent to the u™ 
mode in the m™ waveguide. Hence (assuming 
suitable ordering of mode indices) the nonvanishing 
elements of D(R),,, will appear in smaller square 
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submatrices (‘‘steps’”’) lying on the main diagonal 
of D(R),m. Each step relates equivalent modes 
and is itself a real unitary (=real orthogonal) 
matrix. For, considering any two equivalent modes 
uw’ (distinct or not), one finds with the aid of 
eq 4, 5, 27, and 34, 





Su \ Con mu) e. ° a dS 
' e Ss m 
(Pre) -(Pré,.) dS 


Sk 


(Fo mu) 


Pan Li R),. mul WR my? 


a,« 


SOD) mul WR x. mul, 


where K=R(m), and 7,,, is the common value of the 
wave-admittances of the modes involved. Thus 
each step, hence each D(R),,,, hence D(R) itself, is 
a real unitary matrix. 

Parenthetically it may be noted that even if 
vm iS large the dimension of an individual step 
ordinarily will not exceed two or three. (It can be 
shown that the highest symmetry degeneracy of wave- 
guide modes is two-fold.) Of course, v,,=1 fre- 
quently oecurs in actual problems; the submatrix 
D(R),, then consists of a single one-dimensional 
step and is necessarily equal to +1. 

The magnetic fields of the basis A transform under 
rotation according to the same group-representation 
as do the electric fields of the basis e. Perhaps the 
most easily visualized proof of this is the following. 
From eq 34 and 4 one may obtain 
Pp [he (8) X Bm|=(SSAR i) D (Ry ex. mul X he; 


« 


(r, on S,; k=R(m)) 
where, again, the wave-admittances drop out. The 
left-hand side may be replaced by [Peh?, (rd| X k 
(ef eq 26); and since the fields are transverse, k, may 
be canceled from the resulting equation. Thus one 
obtains for the tangential component of Peh,,, on S, 
PR’, (r.) = SORE. (ty) D (Rx. mus 


« 


(r, on Sy; k=R(m)). (35) 


Since the tangential components of Ph,,, vanish on 
terminal surfaces other than S,, eq 35 implies 


Pehn,(P)= Shale) D (Ray, mu; (rin V), (36) 
Tx 


where D(R),. ».=0, 1#R(m), as above. In view of 
eq 34 and 35 it is clear that to determine ))(f) in a 
concrete case one may consider the transformation 
of either the electric or the magnetic standard termi- 
nal-fields. If the magnetic terminal-fields are con- 
sidered, it should not be forgotten that they are 
pseudovectorial. 
The transformation of an arbitrary electromag- | 
netic state upon rotation by P, is obviously deter- 
mined by eq 33, 36 and is conveniently presented in | 
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matrix form. In matrix notation eq 33 and 36 
become 

P,e=el(R), Pph=hD(R), (37) 
where P,e is interpreted simply as (Pre), . 
Prem, Pry»,); Prh, similarly. For E=ev 
and H=hi, then, 

P,E elD(R)v CUR, P,H ADR) hip, 
where, by definition, D(R)v=vp, D(R)i=irn. That 
is to say, if E has the coordinates pv relative to e as 
basis, then PE has the coordinates rz relative to the 
same basis; similarly for H. The interpretation of 
P, as a matrix operator, defined with respect to v 
and i (which represent electromagnetic states) as 
operands, is immediate: 

Ppv = D(R)v=vrp, Ppi= D(R)i= iz. (38) 

The matrices D(R) are, of course, basis-dependent. 

A unitary change of basis induces the similarity 
transformation 


D(R)—D' (R)=a™"'D(R)a, (39) 


as is readily verified with aid of eq 15 and (choosing 
one of several possibilities) eq 37. The essential 
practical problem in symmetry analysis, briefly 
stated, is to find a basis—that is, to find a transform- 
ing matrix a—such that the representation J) (of the 
symmetry group involved) whose matrices are D(R) 
will be reduced-out ” by the transformation 39. 
Such a basis, in which the matrices ))’(R) of the 
representation J’ appear in reduced-out form, will be 
termed a symmetry basis and the corresponding 
coordinates, symmetry coordinates. The principal 
items to be considered in the remainder of the 
analysis are the construction of D)’, the construction 
of a, and the physical consequences ; it is expedient to 
consider these items in this order. 

A reduced-out representation D’, in which the 
irreducible components of J) are to appear explicitly, 
can in principle be written down as soon as the ir- 
reducible components of D are determined. Let x 
denote the character of )), and let x? denote the 
character of the p™ irreducible representation D? of 
the symmetry group involved.*" Then, by a basic 
theorem, the nonnegative integer c, that tells how 
many times )? must appear in J)” is given by ” 


Cp=(1 9) Bax? (ft) x(h), (40) 


where g is the order of the symmetry group and the 


2° For the concepts and theorems of the theory of group representations that 
will be needed here it will be convenient to refer to Wigner (Gruppentheorie, 
cited in footnote 16), especially chapters [IX and XII 

2! Character tables are given, e. g., by Rosenthal and Murphy (footnote 17), 
and by G. Herzberg, Infrared and Raman spectra of polyatomic molecules (D. 
Van Nostrand Co., New York, 1945 Irreducible representations have appar- 
ently not been tabulated, but many of them are easily found (for one-dimensional 
representations, Dre= x?) 

2 Wigner, p. 95 











summation goes over the group.” (The complex 
conjugate taken in the right-hand side of eq 40, 
which makes the formula correct in general, is of 
no significance here because c, is real in general and 
x is real in the present problem. Let the (distinct) 
nonequivalent irreducible represeatations actually 
contained in D (i. e., those for which ¢, #0) be num- 
bered from 1 to t. Further, let D?(R),, denote the 
element in the 27 row and the p™ column of the 
matrix D)?(R) of the irreducible representation D?, 
which must be unitary and is considered known; and 
let L, denote the dimension of D?. For the element 


in the (pra) row and the (qpb)™ column of the 
matrix ))’(R) we write 

D’(R) pre. ¢o0=SpehunD”(R) en (41) 
where the indices are to have the ranges 
a=1,2 Cy; r=1,2 oe e p=1,2 t; 
6=1,2 Ce; e=1,2 ~ * q=1,2 t 


(Whenever /, #1,, undefined symbols D?(R),, multi- 
plied by zero occur in the right-hand side of eq 41; 
such “products” are defined to be zero. The ele- 
ments of ))’(R) furnished by eq 41 may, if desired, 
be arranged in such a way that ))’(R) takes the form 


(DR)... 0 Ss ves 
D’'(R) " _— af f) ~ —_ 9 (42) 
0 See 0 DR)... 0 








in which the matrix //'(R) appears c, times, 1?(R) 
appears ¢, times, Any consistent scheme of 
ordering elements into matrices may be assumed, 
and it may be noted that an order other than one 
that leads to eq 42 will be convenient later.) 

In the formation of D’ by means of eq 41, the 
choice of the particular irreducible representation 
D”® among the unitary representations equivalent 
to D? is merely a matter of convenience. (This 
question does not arise when /,=1, of course.) How- 
ever, an irreducible representation appearing more 
than once in D’ is to be presented in identical form 
each time it appears, as indeed is insured by eq 41. 
The arbitrariness in J” implies an arbitrariness in 
the choice of a symmetry basis; of more interest in 
this paper, however, is further arbitrariness in a 
symmetry basis that remains after D’ is assumed. 
This arbitrariness will appear in the course of the 
construction of the transforming matrix a@ and will 
be utilized in a brief discussion at the end of this 
section 

Let the row-matrices e’, h’ present the basis 
fields of a symmetry basis in which the matrices of 
D’ are in fact given by eq 41. In this basis eq 37 
become P,e’ e’D’ (R), P,h’ h’D’ (R); where, 


Only finite symmetry groups will be considered explicitly. 


etc 


in accordance with the index notation establishe: 
in eq 41,e’=(e' yea), A’ =(R' yea). The characteristi 
transformation equations of basis fields in a sym 
metry basis follow directly from the form of D’(R 


, 


(43a, b 


lp ly 
Pr; va De) na? (R) ox, AR}, v0 SO) aD? (R 
l 


ep p=l 


The /, members of a set of basis fields (electric or 
magnetic) identified by fixed p and a transform 
among themselves according to the irreducible rep 
resentation D? of {P,»}, and they are said to bi 
partners and to be of symmetry species pP (relative 
to {Pz}). There are ¢, sets of basis fields (electri 
or magnetic) of the species p, in accordance with 
the range 1,2, - ,c, of the index a. More specifi- 
cally, e)., or Ai., is said to belong to the «™ 
of the irreducible representation D?; there are, again, 
c, such electric or magnetic basis fields. The same 
terminology may be applied, of course, to other 
fields or other entities that transform under rotation 
according to an irreducible representation as in eq 
45. 


row 


Let a and _ an which are to obey eq 45, be 
given by 
CO, cc™ > Cugtinn pve Clipse, 


my 


where a,rq denotes the (pxa)™ column of the unitary 
transforming matrix a@ (cf eq 15 and the parentheti- 
cal remark following eq 16). (Note that the column 
matrix a@»rq presents the coordinates, relative to the 
primitive basis, of e’,, and Aj...) From either eq 
43a or 43b we may derive an equation similar to eq 
13 but applying to the ay... Namely, if in eq 43a 
we replace e),,, by Capra (and e),, by @ap..) and then 
invoke the matrix interpretation of Pz (eq 37, 38) 
and abstract e, we obtain 


Patty: SS ap pal)?(R) oe. (44 


i 


Thus arc, as well as e),, and Al.,, may be said to 
belong to the x row of the irreducible representa- 
tion D?. Conversely, it may be seen that if the 
columns of a unitary matrix obey eq 44, then the 
corresponding basis fields will obey eq 43. (Equa- 
tion 44 is, moreover, equivalent to eq 39 and may 
be derived from that equation by equating the 


(pra)™ columns of the matrix products on the two 
sides of D(Rja=aD'(R)—D'(R) being given by 
eq 41.) 

A systematic procedure for the construction of @ 
may be based upon eq 44, assuming that the 


matrices of the irreducible representations involved 
are known. (In a concrete problem various expedi- 
ents, including judicious guessing, may often be 
employed to advantage.) We shall give the con- 
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struction * almost in the form of a recipe, deferring 
the justification to the next paragraph. Define the 


° Pp . ° 
NX1 matrix 7“) as a function of an arbitrary 
(NX1 matrix) ec by means of the equation 


v'*)—(1,/g) SDR). Pee, (45a) 
PR 


which may be written 
p\*) ((l, g)>>D?(R),.-D(R)v Gin) y, (45b) 
P 


thereby defining the NXN “generating” matrix 


} 


G‘**). Corresponding to any set of N linearly in- 
dependent v’s, eq 45 yields a set of exactly c, line- 
: (P)s (?) . 
arly independent v'*’’s; each v'*’, hence any linear 

. . (?), ' 
combination of v\*’’s, belongs to the x™ row of D?. 
Select a set of ¢, such combinations that are mutu- 


ally orthogonal and normalized to unity (in the 

Hermitian sense): these may be taken as a,-), a@p-r2, 
. . y( ? 

"+ % Qpee,- Define the matrix GY*) and construct 

the l, 


1 partners of a, by means of 


A poa= (lp/g) 55D? (RB), 2P rapa; (46a) 
- 


waTTas. ’ p \ . 
[ a g) >> D*(R),.D(R) ap xa= Go | tp zes (46b) 
R 


for p=1,---, r—1, r+1,---, 1, and a=1, 2,--.-, 
c The a@),<’s so obtained will automatically satisfy 
Og Ayer The procedure is to be carried out 
for each symmetry species (p=1, 2, -, t); the 
orthogonality a7,,a.=0, pq, will automatically 
hold. Thus one obtains altogther a complete set of 
column matrices that satisfy eq 44 and are properly 
orthonormalized to serve as columns of a. 

For the most part, the above construction repre- 
sents an immediate application of theorems given by 
Wigner.” Beyond this it is only necessary to show 


5 p0 Dan: 


) 


that (a) G*"’ is in fact of rank ¢, and (b) ay. and 
Gp, & given by eq 46, are orthogonal when a,,,. and 
ax» are. As for (a): since the representations D 
and J” are unitary and equivalent, there is no 
question of the existence of a unitary transformation 
connecting them. Hence the equation defining 
i? 
” a may be transformed formally into 


b 
{ 


Pe,’ dae ‘a (l, g>> D*(R),.D'(R). 
R 


Considering D’(?) in the form 42 (say) and apply- 
ing the orthogonality theorem given below (eq 50) 
; 7 + oe ; : 

one finds that a~'G@**’ a is a diagonal matrix having 

‘The method to be given is essentially equivalent to one given (for the con 
struction of molecular symmetry coordinates) by J. R. Nielson and L. H 
Berry man, J. Chem. Phys. 17, 659 (1949 

* Chap XII, eq 1, 3, 3a, 6, and 8. Compare eq 44, 45a, 46a with eq 1, 6, 3a, 
respectively The reinterpretation of Wigner’s equations with his functional 
P» and f replaced by our matrix Px and vr is immediate 


exactly ¢c, nonvanishing diagonal elements (each 


vee) Awe) s 
equal to 1), so that a“'@ a, hence G _is of rank 
c,. As for (b): from the relation 


(?,) 
(G°""’)*=(1,/9)>> D?(R),. D(R)*= 
R 


—— foal 


(l,/9) >> D*(R-5,, D(R- =a", 
7 


and the observation that, for a set of partners, eq 46 
holds for all x, (in their proper range) ,” one obtains 
* « (avr) * gy\eo) r 

Arcee Ape (GF ) App— A524 Op pdb= Ay zeAped, 
which establishes (b). 

Structural symmetry of a waveguide junction 
places certain restrictions on the form of the matrices 
Z, Y, and S (defined in eq 12, 18): some matrix 
elements may be forced to vanish, and the number 
of independent elements may be reduced materially,” 
as is well known in a variety of cases. One form of 
the governing equations may be found in the follow- 
ing way. The relation v= Zi holds, of course, for any 
given electromagnetic state in the junction. If the 
fields are subjected to a symmetry operation, then 
v—> Vp, i->ipz (eq 38) and vg= Zig must hold, for vg,ig 
are coordinates of a possible state in the given 
invariant structure (cf footnote 19). Now ve=Zir 
is the same as D(R)v=ZD(R)i, v may be replaced by 
Zi, and 7 is arbitrary. Hence 

D(R)Z=ZD(R), every R of {Pe}, 
und we say that Z commutes with the representation 
D of the symmetry group. It immediately follows 
that Z~'(=Y) commutes with D and that (Z7—f£)) 
and (Z+9)"', hence (Z—f»)(Z+%)7'(=S), com- 
mute with ). (In fact any rational matrix function 
of Z commutes with )).) Thus, if \ denotes Z, Y, 
or S, 

D(R)M=MD(R), every R of {Pe}. (47) 
This equation contains the conditions imposed by 
symmetry upon the several matrices characterizing 
the junction. ‘To find the consequences of eq 47 we 
consider it presented in a symmetry basis. 

Upon transformation to asymmetry basis, 1/7 and 
D undergo one and the same similarity transforma- 
tion (eq 15, 19, 39), and eq 47 becomes D’(R)M’ 


M'D'(f). Thismay be written M’= D’(R)* MDL), 
and for the element MM) ce. et of MM’ one obtains, 
using eq 41, 

Mi? oe, b= > DP" (Bee "(Rye M jcc, rv (48) 


o.T 


* Wigne:, Chap XII, eq 3a 
7 Results more detailed than those to be obtaine« 


> 


appendix, 2 
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The g equations of this type (one for each R of { Pe 
may be added to obtain 


gM cc. > 52” R),.D* Bop |M ira op. (49) 
or R 


Now according to the orthogonality theorem for uni- 
tary irreducible representations, * 


2h or D* (RB), = 8 pq Ser Sep (Q/lp) (50) 
so that eq 49 becomes 

M) ec. go, unless q=p and p *) 

, " , 51 

| epee 1/1,) >> M' oa, poo = M?,, \ 
where m?”, is a constant inde ‘pe indent of T. M’ 
determined by the values of the ¢?+ ¢ ---+¢ 
constants m?,, which are arbitrary so ler as sym- 


metry of the junction is concerned. A convenient 
step-matrix form for M’ is obtained by arranging 
elements in dictionary-like order according to the 
values of the indices p,7,a in the sequence pra. 
(The sequence paz yields ))’(2) in the form 42.) If 
t for example, and /;=2, ¢;=2, l=2, e.=1, then 
M’ vill have the particular step form 


Sra sab spa spb tra tpa 
Sra Miva m>, 0 0 0 0 
sab mi, m 0 0 0 0 
spa 0 0 m*, m, 0 0 
spb 0 0 mM, mi, 0 0 
tra 0 0 0 0 mM, 0 
tpa 0 0 0 0 0 mi, 


The general M’ will have, for each p, /, identical 
steps of dimension c,. (If ¢,—1, every p, the steps 
are elements and MW’ is diagonal.) 

Additional restrictions on M’ corresponding to the 
conditions discussed in section II, eq 20, 21, are 
readily imposed. ‘To fulfill I or 11, 1 must satisfy 
one of the conditions \/* M, M*=M, 
M*=M~'; and M’ must satisfy the same condition. 
The reciprocity condition (III) requires M=M, and 
M’ must satisfy 

M’ (aa) M’ (Ga) (52) 
—or simply M’= M’ when a is real. 

The form of M may obviously be found from that 

of M’ by calculating M=aM’'a™'. 


* Wigner, p. 91 


It is of interest to consider briefly the special ca: 
in which M is a normal ™ matrix. Either of the cor 
ditions I and II is sufficient to insure that M | 
normal.” When is normal, M can be reduced t 
diagonal form by transformation to a suitably chos« 
symmetry basis. This statement is easily estal 
lished by considering the transformation needed 


diagonalize M’. If M is normal, so also is .M’; fro: 
the form of M’ (eq 51) it is evident that MM’ can b 
diagonalized by transformation with a unitary 


matrix 8 of the form 


B yxa.q= 0, unless p=q and r=p, 


B pve, pro = B2,, 


where 8? is a suitable unitary matrix of dimension ¢ 

(The argument is trivial if c¢,=1, every p.) Th 
change of basis corresponding to this transformation 
may be denoted e’—>e’’, where e’’ =e’8 =e(a8), and 
we have 

M’—M”" =8-'M’B=(aB)—'M(a8), 

where the diagonal matrix MM” presents the eigen 
values of VM, and the columns of the combined trans- 
forming matrix a8 are eigenvectors of VM. Let us 


put a8B=y and examine the (pra)” column of 
From the form of 6, 
Y pra > fatines; 53 


} 


that is, the eigenvector y,,. is a linear combination 
of the columns of a that belong to the x” row of the 
irreducible representation D?. Equation 53 is in 
fact an expression of the indeterminacy that appeared 
in the construction of @ (following eq 45). Hence 
is, as it were, a possible a, and the basis e” is a 
symmetry basis. 

The degeneracies (of the eigenvalues of /) that 
arise in consequence of symmetry may be made 
wholly explicit by a consideration of the formulas at 
hand. The eigenvector y,.. of M belongs to the 
eigenvalue M’’,.9, x2 presented in M’’; according 
to eq 51 (applied to M’’) M’’y.a,pra=d2, Say, 
where \? is independent of x. Hence the /, eigen- 
vectors Ypra for r=1,2, .. ., J, all belong to the 
eigenvalue \?. These eigenvectors are certainly 
linearly independent (being columns of a unitary 
matrix), and so the eigenvalue \? is at least /,-fold 
degenerate. This degeneracy is necessitated by the 
symmetry of the waveguide junction and may be 
termed symmetry degeneracy. Degeneracy higher 
than that necessitated by structural symmetry, 
when it occurs, may be termed accidental, 
customary in mathematically similar circumstances. 


as is 


* A matrix that satisfies 4 M*= M* M is said to be a normal matrix; MM*= 
M*M is a necessary and sufficient condition that M be reducible to diagonal 
form by a similarity transformation with a unitary matrix. See e. g., Murna 
ghan, Theory of group representations, p. 26 

*® Although M may be and indeed usually will be symmetric (reciprocity), 
this condition is neither necessary nor sufficient to insure that M be normal 


278 








hal 


ron 
L b 
ary 


lc, 

The 
Lion 
and 


tion 
the 
; in 
ired 
ey 
is a 


that 
race 
S al 
the 
ling 
Say, 
ren- 
the 
inly 
tary 
fold 
the 

be 
rher 
trv, 
s is 


ices. 


[M* = 
agonal 
[urna- 


ity), 
al 


V. Examples 


In order to illustrate some of the text material by | 
means of examples, we shall consider the waveguide | 
junctions shown in figures 3, 5, and 6. The para- | 
meters o, wu, € In the interior of the junctions may | 
be the most general consistent with the assumed | 
symmetry in the respective cases. Although the | 
figures indicate specific external geometries, it will 
be recognized that the exact manner in which the | 
waveguides join is of no importance here, provided 
that this also is consistent with the assumed sym- 
metries. Lowest-mode operation is assumed in all | 
waveguides (v,=1, m=1,2, . . ., n); the heavy ar- 
rows in the figures, considered as localized vectors, 
suffice to characterize the electric terminal-fields for 
the lowest mode in rectangular waveguide. In 
applying the notation of the text, unneeded indices 
will be dropped ; for example, e&,, €,,, will be written 
simply as e°, e€,, (since v»,=1). 

Example A. The junction shown in figure 3 has 
the single symmetry plane and the external geometry 
of a junction widely known as a “magic T’’. Let 
the indicated symmetry plane be z=0 of Oxyz. 
The symmetry group Pe, consists of the identity 
P, and the reflection Ps in the plane; the correspond- | 
ing matrices /?}=(/?,,) are 


] 0 0 ] 0 0 
] 0 1 ol, Ss 0 I 0 

0 0 l 0 0 1 
(Det (S) 1; Ps is a particular improper rota- | 


tion.) An inspection of the figure makes it clear 
that 

P,¢=4, P,e=e, Ps& e%, P,e& =e}, 
(cf eq 34); hence 

Ps,e, =e», P.,e,.=e,, P;e,=—e;, Ps@,=,, 


+ adgg y 3), and so the N & N (=44) representa- 


tion D is 
l 0 0 0 0 l 0 0 
0 ] 0 OVO l 0 0 O 
DT) ; TDS 
0 OU ] 0 0 O —] 0 
0 O O l 0 O 0 | 


(DU) is always the N-dimensional unit matrix, of 
course.) J) is real and unitary, as it must be. The 
character x of D is x(J)=4, x(S)=0. 


P,} has two irreducible representations; both 
are one-dimensional: 
D?: D?(IT)=(1), D?(S)=(1). 


Dy: D(I)=(1), D“S)=(—1). 


(The 1X1 matrices are unitary.) Since D? and D* 


| are one-dimensional, y?=D? and x*‘=D*. From 
eq 40, 
c,=(1/2)(1-4+1-0)=2, ¢,=(1/2)(1-4+(—1)-0)=2, 


and a reduced-out representation D’ (cf eq 41, 42 
is given by 


fl 0 O O07) (1 0 0 O) 

S + & 0 ] 0 O 
D'()= , D'S) 

0 0 1 0 0 O-—!1 0 

10 0 0 1 10 0 O- 1] 














For the symmetry basis e’, h’ we write 


e’=(e;,e;,e,, &., (h’=h}, h}, h;, hi.) 

the basis fields must transforms according to eq 43 
Thus, for example, 

Pre, = D*(Ne,.=(le;., Pse.= 


D?*(S)e’,,=(1)e., 


Pye,,= D*{Nhe,,.=(1)e,., Pse.= D(S)e,,=(—1)e,. 
The same equations hold with a replaced by 6, and 
h’ transforms in the same way as does e’. Fields 
of species p and q may, respectively, be termed 
symmetric and antisymmetric (or even and odd) 
with respect to reflection in the plane z=0. Now 
e, is already symmetric and e, is already antisym- 
metric; one may choose e),—e€,, e),=e,. Suitable 
linear combinations to form one further electric 
basis-field of each species may be found immediately 
by inspection in this problem; a simple choice is 








pp = (e, +-e,)/ 2, ef, =(e,—e,)/¥2. The transforma- 
tion e’=€a so determined is 
0 1/y¥2 0. 1/¥2) 
0 1/v2 0 —1/72 
(e,,e,, e., &.,) =(e, €, e; e,) 
0 0 1 0 
U1 0 0 oJ 


and we see that a@ is unitary, as is required, and 
also that a is real, as is convenient (but not always 
possible). (One may verify D’(S)=a~'D(S)a for 
the above matrices.) From eq 51, MJ’ must be of 
the form 


fm? m?, 0 0) 
me? m?, 0 0 
M’ 
0 0 m:, més, 
0 0 mij, mi}, J 
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Upon calculating M=aM’a™ one finds, in partic- 
ular, M/,,=—M,,=0, as may well be expected even 
under the general conditions permitted here. 
Indeed, these zeros can be predicted as soon as it is 
recognized that e, and’ e, are of distinct symmetry 


species. 


The fields e,(=e),) and A,(=h!,) are both of 
species p (even); Pd jn + e,(z,y,2) = Pse,(z,y,2), 
h,(2,y,2)=Pshi(2,y,2 The z,y,z components of 
the ‘se two equations are 
€g(2, Y, 2 €42(—2, Y, 2 h(x, y, 2) = hag! »Y,2 
€y (2, ¥,2)=ea(—Z, y,2) h(x, y,2) hy(—z, y,2), 
Cge(Z, Y, Z) = Cae(—Z, Y, 2), h,.(x, Y, 2) h,(—z, y, 2), 


where for e, eq 24b is used, and for A, eq 24c¢ and 23 
are used (cf also eq 25). The even symmetry with 
respect to the plane z forces ¢g=0 and hy 
hy =0 on the plane. (For e, and Ay, which are odd, 
€y =e. =0 and hy =0 on z=0.) Despite the differ- 
ences in behavior among their components, the tensor 
entities e,, A, are unambiguously classified as of the 
same species by eq 435. 

Example B. For the junction shown in figure 5, 
the symmetry group {2} is of order g=3 and con- 
sists of the identity ?,, a counter-clockwise rotation 


P-, of 120° around the indicated axis, and the inverse 
Pe} of Pe, It is evident from the figure that 
Pee e», Pos e,, etc., and one finds x(/)=3, 
x(C- x(Cy')=0. The characters of the irreduc- 


ible representations of the symmetry group (as well 
as the representations themselves) are furnished in 
the table 


Cei..S ex’ See 
x? l | l 
x! l u u 
x’: ] Ww Ww, 
where w=exp(27j/3). From eq 40, c,=¢e,=¢c,=1: 


in this case there will be exactly one (electric or mag- 





Three-arm junction possessing a threefold axis but 
ae . 


Fiacure 5. 
no plane of symmetry 
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netic) basis field of each of the three possible sy: 
metry species. One may easily verify that t! 
fields 

e, =(e,+-e,+e,)/V3, 


C= (C1 + We2 + We;)/ V9, 


e. =(e,+ we,+ Be,)/ 3, 
and in fact constitute a symmetry basi 


obey eq 43 
aa, and the necessary form of M’. 


We list a, 


+ 
] l l l 0 0 
l - 
a =i ] w wi aa 0 0 l 
Vv ) 
l w Ww 0 l 0 
m? 0 0 
M’=| 0 mt 0 
0 0 m 


It happens that 7’ is necessarily diagonal, so that 
is certainly a normal matrix. In this example the 
reciprocity condition (eq 52) forces the degeneracy 
m’=m'. This degeneracy is technically accidental 
it might well be termed reciprocity degeneracy. 
Example C. The discussion of the waveguide 
Wheatstone bridge (fig. 6) limited mainly to a 
sketch of results. The symmetry group of this junc- 
tion is identical with that of a regular tetrahedron 
and is of order g=24. Of the five nonequivalent 
the two 
are con- 


Is 


irreducible representations of this group, 
say) 


three-dimensional ones (D? and D*, 











Ficure 6. Waveguide Wheatstone bridge. 
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tained in D’. There are accordingly three fields 
(partners) of species p and three of species q in the 
symmetry basis; it can be shown that a symmetry 
basis is furnished by® 


fe) (1 l 0 1 —1 0} fe) 
e’,, 1 0 1 —1 0 l e, 
e', 1} 0 l l 0 1 -1 e; 
Ge 2) 1 l 0 l 1 0 e, 
e\, I 0-1 —-!1 0 —!1 e; 
et. -4 LO 1 —1 0 l 1) LesJ 


























FIGURE 7, Equi alent network. 


Sign convention: m’ is the positive terminal for r,, at the terminal-pair mm’ 
positive im is into the network at m 


(Note that this matrix equation is the transpose 
&’=a@ of the equation e’=ea.) Here, instead of 
considering impedance, admittance, and scattering 
matrices in common, we shall consider the admittance 
matrices particularly. 
the form 


fy? 0 0 0 0 OY) 


0 y” 0 0 0 0 








The eigen-admittances y’, y* are in general complex; 
the condition Re(W)S0 (mentioned on p. 271) 
requires Re(y?)50, Re(y*) 50. A calculation of 
Y=aY’a , yields 


| This transformation represents an adaption of information contained in the 
tables on p. 638 of the reference in footnote 4 


From eq 51, Y’ must be of | 








4 r uM M 0 eet p ) 
rm rm Me 0 im 
: M M A - M 0 
} 
0 ; —- Me M 
M 0 KM M » M 
= “2 0 fs =, J 
where A=(y?+y")/2, w=(y’—y*)/4. The matrix 


Y is both normal and symmetric in virtue of the 
symmetry of the junction. 

The six-terminal-pair network shown in figure 7 
is a Wheatstone-bridge scheme of connections mod- 
ified by the presence of the “parallel”? and “series’’ 
elements Y, and Y,. If Y,=A+2u”=y’ and 
4u=y'—y’, then (with the scheme of signs 
noted under the figure) the admittance matrix of 
the network is the same as that of the junction, and 
the network is an “equivalent’’ network for the 
junction. 

The particular values of y’ and y* for a given 
structure depend upon the (common) distance of 
the terminal surfaces from the center of the junction. 
If the structure is nondissipative, it is always possi- 
ble, and sometimes convenient, to assume terminal 
surfaces so located that y?=0 or, alternatively, 
y*—>@. In the first alternative the equivalent net- 
work is simplified by the absence of Y,; in the 
second, by the virtual absence of Y,()’, being re- 
placeable by a perfect conductor). (This second 
case is the one presented in the reference in foot- 
note 4.) 


The author is greatly indebted to K. F. Herzfeld 
for valuable discussions and guidance in the course 
of this work. The writer is also very grateful to 
C. C. J. Roothaan for helpful discussions, and to 
H. Lyons for his wholehearted support of this work. 


VI. Appendix 
1. Transformation of Equation 9 


The integral of (1/2)(E'* H’)-k taken over the complete 
boundary of V (assuming an infinite domain) is equal to 


[E', H?]+- (1/2 (E'XIP)-kdS 1/2 [_v. EX )dv, 


Ja . 


where in the left-hand side the definition of the bracket (eq 9 
and the boundary condition on the surfaces S)(p. 268) are 
used, and the right-hand side is given by the divergence 
theorem of vector analysis. From the above equation the 
desired express on, 


- 


(E!, He) n/2) | E..F ds 4 


Sa 


9 1. ya , 
1/2 \F o-E'dl ) 
: 54) 
jw/2 | Peet Ee E)dV, ) 


may be obtained with the aid of further vector identity, 
Maxwell’s equations (eq 1), and the following relations hold- 
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ing (for both electromagnetic fieids 
k-E°=0. where »=ye/u and vacuum values of ¢,4 may be 
assumed. Equation 54 applies also in the case of a finite 
domain provided merely that the S..-integral is omitted. 


on So:H?=nE? Xk and 


2. Additional Symmetry Properties 


The complex conjugate of the element Y’ pra,q.5 of the 
admittance matrix Y’ is, from eq 12, equal to [ebrs, (Eos) ], 
so that, making appropriate substitutions in eq 54 and taking 
the complex conjugate, one obtains 


»’ dS 1/2 Obra'oe 4 dV ) 
1 


. 


(Oo) 


- 


T Ye 2) Goce € e',..0dV je 2 H) C's b) wHerav.) 
\ Jt 


We wish to show here that an equation of the form 51, which 
holds for the elements of Y’, holds also for each of the four 
integrals making up the right-hand side of eq 55. It will 
suffice to consider the third and the fourth of these integrals. 
Since the structure is invariant with respect to Pr, and since 
Pr is unitary in the sense of eq 27, one may write 


| Brew dV = [ (PrPipea)-€-(Pre’,s)4V, (54 


|, 21e'sos)-w MlepeddV= | [Pat (e'sps) ul PaBleoea) AV 


(57 


The right-hand side of eq 56 may immediately be expand 
with the aid of eq 43a: 


[) Frewee'yndV SDR ..D(R) ey [ BheaeesndV. (58 
V Vv 


o.T e 


The right-hand side of eq 57 may be expanded similarly 

soon as it is observed that Pr (ebea)—Al(Preiec) (as 

directly implied by eq 32) and that ® is a linear function « 
its argument. Equation 58 is of the form 48, and it is thu 
clear that an equation of this form holds for each of th 
integrals involved. The derivation in the text leading t 
eq 51 obviously applies. An equation of the type 55 may 
of course be written for 7 $e4,9.5, and a similar argument 
leads to the same results for the integrals making up 7! 


pTay ap 


WaAsHINGTON, June 7, 1950. 
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Comparison of Viscosities of Rubbers from the McKee 
Worker-Consistometer and from the Mooney Viscometer ' 


A. B. Bestul, G. E. Decker, and H. S. White 


Viscosities of three GR-S rubbers and three GR-I rubbers were measured at 100° C 
with the McKee worker-consistometer, a piston-type capillary viscometer, and with the 
standard Mooney viscometer, a rotating disk instrument. Rates of shear from several 


(58 


ly a : ° P 
as i tenths to several hundred sec—' were covered with the consistometer and from 0.1 sec— 
on < to several sec—' with the Mooney viscometer. The procedure for preparing the samples 


thus and taking the measurements with the Mooney viscometer was selected as giving the best 
the possible simulation of the conditions of the standard Mooney measurement, consistent with 
gz to a valid indication of the temperature of the sample and as nearly as possible a like state of 
may degradation at all rates of shear. The viscosities from the consistometer were calculated, 
ment using the Weissenberg-Rabinowitsch differentiation method for reducing observed flow 
data to fundamental quantities. The agreement of the results from the two instruments in 
— their common range varies for the different rubbers, being better for the GR-S’s. The best 
agreement is within about 2 percent, the uncertainty of the measurements, and the worst 
is within about 20 percent. The results from both instruments for X-518 GR-S agree with 
results obtained by Treloar for a “Standard GR-S” with a Piper and Scott biconical rotor 
modification of the Mooney viscometer. As the disagreements shown are not the same for 
all the rubbers they cannot be the result of constant instrumental differences, but must 
involve factors such as hydrostatic pressure, elasticity, and slippage at instrument surfaces. 


1. Introduction TABLE 1. Rubbers inve stigate d 


The McKee worker-consistometer [1],?> which was | 
recently adapted to measurements of the viscosities Rubber Origin 
of raw rubbers [2, 3], is a piston-type capillary vis- ee 
cometer. Probably the most commonly used instru- | te 
ment for measuring such viscosities is the standard X-518 GR-S Standard reference lot 50 49 
Mooney viscometer [4], which is based on a rotating x48 ORS ae 5 64 
disk. The results obtained with the two instru- GRi-s Bais. Preduction sample. os 
ments for a given rubber should agree if, for both GR-L-18 D do SI 
instruments, the method used for interpreting the 
observed data and for reducing them to fundamental 
quantities is truly applicable to the situation con- 
cerned. In order to check the applicability of the | For this reason it was felt that the standard Mooney 
methods used it is necessary to compare results from | viscometer was satisfactory for the present: com- 
the two instruments. Such comparison has been | parisons, and that it was not necessary to use the 
made with the results for three GR-S and three | biconical instrument, which has a more nearly 
GR-I rubbers at 100° C. uniform rate of shear throughout the flow field. 


Polymeri- | Standard 
zation tem- Mooney 
perature reading * 


* 100° C, 2rpm, large rotor, 4-min readings, GR-S samples prepared according 
to Taylor [5], GR-I’'s tested as received. 


II. Materials and Apparatus 


III. Measurements and Treatment of Data 
The six rubbers for which viscosities were meas- 


ured are listed, with descriptive data, in table 1. 1. Worker-Consistometer 


The worker-consistometer used and its operation 
are described in reference [2, 3]. Rates of shear 
lower than those mentioned in reference [2, 3] were 
obtained by increasing the number of gears between 
the driving motor and worm gear. 

The standard Mooney viscometer and its opera- 
tion are described in reference [4]. However speeds 
other than the standard 2 rpm are reported here. 
In unreported measurements on the GR-S rubbers 
with this apparatus and with a biconical rotor 
Mooney viscometer [6], these two instruments gave 
identical viscosities for a given rubber if the samples 
were treated similarly in the two instruments. 

The work reported here was sponsored by the Reconstruction Finance 


‘orporation, Office of Rubber Reserve 
? Figures in brackets indicate the literature references at the end of this paper 


Measurements were made in the worker-con- 
sistometer with a 50-hole shearing element at 
reciprocal rates of piston travel from about 13.6 
hr/in. to 1 min/in. These rates correspond to rates 
of shear from about 0.3 to 500 sec™'. The pro- 
cedure in these measurements is as follows: A fresh 
sample of rubber (not mill-massed) is placed in the 
instrument and forced through the shearing element 
10 times at a nominal rate of shear of 100 sec™'. 
During these passes the apparent viscosity of the 
rubber increases. Previous measurements [2] show 
that the viscosity of the rubber at a given rate of 
shear remains nearly constant for the second 10 
passes. The data at all rates reported are obtained 
between the tenth and the twentieth passes on one 
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filling of the instrument. In this instrument the 
differences between results for GR-S samples mill- 
massed for 25 passes and samples that have not been 
mill-massed are small enough to be ignored. 

The observed pressure drop across the capillary 
(P) and volume rate of flow (Q) are converted to 
shearing stress at the capillary wall (s,.) and nominal 
rate of shear (1)) by the usual relations 

s,=RP/2L, and D=4Q/rR’, (1 
where R and LZ are the radius (0.0188 em) and length 
(0.635 em), respectively, of the capillary. The 
actual rate of shear at the capillary wall (y,) is 
obtained from the relation 

Y= (D/4) (3+d log D/d log s,, : (2 
This relation is the result of the differentiation 
method for reducing observed data to fundamental 
quantities This method is described by Mooney 
[7] and others The term d log D/d log 8» 1s obtained 
for a given value of D by plotting log D against log 
the range available and measuring the 
tangential slope at the point concerned. In this 
paper viscosity (») is defined as the ratio of shearing 
stress to rate of shear, regardless of whether the 
flow is Newtonian. According to this usage the 
quantities called viscosities are not material con- 
stants for non-Newtonian substances, but are 
merely helpful ratios, the numerical values of which 
vary with the rate of shear or the shearing stress. 
This usage has been adopted by Burgers and Scott- 
Blair [8] and by Reiner [9], and has been used by 
numerous other authors. The ratio s,/y, is thus 
the v iscosity corresponding to the rate of shear V~- 


&y over 


2. Mooney Viscometer 
When a fresh sample of any of the GR-S rubbers 
used (mill-massed according to Taylor [5]) is placed 
in a Mooney viscometer and the instrument started, 


the viscosity first rises to a maximum and then 
decreases continually with further shearing. This 
decrease results from molecular degradation. Vis- 
cosity results varying by as much as a factor of 2 
can be obtained, depending on the procedure used 
in taking the measurements [10, 11]. The data 


reported here for the GR-S rubbers in the Mooney 
viscometer were obtained by a procedure that meets 
the following three requirements: (1) the rubber 
should be allowed to come to the temperature of the 
instrument before data are taken; (2) the data at all 
rates of shear should be taken on material in as 
nearly as possible a uniform state of degradation; 
and (3) conditions should simulate as nearly as pos- 
sible those of the standard Mooney measurement. 
The procedure selected is as follows. A fresh, mill- 
massed sample is placed in the instrument, which 
has previously been heated to 100° C. Ten minutes 
is allowed for the rubber to reach the temperature of 
the instrument. The viscometer is then started at 
2 rpm. At the end of 4 min the gage reading is 
recorded, and the speed is reduced to 1.5 rpm. As 
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soon as the gage reading becomes steady at the ney 
speed (less than 1 min) it is recorded and the spee 
reduced again. Further readings are made at su 
cessive speeds of 1.0, 0.5, and 0.2 rpm. The speed 
used correspond to rates of shear from 0.103 to 1.0 
sec™', As GR-I rubbers do not show the aboy 
molecular degradation [5, 12], this procedure is no 
followed so exactly with them, and measurement 
are made at higher rates of shear. 

For the Mooney viscometer the observed rag 
reading (@) and rate of revolution (W) are converte: 
to a volume average viscosity (y) and rate of shea 
(y) by the approximate method given by Moone, 
[4]. Mooney found that a purely theoretical anal 
ysis had to be empirically corrected to give agree- 
ment with results from the “rubber rheometer” [13] 
a viscometer of more nearly uniform rate of shea: 


developed by him. This correction is used here 
The numerical relations used here are 
7=34,000G W’, and 7 0.516W 3 


IV. 


Results on the six rubbers with the worker-consis 
tometer and with the Mooney viscometer, used 
according to the procedures described above, are 
given in figures 1 through 6 as logarithmic plots of 
viscosity at 100° C versus rate of shear. These 
figures also show the viscosities calculated from the 
standard Mooney readings. 


Results 
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within 10 percent over their common range. The 
standard Mooney value indicated by the x agrees 
with the Mooney curve within 3 percent. The 
two curves have slightly different slopes in their 
common region. However, an extension of the 
Mooney curve to slightly higher rates of shear 


appears to show better agreement with the worker- 
consistometer curve. 
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For X-558 GR-S (fig. 2) all comparable data agree 
within the uncertainty of the measurements. 

For X-478 GR-S (fig. 3) the only Mooney obser- 
vation is the standard Mooney value. This value 
agrees with the worker-consistometer curve with- 
in the uncertainty of the measurements. 
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FiGuRE 2 Viscosity of X-558 GR-S versus rate of shear 


X, Standard Mooney value , worker-consistometer values; @, Mooney 


The results for the three GR-—I rubbers (figs. 4, 5, 
and 6) are sufficiently similar so that they may be 
mentioned collectively. For these rubbers the 
standard Mooney values are in complete agreement 
with the present Mooney curves. This result would 
be expected since no degradation occurs under the 
conditions of the experiments [5, 12], and it is pre- 
sumably degradation that causes discrepancies here 
for the GR-S rubbers. In the region around 1 
sec~' the Mooney values are all about 20 percent 
higher than the corresponding worker-consistometer 
values. A Mooney curve and a corresponding work- 
er-consistometer curve do not have identical slopes 
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over their common region in every case. However, 
from figure 5 it appears that an extension of the 
Mooney curve to higher rates of shear may coincide 
with the worker-consistometer curve. 
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FIGURE 4 Viscosity of Y-105 GR-I-17 versus rate of shear 
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Fiaure 9. Flow curves for GR-S and GR-] rubbers. 
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Figure 7 is a composite plot of the worker-cor 
sistometer data for the GR-Srubbers. It is include 
as a demonstration of the differences in the flo. 
curves for GR-S rubbers. The plots for X-—51 
GR-S and for X-478 GR-S have approximately th 
same shapes, but different locations. The plots fo 
X-558 GR-S and for X-478 GR-S have roughly th 
same locations, but different shapes. 

Figure 8 is a similar composite plot for the GR 
rubbers. It is interesting that although the flow 
curves for these three rubbers are separated at lowe 
rates of shear they are nearly coincidental at highe 
rates of shear. 

In order to give a direct comparison of the data fo 
a GR-S rubber with that for a GR-I, the results fo: 
X-558 GR-S and for Y-105 GR-I are plotted to 
gether in figure 9. The viscosity of the GR-] 
rubber shows the most rapid decrease at lower rates 
of shear than does that of the GR-S. 


V. Discussion 


Since the differences shown in figures 1 to 6 be- 
tween the results from the two instruments for 
identical rubbers are not the same for all the rubbers, 
there is no constant instrumental difference between 
the instruments. 

In the case of the GR-S rubbers the variation of 
the differences may represent varied degrees of deg- 
raduation produced in different rubbers by the same 
procedure. Treloar [10] shows what large variations 
may be produced in the location of a flow curve by 
procedural differences. Unreported data of our own 
confirm this fact and also show that large differences 
in the slopes may be thus produced. 

The differences between the results from the two 
instruments for the GR-I rubbers show a large scale 
similarity within the group. In all three cases the 
Mooney values are about 20 percent higher than the 
worker-consistometer values. Degradation cannot 
be a factor in these differences. However, since the 
differences do not show, sufficient uniformity to be 
attributed completely to constant instrumental dif- 
ferences, they are apparently also contributed to by 
properties of the rubbers other than viscosity, such 
as elasticity, and adhesion to the viscometer surfaces. 
These properties may vary for the different rubbers 
and thus produce varied discrepancies. Known 
differences and variations in the hydrostatic pres- 
sures, and possible ones in the temperatures in the 
two instruments, would also result in differences 
that would not be related to the flow mechanics and 
that would vary for different rubbers having unlike 
pressure and temperature coefficients of viscosity. 
Another possible source of difference is the flow 
analysis for the Mooney viscometer, which is only 
approximate and may not take adequate account of 
the change of viscosity with rate of shear. Therefore 
the Mooney results on two rubbers having different 
flow curves may well compare differently with 
worker-consistometer results for the same _ two 
rubbers. 
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In the case of figure 1, a tentative comparison 
with literature can be made. Treloar [10] gives 
flow data obtained with a biconical rotor Mooney 
viscometer |6] for a rubber described as a “standard 
GR-S”. Although this rubber is not described 
further, and the preparation of the sample is de- 
scribed only as “a preliminary light milling’’, the 
results probably represent a material similar to that 
treated in figure 1. Treloar’s data have been ex- 
pressed in the form of figure 1 and are plotted on 
that figure (@). Good agreement with the present 
results is shown, but it may be entirely fortuitous in 
view of the extreme dependence of results on the 
preparation of the sample. 


VI. Conclusion 


It can be concluded that in the worst cases there 
is agreement within about 20 percent between results 
from the worker-consistometer and those from the 
standard Mooney viscometer used according to the 
procedure described here. For two of the rubbers 
investigated the differences are within the limits of 
uncertainty of the measurements. The procedure 
used was selected as giving the best possible simula- 
tion of the conditions of the standard Mooney 
measurement, consistent with a valid indication of 
the temperature of the sample. The degree of agree- 
ment shown is within reasonable expectations in 
view of the confusing factors mentioned above. 

It has also been demonstrated that the flow curves 


for different rubbers of a given type can be of notice- 
ably different characters. 


Many of the calculations and all of the drawings 
for this paper were made by Miss H. V. Belcher. 
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Spectrophotometric Determination of Carboxy]l in 
Cellulose’ 


Florence H. Forziati, John W. Rowen, and Earle K. Plyler 


{ method that uses the intensity of the absorption of the oxidized celluloses at 5.8 


microns as a measure of their carboxyl content is described. 
on mineral oil slurries of 1:1 mixtures of finely divided polystyrene and cellulose. 


The analyses are carried out 


The 


absorption of the polystyrene at 6.2 microns is used as a reference band, thus obviating the 


necessity of accurately controlling the thickness of the absorption cell 


The method gives 


results that agree fairly well with those obtained by the calcium acetate method and can be 
used to detect as little as 0.1 millimole of carboxy! per gram of cellulose. 


I. Introduction 


Cellulose products are frequently exposed to oxi- 
dative environments during manufacture and use. 
The oxidation of cellulose, commonly represented 
by the formula below is, therefore, of great technical 
importance. 


The types of oxidation theoretically possible for 
cellulose are: 

Oxidation of primary alcohol groups (carbon 6 of 
the anhydroglucose units) to aldehyde or carboxy! 
groups. 

Oxidation of secondary alcohol groups (carbon 2 
or 3) to ketone groups 

Oxidation of glycol groups (carbons 2 and 3) to 
aldehyde or carboxyl groups with attendant cleavage 
of the carbon-to-carbon bond. 

Oxidation of hemiacetal groups (carbon 1) to 
carboxyl groups. Although this type of oxidation 
is unimportant in undegraded celluloses where the 
number of terminal groups is extremely small, it 
may be of importance in certain degraded celluloses. 

With few exceptions (] to 5|,7 oxidants appear to 
be nonselective in their action on cellulose, giving 
rise to two or more of the above types of oxidation 
simultaneously. Accordingly, for the characteriza- 
tion of oxidized celluloses, methods capable of detect- 
ing, estimating, and allocating ketone, aldehyde, and 
carboxyl groups in the presence of each other are 
needed. Available chemical methods are not entirely 
satisfactory for this purpose 

The presence of carbon attached to oxygen by 
means of a double bond in the aforementioned groups 
suggests that oxidized celluloses will show the char- 


A report of work done under cooperative reement with the United States 
Department of Agriculture and authorized by the Research and Marketing Act 
he work is being supervised by the Southern Regional Research Laboratory 

of the Bureau of Agricultural and Industrial Chemistry F 
? Figures in brackets indicate the literature references it the end of this paper. 


acteristic absorption of the C=O stretching vibra- 
tions in the 5.5- to 6.5-u region of the spectrum 
As the frequency and intensity of this absorption 
are normally influenced by the adjacent molecular 
structure [6, 7], ketonic, aldehydic, and carboxylic 
carbonyl groups of oxidized celluloses might reason- 
ably be expected to absorb at different wavelengths 
If such is the infrared spectrophotometry 
should be useful in the analysis and characterization 
of oxidized celluloses. 

During earlier work in this laboratory, an absorp- 
tion band at 5.8 uw was observed in the spectra of 
films of regenerated cellulose that had been sub- 
jected to moderate oxidation with nitrogen dioxide 
[8, 9]. Under the conditions used, this reagent is 
selective in its action on cellulose, resulting only in 
the conversion of primary alcohol groups to car- 
boxyl groups [1, 2]. The absorption of these films 
at 5.8 w was, consequently, attributed to carboxylic 
carbonyl. Subsequent examination of spectra of 
mineral oil slurries of cottons that had been treated 
with a vari¢ty of oxidants showed absorption at 
5.8 uw to be of frequent occurrence. Although it was 
recognized that noncarboxylic carbonyl might be 
responsible for this absorption in some cases, the 
possibility of using the intensity of the absorption 
of the oxidized celluloses at 5.8 uw as a measure of 
their carboxyl content was investigated. The re- 
sults of this investigation are presented in this 
report. 

The use of the intensity of an absorption band 
as a measure of the concentration of the group or 
compound responsible for the band is based on 
Beer's law. This law states that log J,//=ked, 
where /, and J are the intensities, respectively, of 
the radiant energy of a given wavelength falling 
upon and passing through the sample, & is an absorp- 
tion coefficient characteristic of the sample, ¢ is 
concentration, and d is the cell thickness. In the 
analysis of discontinuous solids, such as the oxidized 
celluloses, the necessity for measuring cell thickness 
is usually avoided by the addition of a properly 
chosen solid reference material in a fixed ratio to 
the total sample [7, 10]. Spectrophotometric meas- 
urements on mixtures of the reference material with 
samples containing varying and known amounts of 
the component of interest supply the data required 


case, 
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for the construction of a working curve to be used | 
in the analysis of unknowns. The curve is obtained 
by plotting the concentration of the component of 
interest against the ratio of log J,/J for the wave- 
length of interest in the sample to log J,// for a 
wavelength at which the reference material absorbs. 


II. Construction of an Analytical Working 


Curve 


In the present work, data for the construction of 
an analytical working curve were obtained from | 
spectra of mineral oil slurries of 1:1 mixtures of 
cellulose samples of known carboxyl content with 
finely divided polystyrene, which provided a refer- 
ence absorption band at 6.24. The polystyrene was 
obtained from the manufacturer, the Dow Chemical 
Co., in granular form. Prior to use, it was ground 
to a fine powder in a vibratory ball mill [11], washed 
with 95 percent ethanol and dried in air. 

Empire cotton purified by a modification of the 
procedure described by Worner and Mease [12] and 
milled for 30 min in the vibratory ball mill, was used 
as a diluent in the preparation of cellulose samples 
of known carboxyl content. A_ nitrogen-dioxide- 
oxidized cotton prepared by the Tennessee Eastman 
Corp. was used as the source of carboxyl. This 
oxidized cellulose had been ground in a Wiley mill 
to pass a 100-mesh screen, freed of water-soluble 
material by repeated extraction with water, and 
dehydrated by successive treatments with absolute 
methanol and benzene. The benzene was removed 
by prolonged storage over paraffin in vacuum. 

Before the cotton, oxidized cotton, and 
polystyrene were conditioned at 65-percent relative 
humidity at 21°C. The carboxyl contents of the 
conditioned cotton and nitrogen-dioxide-oxidized 
cotton as determined by the calcium acetate method 
[1] were 0.01 and 1.78 millimoles per gram. This 
sample of nitrogen-dioxide-oxidized cotton was used 
in preference to a more highly oxidized sample in 
order to keep the concentration of noncarboxylice 
carbonyl to a minimum [1,13]. The quantities of 
cotton and oxidized cotton used were such that the 
total weight of each cellulose sample was 100 mg. 

The cellulose samples were dry-mixed with an 
equal weight of polystyrene and were then mulled 
with 0.5 ml of mineral oil between weighted ground 
glass plates. For the absorption measurements, the 
resulting slurries were injected into demountable 
rock-salt absorption cells equipped with silver or 
cellulose acetate spacers of approximately 25, 50, or 
100 w in thickness. The spacers were chosen in such 
a way that the values for the transmission at 5.8 u 
ranged from 30 to 60 percent. 

Spectra were obtained with a double-beam auto- 
matic recording Baird spectrophotometer [14] 
equipped with a sodium chloride prism. The per- 
formance of the instrument was checked frequently 
against a polystyrene film of known spectral charac- 
teristics. The single-cell rather than the double- 
cell method was used in securing spectra. In order 
to compensate for variations in the behavior of the 


use, 


spectrophotometer, four to six spectra were recorded 
for each absorption cell. Values for J) at 5.8 and 
6.2 w could not be obtained directly from the spectro- 
grams. They were, therefore, assumed to be equal 
to the percentage transmission at nearby wave- 
lengths at which neither oxidized cellulose nor 
polystyrene has absorption bands—namely, 5.5 and 
6.4 4. Average values for the ratio, 


(To/D)s. 0 
(1h I 16 2u 


were used in the construction of the working curve 
and in spectrophotometric analyses. 

In figure 1 are shown the 5.2- to 6.5-» region of 
the spectra of mineral oil slurries of three of the 
polystyrene-cellulose mixtures. For comparison, 
spectra of mineral oil slurries of the polystyrene, 
cotton, and nitrogen-dioxide-oxidized cotton are 
given in figure 2. In table 1, the values of the ratio, 
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lulose samples. A comparison of the plot of these 
these data with the straight line obtained by the 
application of the method of least squares to the 
experimental data, both presented in figure 3, shows 
that 


Tol Ds. 
To/D) 6.2 


log 


log 


is essentially a linear function of the concentration 
of carboxyl in the cellulose samples. The actual 
values of the points plotted have an uncertainty of 
about one-twentieth of a unit 

The failure of the straight-line plot to pass through 
the origin as required by Beer’s law is the result of 
mixtures 


TABLE 1 measurements on 1:1 


of polystyrene with cellulose of known carboryl content 


S pe ( trophotome frie 


Carboxvl cor Carboxyl cor log (lo/T)s 4» 
ent of cellulose / ‘ t of cellulose log (Jo/1)e.2 » 
Millimoles Millimol 

0.01 0.1 7a 1. 22 

ol 2 su 1, 32 

ol il xu 1.18 

" mS Tr) 1.15 

i 4! 16 1. 38 

f ou 1. 52 1.84 

‘ fis 1.71 2. 26 

4 70 7 2. 33 

i s4 1. 78 2.12 
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the arbitrary procedure employed in obtainin 
values for /) from the spectrograms. It should n 
affect the usefulness of the plot as an analytic: 
working curve 


III. Analysis of Oxidized Celluloses 


As a check on the general utility of the analvtica 
working curve, it was used in the analysis of | 
oxidized Using the procedure outline 
above, values for the ratio, 


celluloses 


log 


€ 
log CU, 


} > 
2M 


D) 5s 
I 
were obtained for 1:1 mixtures of polystyrene with 
the oxidized celluloses The concentrations of 
carboxyl in the oxidized celluloses were then read 
directly from the working curve. 

The oxidants used in the preparation of the oxi 
dized celluloses were air at 145° C, chromic oxide in 
a mixture of glacial acetic acid and acetic anhvdrid 
115], hvdrogen peroxide [16, p. 114], potassium 
dichromate in oxalic acid solution [16, p. 114}, 
potassium permanganate in sulfuric acid solution 
|16, p. 112], bromine in the presence of excess barium 
carbonate [17], buffered periodic acid [18, 19], and 
nitrogen dioxide [1!. With the exception of the two 
nitrogen-dioxide-oxidized and one of the 
periodate-oxidized cottons that were prepared, re 
spectively, by the Tennessee Eastman Corp. and 
R. E. Reeves of the Southern Regional Research 
Laboratory, all of the oxidized cottons were prepared 
from purified Empire cotton that had been ground in 
a Wiley mill to pass a 20-mesh screen. The fact 
that the cotton was added to the oxidant without 
previous swelling probably accounts for the low 
degree of oxidation observed in most samples 

The results of the spectrophotometric analyses are 
riven in table 2. Also given are the values obtained 


cottons 


for carboxyl and noncarboxylic carbonyl by the 
calcium acetate [1] and hydroxylamine hydro- 
chloride [15] methods. In table 3 are given the 


results of 13 analyses for carboxyl carried out on a 
single oxidized cotton over a 3-month period by 
2 different analysts as a check on the precision of the 
spectrophotometric method. The values obtained 
for carboxyl by the spectrophotometric and calcium 
acetate methods agree fairly well. The greatest 
difference observed, 0.34 millimoles per gram for a 
nitrogen-dioxide-oxidized cotton, is of doubtful 
significance, as the yellow color of this oxidized 
cotton introduced considerable uncertainty in the 
determination of carboxy! by the calcium acetate 
method. Difficulty was also encountered in applying 
the calcium acetate method to the periodate- 
oxidized cottons. Thus, here also the value obtained 
by the spectrophotometric method may be the more 
reliable of the two. In general, the results indicate 
that the intensity of the absorption of the oxidized 
celluloses at 5.8u constitutes a fairly satisfactory 
measure of their carboxy! content even in the 
presence of appreciable amounts of non-carboxylic 
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carbonyl. When measured according to the procedure 
described in this report, as little as 0.1 millimole of 
carboxyl per gram of cellulose can be detected. 


TABLE 2. Analyses of oxidized celluloses 


Noncar- 
boxylic 
Carboxyl Carboxyl carbonyl 
»¥ cal- Vy spec- ) ° 
Oxidant used in ooo y &-, ek meat Devia- én. 
of oxidized cellulose acetate metric tion mine 
method | method hydro 
chloride 
method 
Milli- \Milli- Milli \Milli- 
mole/g mole/g mole/g mole/g 
Air at 145° C 0. 02 0. 04 +). 02 0. 40 
CrO;/acetic acid-acetic anhy- 
dride 00 02 +. 02 . 95 
H.O 4 06 +. 02 as 
K»Cr?O7/oxalic acid OF 02 —. 02 44 
K Mn04/H2S80, wy 56 +. 12 18 
Bre/BaCoO 38 40 +. 02 1.19 
Bre/ BaCOo 66 6S +. (2 0. 38 
H1O, 02 19 +.17 7. 37 
HIO, 09 19 +. 10 
NO 1. 40 1.39 ol 
NO 2. 65 2. 31 34 
TARLE 3 Spectrophotometric analyses of potassium-perman- 


ganate-oxidized cotton 


Deviation from 


Sample Carboxy! 
mean 
Millimole/g Millimole/g 
0. 46 0.10 
2 59 + .03 
ci 00 
1 fil 0 
55 ol 
‘ il O5 
52 04 
x 53 0 
Qu »S ths 
TT 5s rd 
11 53 o 
12 ii 10 
13 0) 0 
0. 46 mean +(). 04 average de- 
viation from 
mean 


IV. Summary 


Details of a spectrophotometric method for the 
quantitative estimation of carboxyl in cellulose are 
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given. The method gives results that are in fair 
agreement with those obtained by the calcium 
acetate method and can be used to detect as little 
as 0.1 millimole of carboxyl per gram. 


The authors express their appreciation to Walter 
K. Stone, Edith T. L. Voong, and Eva Estermann 
for technical assistance. 
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Static Friction Tests with Various Metal Combinations 
and Special Lubricants 


H. S. White and Dino Zei 


An inclined-plane apparatus was used in obtaining static friction data for various 


lubricants with different 


tives, 


petroleum compound, was investigated 


metal combinations 
fluorinated hydrocarbons, were compared with a reference mineral oil 
such as graphite, molybdenum disulfide, zinc oxide, boron nitride, and an oxidized 
Metal 


such as chlorinated and 


The effect of addi- 


Special lubricants, 


combinations involving stainless steels, 


carbon steel, cast iron, aluminum alloy, and chromium plate were used 
It was found that the special chlorinated and fluorinated lubricants without additives gave 


higher friction than the Navy symbol 2110 reference oil 


A lubricant containing the oxidized 


petroleum compound and one containing a silicone grease gave slightly less friction than the 


reference oil 


With a mineral lubricating grease, with dry molydenum disulfide, and with 


lubricants containing graphite or molybdenum disulfide, the friction was about 50 percent 


less than with the reference oil 
reference ol 


The remaining lubricants gave higher friction than the 2110 


In general, the heat-treated stainless-steel combinations gave the lowest friction, whereas 


the combinations involving cast iron or aluminum alloy gave relatively high friction 


Fric- 


tion was less with lapped and polished surfaces than with ground surfaces 


I. Introduction 


Static friction data were obtained for making a 
practical comparative evaluation of the oiliness of 
special lubricants when used with various metal 
combinations. Navy symbol 2110 oil was used as 
a reference lubricant for comparing the friction 
obtained with the special lubricants. The work was 
done at the request of the Naval Ordnance Labora- 
tory, which furnished most of the lubricants and the 
unmachined metals and specified the metal combina- 
tions to be used. ‘The time allotted did not permit 
a comprehensive investigation, but it is believed that 
the data may be helpful to those interested in the 
use of these types of lubricants and metal combina- 
tions. 


Il. Apparatus 


An inclined-plane apparatus was used for making 
the static friction tests. The coefficient of friction, 
f, is defined as F/W, where F is the tangential force 
at the sliding surface, and W is the force normal to 
the surface. As is customary with this type of 
apparatus, the tangent of the angle with the hori- 
zontal when a rider on the incline starts to slide is 
used as a measure of f. For a given vertically 
applied load, the magnitude of W changes with the 
angle of the inclined plane, and the distribution of 
the load at the sliding surface varies with the angle. 
However, these variations are not important, as the 
coefficient of friction was found to be practically 
independent of load over a wide range. 

The testing machine is shown in figure 1. An 
electric hot plate is mounted on transite between two 
journals. A handwheel acts through a_ 100-to-1 
reduction gear to give slow tilting of the hot plate 
A pointer, perpendicular to the hot plate, indicates 
ona calibrated scale the tangent of the angle between 
the surface of the hot plate and the horizontal. The 
heavy cylindrical weight at one corner of the transite 


mounting is used to load the reduction gear in on 
direction, thus eliminating backlash. The plate 
specimen is fastened to the hot plate with two screws 
A loading yoke with a weight pan and weights gives 
the desired load on the rider specimen. ‘Two steel 
bars (thicker than the plate specimen), fastened to 
the hot plate, have several pairs of tapped holes to 
allow fastening of cross bars at different positions 


One of these cross bars (at the left in fig. 1) is pro- 


Figure 1. Static friction machine. 
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steel 
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pro- 


vided with an electrically insulated screw and needle 
that are connected in series with a dry cell, lamp, 
and rider. Connection from the lamp to the rider 
is made with a No. 38 copper wire in such a way as 
not to influence movement of the rider. When the 
rider contacts the needle, the completed circuit is 
indicated by the signal lamp. The temperature, 
controlled manually, is measured by a thermocouple 
fastened to the plate specimen. 

Typical test specimens are shown in figure 2. The 
plate dimensions are * in. by 2 in. by 7 in., and the 
rider dimensions are % by 1 by 2in. Each rider has 
a cone-shaped depression formed in the center of its 
upper surface for locating the point of the loading 
voke. The central portion of the underside of each 
rider is undercut so as to leave a bearing surface 
y by lin. ateachend. The tapped holes for fasten- 
ing the thermocouple to the plates and the screws for 
fastening the signal light wire to the riders are shown 
in figure 2. The rider and plate shown inthe back- 
ground were finished by grinding. After the plate 
shown in the foreground was ground, it was lapped, 
and the contacting surfaces of the rider were finished 
by lapping and polishing. 


III. Metal Combinations 


The code numbers and descriptions of the plate 
specimens are given in the following tabulation. 


Remarks 


< 
~ 


( 
Roughness 


Dry-ground 
Dry-ground and lapped 

Do 
Kerosine-ground 
Kerosine-ground and lapped 
Oil-ground and lapped 


The code numbers and descriptions of the rider 
specimens are given in the following tabulation. 


Rock well 
rms 


( 
Roughness 


Hardness, 


» st Dry-ground 
do.* ¢ Dry-ground, lapped and pol 
ished 
Chromium plate Do 
4 “4 stainless steel 5 | Dry-ground 
o4-LP do Dry-ground, lapped and pol 
ished 


sinless stee] ® 


1050-LP 1050 s Do 


® Sulfurized for free machining 


FiGurRE 2. Static friction machine specimens, 


The specimens of stainless steel 431, 440C, and 
416 were hardened by heating to 1,850° F and quench- 
ing in oil, and tempered by heating at 500° F for 1 hr 
and cooling in air. 

The sharp edges of the by 1-in. contact surfaces 
of the riders were stoned and polished to prevent any 
shaving effect. 

The values of roughness given are in microinches 
(rms) as determined by an Abbot Profilometer. 
With the lapped cast iron (2—L) the open pores gave 
erratic high readings above 2 yu in., but examination 
with an interference microscope indicated a rough- 
ness of about 2 wm in., if the pores were neglected. 
With the ground riders, roughness was measured on 
the backs of the specimens, which were ground in 
the same manner as the contacting surfaces. It is 
assumed that the lapped and polished surfaces of 
the riders had a roughness less than 1 win. In the 
code numbers, L means lapped, and LP means 
lapped and polished. 

A cast iron plate with grids was used to charge a 
similar lead-tin-alloy plate, which was then rinsed 
and used to lap the test specimens. This method 
prevented charging of the specimens with abrasives. 
After lapping the riders were polished on 4/0 polish- 
ing paper resting on a surface plate. 

The specimen combinations used in these tests are 
listed in the following tabulation. 


$1 

14 

s1-~Land 440C-L 
s1-L and 214-L 
14-L 

-L 





IV. Lubricants 


The lubricants used in this investigation are listed 
in the following tabulation. 


Lubricant 


Navy symbol 2110 oil (SAE 10 

Chlorinated hydrocarbon 

Organic phosphate ester 

Fluorinated hydrocarbon 8 (medium viscesity 

Fluorinated hydrocarbon 8+1° commercial 
uiditive, oxidized petroleum compound 

Fluorinated hydrocarbon 8+1% molybdenum 
disul fice 


Fluorinated hydrocarbon FS (light viscosity 

Fluorinated hydrocarbon 8+1°) graphite 

Perfluorokerosene 

Fluorinated hydrocarbon 
phthalocvanine 

Silicone DC-200 (100 centistokes at 25° C 


FS+0.5% copper 


Fluorinated hydrocarbon FS+2°; silicone grease 
DC-—44 (medium 

Fluorinated hydrocarbon KF L-F 

Fluorinated hydrocarbon FS+5* 
grease (light 

Cirease G (23.5¢ 

Cirease G+-5°, 


fluorinated 


sodium soap 
molybdenum disulfide. 


Molybdenum disulfide (dry 

2110 oil+-1°% molybdenum disulfide 

Fluorinated hydrocarbon FS+1°% zine oxide. 

Fluorinated hydrocarbon S+2°, molybdenum 
disulfide concentrate 

Fluorinated hydrocarbon 8+0.75* 
amorphous 


> boron nitride 


was selected as the reference oil for 
frictional characteristics of the 


Lubricant 7 
comparing the 


samples, most of which were furnished by the Naval 


Ordnance Laboratory. The lubricants were received 
in groups and each group was assigned a series 
number. 


V. Test Procedure 


For cleaning the test specimens, five solvents were 
used in the following order: Stoddard solvent; 
carbon tetrachloride; benzene, cp; acetone, cp; ethyl 
alcohol, 95 percent. Absorbent cotton was used for 
scrubbing the specimens with the first two solvents. 
This was followed by rinsing with the other three and 
drving with clean cotton. 

The liquid lubricants were applied to the plates 
and riders by means of small clean cotton swabs. 
Grease samples were applied with a finger tip and 
dry molybdenum disulfide by dusting. 

With the plate horizontal, the rider was placed in 
position, contacting the insulated needle, and the 
load was applied. The needle was retracted 0.006 
in. (4% turn of the screw), and the plate was tilted 
slowly until the rider moved sufficiently to contact 
the needle and light the signal lamp. The tangent 
of the angle of tilt was recorded. <A second reading 
was obtained by leveling the plate, retracting the 
needle another 0.006 in., and again tilting the plate 
until the signal lamp was lighted. This was repeated, 
usually until eight readings were obtained, giving a 
total movement of the rider of 0.048 in. This pro- 
cedure and series of readings were considered as one 
test. Between tests the specimens were cleaned and 
the needle was returned to its original position so 


that the next lubricant could be tested on the sar 
area of the plate specimen. Because of the variati 
in static friction resulting from such things as t 
presence of the oxide films, microscopic irregularit; 
of the metal surfaces, and contamination from t| 
atmosphere, tests with the reference oil (Navy sy: 
bol 2110) were made in conjunction with each seri 
of special lubricants. For example, with series No 
and one needle bar position, tests were made wit 
lubricants in the following order of Code Nos.: 7, 
13, 7, 14,7, 15,7. After all the lubricants of a seri: 
had been tested at one plate area, the needle bar w: 
moved to another position, the rider was lapped an 
polished again, and the series was tested at the ne 
plate area. In general, these comparative tests were 
made at three different areas of the plate specimen 
With each new surface area, two or more tests were 
made with the 2110 oil to assure reasonable repro 
ducibility before testing the series of lubricants 
After three areas had been used on a plate specimen, 
it was lapped again before being used with the next 
series of lubricants. 


VI. Test Results 


Some preliminary testing was done in order to 
study reproducibility and the effects of temperature 
and load. The effect of temperature was small and 
was within probable experimental error for the range 
of 80° to 120° F. Tests with loads of 5, 10, and 15 
lb. indicated that variation in the load had no effect 
on the coefficient of friction, that is, the tangent of 
the angle when sliding occurred. This was checked 
later with the 2110 reference oil and special lubricants 
16 and 17, with the 416S—-LP rider on the 431-—L 
plate. All the values of coefficient of friction report- 
ed herein were obtained at a nominal temperature of 
100° F and with a 10-lb. load (40 lb/in.? average). 
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Representative friction data at 10-lb load and 100 
F, with ground surfaces. 


Figure 3. 
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LUBRICANT CODE NUMBER 
FicurE 4. Representative friction data at 10-lb load and 100° 
F, with 416S-LP rider on 431-L plate (first series of 
lubricants). 


Some tests were made with ground surfaces on 
the metal specimens. Typical data for the 2110 
reference oil and lubricants 13, 14, and 15 are shown 
in figure 3. The horizontal lines represent the 
averages for each group of data obtained in a test. 
The averages for the reference oil are connected by 
broken lines from which vertical lines are drawn 
to the averages for the special lubricants. 

Due to the large spread of data on ground surfaces 
and the fact that lapped or lapped-and-polished 
surfaces are more easily reproduced and are more 
commonly used for machine parts involving small 
clearances, subsequent tests were made with lapped 
plates and lapped-and-polished riders. The lapping 
and polishing processes have already been described. 

Figure 4 shows typical data for the first series of 
lubricants with the 4168 lapped-and-polished stain- 
less-steel rider on the 431 lapped stainless steel plate. 
Where a test was obviously not representative it was 
ignored, as indicated by the question mark for one 
of the tests with oil 7. Otherwise, the spread of 
data in figure 4 is much less than in figure 3. 

Representative data for the second and fifth series 
of lubricants are shown in figures 5 and 6, respec- 
tively. 
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Ficure 5. Representative friction data at 10-lb load and 100° 
F, with 416S—LP rider on 431-—L plate, at two locations on 
the plate (second series of lubricants). 


—| | metal combination are given in table 1. 


| 
| 
| 
| 
| 


The data for all the lubricants with this particular 
These 
were obtained from plots similar to figures 4, 5, and 
6. The position number refers to the position of the 
needle bar, and hence each different position means 
a different location of the rider on the plate specimen. 
The spread of data is indicated by giving the maxi- 
mum and minimum readings for each test, and the 
average corresponds to the horizontal mean lines 
as in figures 4, 5, and 6. The percentage increase 
in friction corresponds to the vertical lines connecting 
the means to the broken lines for the 2110 reference 
oil. The percentage increase is indicated as nega- 
tive when the friction is less than with the reference 
oil. The reference values of f used for the 2110 oil 
in connection with each special lubricant test is 
given in the last column. 
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Ficure 6. Representative friction data at 10-lb load and 100° 
F, with 416S—LP rider on 431-—L plate (fifth series of lubri- 
cants). 


TaBLeE 1. Static friction data at 10-lb. load and 100° F. for 
4i6 stainless steel (416S—LP) on 431 stainless steel (4,31-L) 


Lubri f (tangent) , 
Percent 


Test Position cant age in- f 2110 
No No Code ntact. | Mini- Oil 


crease ® 
verage 
mum mum Average 


0. 210 0. 192 
190 Is4 
146 


221 

218 

178 

191 

164 

138 

162 

150 130 
149 139 
OSY O78 
O49 O34 
. 060 043 
219 194 


* Percentage increase in average f relative to tests with 2110 oil 
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TaBLe 1. Static friction data at 10-lb. load and 100° F. for Tables 2, 3, 4, 5, and 6 present similar data 
stainless stee iS *\ or tT ste less steel (438 . . : 
Go { ating sane el (431-L the other metal combinations. 
ontinues 7 ; 
The data in the above tables are summarized 
' Pee, - 
figure 7. The average values of f for lubricant 
Lat Pere (2110 reference oil) and special lubricants 13 to 
osit ‘ t 2110 ° . . . . rr 
Ni ( - * age | 0 are given for six metal combinations. The averag 
N mt . for lubricants 21 to 25 and 31 and 32 are given | 
the two metal combinations used with them. Lub 
: : 4 cants 26 to 30 were prepared at this Bureau and we 
tested with one metal combination only. 
) ~ ) 2 14 : 
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Static friction 
304 stainless steel (304 


N 

nm 
108 7 13 { 
116 5 13 


148 } l 


110 14 
118 14 
150 4 
112 7 l 
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VII. Conclusion 


For static friction under the conditions used for 
these tests, the following indications have been 
obtained: 

1. Polished surfaces sliding on lapped surfaces 
tend to give lower friction and less variation than 
ground surfaces sliding on ground surfaces. 

2. For the metal specimens tested, with most of 
the lubricants, the heat-treated stainless steel com- 
binations gave the lowest friction. The chromium- 
plated rider on the heat-treated stainless steel plate 
also gave correspondingly low friction with lubri- 
cants 16,17,and 19. For the chromium-plated rider 
and the 304 stainless steel rider on the Alcoa 214 
alloy plate, and for the 1050 steel rider on the cast 
iron plate, the relative friction varied greatly with the 
lubricant. Each of these combinations gave high 
friction with one or more lubricants while giving low 
friction with other lubricants. 


3. The lubricants containing molybdenum disulfic 
gave the lowest friction, whereas mineral lubricati 
grease, dry molybdenum disulfide, and the lubrica: 
containing graphite gave friction nearly as loy 
With the above lubricants the friction was about 
to 60 percent less than with the 2110 reference o 
The lubricant containing an oxidized petroleum con 
pound and the one containing a silicone grease gay 
slightly less friction than the 2110 oil. The chloriy 
ated and fluorinated lubricants without additiv« 
and the remaining special lubricants gave high: 
friction than the 2110 reference oil. 

These relative evaluations of friction apply onl 
for static conditions as used in these tests. Any sig 
nificance with respect to kinetic boundary lubricatio: 
is questionable, except possibly in those cases wher 
the velocity approaches zero. 


WasHINGTON, October 9, 1950. 
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Apparatus for the Preparation of Anhydrous Titanium (III) 
Chloride and Titanium (III) Bromide 


J. M. Sherfey 


{n improved method is described for the preparation of titanium (III) chloride or 
titanium (III) bromide from the corresponding tetrahalide. The apparatus consists essen- 
tially of a reaction kettle, in the upper half of which is suspended a hot tungsten filament. 
\ mixture of hydrogen and the vapor of the refluxing tetrahalide reacts at the surface of 
the filament, and the resulting solid product collects at the bottom of the vessel 


The laboratory preparation of anhydrous titanium 
(111) chloride has been described by several workers. 
The method described by Georges and Stahler' has 
been perhaps the most satisfactory. A stream of 
hydrogen is saturated with titanium (IV) chloride, 
and the mixture is passed through an annular space 
formed by two concentric tubes. The inner wall is 
electrically heated, and the outer wall is cooled by 
a water jacket. The mixture reacts at the hot sur- 
face, and the resulting trichloride collects on the 
cold surface. It was found that the rate of produc- 
tion of titanium (III) chloride in the apparatus 
described was only a few grams per hour and that 
the yield, particularly with regard to the hydrogen 
used, was very poor. The product was contami- 
nated with tetrachloride, especially if the tempera- 
ture of the tetrachloride saturator was high. 

The apparatus described below, utilizing the same 
reaction (2TiCl, t H.: “2TiCl; +-2HCl), has been 
found very satisfactory, and 150 to 200 g can be 
produced in| day. The product is about 98 percent 
pure and can be further purified 

The essential parts of the apparatus are shown in 
figure 1. The 4-liter reaction kettle, A, has a flanged 
top with four female standard taper openings, one 
34/45) centrally located with three (24/40) joints 
symmetrically spaced around it. For simplicity 
these are shown in the drawing as if they were in a 
straight line. 

The center joint, B, of the reaction kettle is fitted 
with a chimney, C, in the top of which is a rubber 
stopper, D. Two tungsten rods, E, about 6 mm in 
diameter, pass through this stopper and support a 
tungsten filament, F. The diameter and length of 
this filament are fixed by the power supply available. 
A filament 1 mm in diameter and 30 em long re- 
quires 38 amperes at 8.6 volts to heat it to 1,000° 
to 1,100° C, a satisfactory temperature range. Thin- 
ner filaments should not be used because they are 
apt to burn out during a run. 

A run is started by flushing the entire apparatus 


Ber. deut. chem. Ges. 42, 3200 (1909); see also, Stahler and Bachran, Ber. deut 
hem. Ges. 44 2006 (1911 





Figure l. Apparatus for the preparation of anhydrous titanium 
(111) chloride and titanium (111) bromide. 


A, Pyrex reaction kettle; B, central standard taper joint of reaction kettle; C, 
chimney; D, rubber stopper; E, 6-mm tungsten rods; F, 1-mm tungsten filament 
G, hydrogen inlet; H, initial hydrogen outlet; J, final hydrogen outlets; K, 
burette 


with dry oxygen-free hydrogen, purified by passing 
it over hot copper turnings and then over a desiccant 
such as magnesium perchlorate. After purification 
the hydrogen should pass through copper tubing 
or, if necessary, a minimum length of sound rubber 
tubing. The purified hydrogen is introduced at G 
and passes through the distillation assembly and out 
through the stopcock, H, until all the air has been 
removed. This will also remove moisture and thus 
prevent the formation of hydrolysis products. 

The hydrogen stream is continued while stopcock 
H is closed and that at J is opened. About 1,500 
ml of titanium (IV) chloride is run into the distilling 
flask through the burette, K, and distillation is 
started. About 1,000 to 1,500 ml is distilled over 
into the reaction kettle. The distilling apparatus 
ean then be removed from the reaction kettle and 
the opening, L, quickly stoppered. The hydrogen 
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stream should be very rapid during this operation. 

A hemispherical heating mantle is now used to 
reflux the tetrachloride over the tungsten filament. 
A heating mantle that entirely encloses the reaction 
kettle must not be used, as this will cause the tetra- 
chloride to reflux too rapidly and flood the condenser. 
This situation is detrimental because, after the re- 
duction is started, solid trichloride will be carried out 
of the apparatus by the hydrogen stream. Most, if 
not all, of the tetrachloride should be condensed on 
the walls and top of the reaction flask. Jets of air 
directed against the upper part of the flask will aid 
in this. 

When the boiling of the tetrachloride has reached 
equilibrium, the filament is turned on and brought 
to a bright red heat. Reduction begins with the 
formation of a dark purple smoke of titanium (IIT) 
chloride which collects on the walls of the reaction 
vessel. Most of these particles are washed down to 
the bottom by the refluxing tetrachloride. 

If for any reason the vessel is cooled, for instance 
by turning off the filament, the hydrogen flow must 
be momentarily increased to prevent influx of air as 
a result of the partial vacuum formed. This influx 
would cause hydrolysis and oxidation. A length 
of rubber tubing on outlet J minimizes this danger. 

The reaction is continued until the tetrachloride no 
longer refluxes freely down the sides of the vessel. 
The filament is then turned off, and the whole reac- 
tion kettle is heated to drive off the residual tetra- 
chloride, which tends to be absorbed in the cake of 
solid trichloride. To and remove the 
tetrachloride, a still head and horizontal condenser 
must be substituted for the reflux condenser. 


condense 


The hydrogen stream is continued while the appa- 
ratus cools. The flask should be opened and the 
product bottled in an inert atmosphere. Titanium 
(IIL) chloride thus formed often ignites spontaneously 
when exposed to moist air and invariably does so 
when gently heated. 

In a typical run, | liter (1,700 g) of titanium (IV) 
chloride was distilled into the reaction flask and 
5 ielded 154 g of titanium (I11) chloride, a conversion 
of 10 percent. Based on the tetrachloride consumed, 
the yield of trichloride was about 90 percent. 


Analysis of the product gave 


Ti Cl |Reducing powe 


31.9 
31.0 


67. 6 
68. 9 


97. 2, 98. 7 


100. 0 


Found (percent) 
Calculated (percent) 


The reducing power was determined by titratio 
with potassium dichromate. The crude produ 
was heated in a vacuum at 200° C, and a sma 
quantity of a volatile impurity was driven off 
After this treatment a reduction value of 101 
percent was obtained, indicating the presence of ; 
small amount of titanium (II) chloride. 

It is possible to obtain a much higher percentag 
conversion, particularly if the whole flask is heated 
but two factors must be borne in mind. First, 
solid reaction mass at the bottom is no longer being 
stirred by its own ebullition, and second, the solid 
cake on the sides, if no longer washed down by the 
refluxing tetrachloride, tends to build out toward the 
filament. Either of these factors could result in 
localized overheating and consequent disproportion- 
ation of trichloride into dichloride and tetrachloride 

The rate of formation of titanium (III) chloride 
is largely determined by the size of the filament 
This limitation could be circumvented if a sheet of 
tungsten, preferably in the shape of a cylinder, 
were heated by induction. 

Titanium (III) bromide has been prepared in th 
same apparatus. As the tetrabromide is solid at 
room temperature, hot water or steam is used as a 
coolant for the condensers (instead of the cold water 
used in the preparation of trichloride). The tetra- 
bromide is poured hot into the distilling flask and 
allowed to solidify before flushing with hydrogen 
The yield and purity of the product are very simila: 
to those reported for the trichloride. 

In some runs a multiple collecting device was used 
which made it possible to discard the first fraction of 
distilled tetrachloride without opening the apparatus 
It has not been shown that this is either necessary 
or beneficial. 


WasHINGTON, November 17, 1950 
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An Apparatus for Studying Autoignition of Engine Fuels: 
Results with Normal Heptane and Normal Hexane 


William J. Levedahl and Frank L. Howard 


\ single-cylinder AST M—CFR fuel-testing engine with a variable compression ratio has 
been modified to facilitate the study of compression-ignition of homogeneous fuel-air mixtures. 
An investigation of this type of combustion may lead to a better understanding of engine 
knock. A cylinder-pressure limiter is described that permits firing a small number of cycles, 
the first of which is uncontaminated by exhaust gases. Pressure-crank angle oscillograms 
of the process are shown. The combustion of n-heptane and n-hexane was found to take 
place in two distinct stages over a wide range of fuel-air ratios. The pressure at which the 
first-stage reaction begins is relatively unaffected by the mixture composition throughout 
the firing range. The second, or firing, stage occurs most readily at fuel-air ratios near 
stoichiometric, and ionization during firing is highest in this range of mixture compositions. 
No ionization was observed during the first stage of combustion with either fuel. 


I. Introduction 


Many investigators have recognized two major 
stages in the combustion of hydrocarbons. The first 
stage is characterized by a relatively slow reaction 
following an induction period whose duration depends 
on experimental conditions such as temperature, 
pressure, fuel-air ratio, presence of additives, and the 
condition and area of the surface exposed to the mix- 
ture. Cold flames, blue luminescence, and relatively 
small changes in pressure are some of the physical 
phenomena associated with the first-stage reactions. 

From a chemical standpoint the first isolable ma- 
terial from this reaction is peroxidic in character, and 
hydroperoxides have been identified. The exact 
mechanism of peroxide formation is unknown, but 
it is presumed to result first from the formation of 
short-lived free radicals from the hydrocarbon mole- 
cules, which in turn react with oxygen to form perox- 
ides. Decomposition of these hydroperoxides gives 
rise to new oxygenated free radicals that act as 
carriers for the propagation of a chain reaction, some 
of whose intermediate products have been identified 
as aldehydes, ketones, and acids. 

The progress of the first-stage reaction is governed 
by experimental conditions and by the molecular 
configuration of the fuel. It is unlikely that it ever 
reaches a condition of equilibrium in an engine, or 
even that a concept of equilibrium could be applied 
to the reaction. Possibly the first-stage reaction in- 
fluences the second or “firing’’ stage of combustion 
by providing certain essential chain carriers. It may 
be that the second stage will not occur until the 
first stage has progressed to a certain point. 

In a spark-ignition engine, normal combustion is 
initiated at the spark plug, and a flame front travels 
through the mixture from this point of origin, com- 
pressing the unburned mixture of fuel and air ahead 
of it. This unburned portion of the charge is fre- 
quently called the “end gas.’’ The temperature of 
the end gas is increased due to adiabatic compression 
by the burning gases in and behind the flame front 
and by motion of the piston. This temperature rise, 





accompanied by increases of pressure and density, 
sometimes causes first-stage reactions to occur, often 
followed by a very rapid second-stage reaction, which 
may cause combustion knock if it takes place with 
sufficient rapidity. It is thought that investigation 
of the characteristics of these reactions may prove 
helpful in obtaining a better understanding of com- 
bustion knock. 

In 1947-48 Levedahl and Sargent [1] ' studied the 
autoignition characteristics of diethyl ether, n-hep- 
tane, and other fuels in a single-cylinder CFR ? 
engine in which a well-mixed charge was ignited by 
compression in the cylinder. Diagrams of this 
indicated pressure versus crank angle showed very 
marked and reproducible evidence of a two-stage 
combustion reaction. Diagrams made at lower 
compression ratios showed results similar to those 
obtained by Peletier [2] and Muhlner [3] in that 
only the first stage occurred. Peletier made a 
number of runs at compression ratios low enough to 
preclude firing and measured the extent of the first 
stage reaction by measuring the temperature of the 
exhaust gas. Mulner went further with this method, 
making graphs of the rise of temperature and the 
torque obtained at many fuel-air ratios and compres- 
sion ratios with several gasolines, and obtaining some 
pressure diagrams. 

Retailliau, Ricards, and Jones [4] extended this 
technique by making analyses of the exhaust-gas 
and other refinements. Pastell [5] carried out ex- 
periments with the same general type of apparatus 
and added a photomultiplier tube to detect the radia- 
tion from cool-flames during the (first-stage) re- 
actions. He also measured the energy liberated by 
these reactions. 

A number of single-cycle rapid-compression ma- 
chines [6 to 15] have been built for the purpose of 
studying autoignition. Most productive of these 
has been the MIT rapid compression machine with 
which Leary and coworkers [12 to 15] have obtained 
significant results from measurements of the delay 

! Figures in brackets indicate the literature references at the end of this paper 


?The Motor Fuels Division of the Cooperative Fuel Research Committee 
develope i the CFR engine for knock testing of engine fuels. 
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time in the ignition of several fuels under 


many 
conditions. 


In order to extend and augment work previously 
reported on compression ignition studies, a program 
has been initiated to determine how the structure 
of the fuel affects the duration of preflame reactions 
and the pressure developed by these reactions. It 
is planned to use several six-carbon hydrocarbons 
as fuels. These compounds have a variety of struc- 
tures; paraffins, naphthenes, olefins, and an 
matic being included. This paper describes the 
apparatus and the methods of operation, and 
presents the results of preliminary experiments in 
which n-heptane and n-hexane were used as fuels. 


aro- 


Il. Apparatus 


Engine Figure 1 shows the ASTM supercharge 
method fuel-rating engine |16], modified for the 
present work, together with the accessory equipment 
used. A special cylinder was obtained with three 
holes for 18-mm spark plugs around its circumference 
and one %-in. hole in the top. The springs beneath 
the cylinder were removed and the zero setting 
on the micrometer used to measure the cylinder 
height was changed 0.3 in. so that compression 
ratios in excess of 18:1 could be obtained. An 
alternating-current dynamometer was used to main- 
tain the engine speed at 600, 900, 1,200, or 1,800 
rpm. 


Cooling system. Immediately before entering the 


Figure 1. Engine and auxiliary equipment. 
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cylinder jacket, the cooling water passed through t! 
central tube of a heat exchanger consisting of \-i 
copper tubing within a %-in. tube. Either stea: 
or water could be passed through the annular spac 
In this way the temperature of the cylinder jack: 
could be maintained constant at any desired ten 
perature between 150° and 212° F when the engi: 
was in operation under either motoring or firi: 
conditions. A motor-driven vane pump circulat 
the cooling water through the jacket from a coolar 
tower (A, fig. 1) used as a reservoir. The exhaus 
system was cooled by passing water through coi! 
wound on the exhaust pipe, B, and by water spraye: 
directly into the gases. 

Induction system. Air for the mixture was sup 
plied at high pressure, passed through a water tray 
and was then throttled to 40 psig. <A pair of surg 
tanks, C, and a thermostatically controlled air heate: 
insured a steady flow of air through an ASME sharp 
edged orifice plate, D. The drop in pressure across 
the orifice was read directly on a manometer, E 
calibrated in pounds of air per hour. Following 
another throttle valve, which automatically main- 
tained inlet pressure at the desired level, the air 
passed over a 3,000-w electric heater, F. Fuel was 
injected into the airstream as it left the heater ; vapori- 
zation and thorough mixing took place in a large 
steam-jacketed mixing tank, G. The charge then 
flowed over a thermometer and through a copper- 
tube flame trap into the engine. The fuel was 
metered by one of two calibrated flowmeters, H 
with a combined range of approximately 0.5 to 20 
lb of fuel per hour, and the rate was controlled 
manually with a needle valve. The composition 
of the fuel-air mixture is expressed throughout this 
paper as fraction of stoichiometric fuel, and is given 
the symbol 5 in the figures. 

Pressure-measuring equipment. In order to obtain 
records of the variation of pressure with time during 
combustion cycles it was necessary to use a sensitiv: 
and sturdy pressure pickup with a rapid response and 
a stable calibration. A Draper-Li strain-gauge unit, 
[, with an extra-heavy catenary diaphragm and 
strain-generating tube was selected. A direct-current 
bridge circuit, of which the pickup formed two legs, 
and an electronic amplifier fed a voltage that varies 
linearly with cylinder pressure to an oscilloscope, J 
with direct-current amplification on the deflection 
plates. A resistor in parallel with one leg of the 
bridge was automatically cut into the circuit once 
each cycle by a microswitch and caused a calibration 
mark to be made on the screen at a known pressuré 
As the indicator operated at a voltage well above 
ground, a 22':-v battery-potentiometer circuit was 
used to put a steady opposited bias on the y-axis 
input, and to permit use of the direct-current ampli 
fier on the oscilloscope. 

Tonization-detecting equipment. Ionization of the 
gases in a spark ignition engine has been observed in 
both knocking and nonknocking combustion [17,18] 
In order to investigate the relationship betwee 
ionization and time in autoignition cycles, a standard 
18-mm spark plug, K, was used as the ionization ga] 
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in a circuit formed by a 135-v battery and a 250,000- 
ohm resistor 
was proportional to the ionization current and was 
fed to the vertical amplifier of the oscilloscope. 

Timing and recording equipment. In order to 
correlate the pressure-time and ionization-time traces 
on the oscilloscope with crank angle, it was desirable 
to have the beam sweep across the face of the tube 
once each revolution and have the trace marked to 
indicate known crank positions. These two fune- 
trons were performed by a photoelectric cell, L, which 
reacted to a beam of light reflected from fine strips 
of tinfoil located on the black flywheel so as to register 
top center, 30°, 60°, and 90° before top center, The 
strip at 90° before top center triggered the sweep, 
and the others produced small vertical marks. The 
oscilloscope traces were photographed with a camera, 
\I, whose exposure time was set to expose the desired 
number of consecutive cycles. The shutter was 
tripped at a selected time in the cycle by a micro- 
switch and solenoid actuated by a projection on the 
engine camshaft 


Cylinder-pressure limiter. It is often desirable to 


FicuRE 2 ( jlinder-pressur limiter 

record combustion cycles in the absence of residual 
exhaust gases and hot spots in the combustion cham- 
ber. This can be done without disturbing the flow 
of fuel and air to the engine by keeping the pressure 
in all except a few chosen cycles from rising enough 
to cause appreciable chemical reaction. A pressure 
limiting device was developed for this purpose and 
mounted in the upper cylinder hole (N, fig. 1). It 
s shown in assembly and layout in figure 2. 

When the pressure 1n the evlinder exceeds 10 psig, 
the valve (A, in fig. 2) leaves its seat in the body, B, 
and moves upward until stopped by the piston, C 
The gases in the evlinder can then escape through 
the annular space between A and B, through holes 
to the exhaust collector, D, and thence to an exhaust 
On the power stroke, when the pressure 
again drops below 10 psig, the limiter valve, A, is 
losed by the spring, E, and remains closed through- 
out the exhaust and intake strokes, thus permitting 
the induction of a full charge of combustible mixture. 


system. 


The voltage drop across the resistor 


PRESSURE PSIA 


_ 


-30° Toc 30° 60° 
CRANK ANGLE 


Fr URE: Pressure crank-angle diagram made with pressure 


imiter open ata compression ratio of 6 59 


In this way the mixture can pass through the engine 
without being burned. 

When it is desired to fire the mixture the switch, 
F, is closed, and at top center on the next exhaust 
stroke a microswitch actuated by the engine camshaft 
completes the circuit to the solenoid, G, which pulls 
the trigger, H. The lever, I, which is attached to 
the camshaft, J, is released, permitting the torsion 
spring, K, to rotate the cam against the piston, C 
After the cam has rotated about 60° the piston is 
locked firmly against the valve, A, by the cam, 
permitting the pressure in the cylinder to rise to that 
of full compression. 

When the lever, I, rotates, it operates the micro- 
switch, L, and when bottom center on the inlet stroke 
is reached another microswitch is closed by the cam- 
shaft to complete the circuit to the recording camera. 
Thus the film is exposed during the first cycle after 
closing the pressure limiter and for subsequent cycles 
until the shutter is closed. Figure 3 is a typical 
pressure record made at 900 rpm with a compression 
ratio of 6.6 and with the pressure limited in the man- 
ner just described. The pressure reached a maximum 
of about 47 psig at approximately 20° before top 
center and dropped off rapidly thereafter. 


Ill. Experimental Procedure 


Steady firing conditions. The vaporizing tank was 
preheated by steam, and the engine was motored 
(driven by the dynamometer) until the inlet air 
temperature reached the desired value. With the 
compression ratio set at about 8:1, the fuel was 
introduced at a rate to give mixtures from 0.5 to 0.7 
of stoichiometric. Firing usually started imme- 
diately; if not, the compression ratio was in- 
creased until it did occur When the rates of flow 
of the fuel and air and the temperatures had been 
adjusted to the desired values, the oscilloscope trace 
of one cycle was photographed. 

Cylinder-pressure limiter. The engine was mo- 
tored at full compression without fuel, and hot 
water was circulated through the cooling jacket 
until the temperatures were stabilized. The pressure 
limiter was placed in operation, the desired flow of 
fuel was established, and all temperatures were 
again allowe? to stabilize. The exposure time 
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of the camera was set to expose three cycles, and | temperature equilibrium is virtually undisturb 
the limiter switch, F, was closed. The limiter lever | bv three firings, series of observations can be ma 
was raised manually after three firings. As the | mm rapid succession 
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PRESSURE IN POUNDS PER SQUARE INCH ABSOLUTE 





: 
_ 
7 . 60° 
Ficure 4 Pressure-crank angle diagrams of n-heptane combustion with several mixture compositions 
at 600 rpm, 220° F inlet te m perature, 180 F jacket tem pe rature, a consumption 18 tb hr. and 
compression ratio ,.88 
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IV. Results and Discussion 
1. Steady Firing of n-Heptane 


Effect of mixture composition. The engine was 
operated on n-heptane over a wide range of mixtures 
and compression ratios under the following condi- 
tions: inlet temperature, 220° F; jacket temperature, 
180° F; air rate 18 lb/hr; speed, 600 rpm. At a 
compression ratio of 7:1, the engine ran steadily at 
all fuel mixtures between 0.3 and 2.5 times stoichio- 
metric. No knocking was observed with the leanest 
mixtures, but all those richer than about 0.5 of 
stoichiometric knocked violently. Decreasing the 
compression ratio reduced both the knock intensity 
and the range of fuel-air ratios that weuld fire. 
Figure 4 shows one set of typical diagrams of pressure 


versus crank angle taken at a compression ratio of | 


5.68 and at various fuel compositions in the range 
0.49 to 1.76. The occurrence of two reactions, as 
indicated by sudden changes in the slope of the 
pressure curve, are first evident in figure 4, B. The 
coordinates of pressure and crank angle at which 
these changes in slope occur are designated by P, 
and @, and by P, and 6 for the first and second 
stages, respectively, in figure 4, C. 

Firing was steady and reproducible from cycle to 
cycle at all mixture ratios. In figure 4, A, only the 
first-stage reaction took place. In figure 4, B, the 
mixture was at the “critically lean’’ point, that is, a 
very small decrease in the fuel concentration would 
prevent the second- or firing-stage of the reaction 
from occurring. Conversely, a very small increase 
in fuel would cause steady firing. 

Observations of pressure and crank angle made 


from figure 4 are plotted as functions of mixture 


composition in figures 5 and 6, respectively. Figure 
5 shows that the pressure, P), necessary to start the 
first-stage reaction is nearly constant over a wide 
range of mixtures above the critically lean value of 
0.52. The pressure at which the second stage begins 
is lowest near stoichiometric and highest for lean 
mixtures. Figure 6 shows that the time 6, (in terms 
of crank angle) at which reaction begins is nearly 
constant, except for very lean and very rich mix- 
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Figure 5. Plot of data from figure 4: Effect of mixture com- 
position on pressures for beginning of first (P,) and second 
\ P)) stage 8 of combustion. 


Mixture composition is expressed as fraction of stoichiometric fuel 
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Ficure 6. Plot of data from figure 4: Effect of mixture com- 
position on crank angles at which the first stage (@,) and second 
stage (02) of combustion begin, and duration of the first-stage 
reaction (6,.—;). 


Mixture composition is expressed as fraction of stoichiometric fuel 


tures. It is possible that the reaction starts later in 
rich mixtures because of their lower specific heat 
ratio; the temperature and pressure required to 
start the reaction are not reached until later in the 
piston stroke. The crank angle at which firing 
occurs (@,) decreases sharply as the fuel-air ratio is 
increased in the very lean range, flattens out be- 
tween about 0.6 and 0.8 stoichiometric then drops 
more gradually as S is increased. The duration of 
the first-stage reaction (@,—6,) decreases sharply in 
the lean range, then decreases more gradually as the 
mixture becomes richer. 


2. Ionization as an Indication of Firing 


Diagrams of ionization versus crank angle. The 
engine was operated under several combinations of 
mixtures, engine speed, and inlet temperature to 
determine the relation between the ionization at one 
location in the cylinder and the corresponding pres- 
sure. Innocase was any ionization detected during the 
first-stage reaction, but a sharp rise in ionization was 
always noted exactly at the onset of the second stage 
of combustion. Typical ionization curves for ex- 
periments in which n-hexane was used are shown in 
figure 7. Figure 7, A, shows a low wide peak charac- 
teristic of lean mixtures, and figure 7, B, shows a 
curve typical of richer mixtures. Figure 7, C, shows 
three consecutive cveles, in each of which the sudden 
rise in ionization current occurred at the same crank 
angle. The maxima (not shown in this curve) show 
cycle-to-cycle variations as high as 20 percent. After 
passing the peak, the traces show considerable 
oscillation. 

Figure 8 shows the effect of mixture composition 
on peak ionization current, using n-hexane as fuel. 
In every case the compression ratio was set to cause 
the onset of ionization at exactly top center. When 
the fuel content of the mixture was increased from 
0.5 to 1.1 stoichiometric, the ionization increased 
about 200 times. Points taken at three engine 
speeds and two inlet temperatures lie very nearly 
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on the same curve, and data taken with n-heptane 


under comparable conditions have shown similat 


results. From this it would appear that peak ioni- 
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Mixture ens 
zation in autoignition is primarily a function of 
mixture composition, whereas the changes in other 
variables cause only secondary effects. 

Effect of m prture composition on the compre S810? 
ratio required to cause fi ing at top cente The com- 
pression ratios used to obtain the points in figure 8 
are plotted against mixture composition in figure 9 
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Ficure 9 Effect of mixture com position on compression 


req wired to 


maintain continuous firing al top cente 


n-hexane as fue 


The effect of engine speed on the required compres- 
sion ratio is slight, decreased inlet temperature con- 
sistently raises the requirement throughout the 
mixture range. The effects of fuel composition are 
great, however, with the minimum requirement being 
at a point slightly richer than stoichiometric. With 
mixtures between 0.6 and 0.8 stoichiometric there is a 
reversal of curvature that is comparable to the re- 
versal in the firmg angle curve (@,) of n-heptane in 
figure 6. The reasons for this are not yet understood 
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3. Experiments with n-hexane with the cylinder 
pressure limiter 


The cylinder-pressure limiter was installed for t 
rest of the work with n-hexane. Figures 10 and 
show diagrams of pressure versus crank angle ma 
at 900 rpm, 212° F jacket temperature, 220° F in! 
temperature, and 30-lb/hr air consumption. A 
each mixture the first set of observations was ma 
at the compression ratio that caused firing at to 
center under steady-firing conditions (see fig. 9 
Subsequent observations were made at progressive! 
higher compression ratios, and each was repeated t 
determine reproducibility. Average  differenc: 
amounted to less than 1° of crank angle and 1 psi fo 
6, and P,, respectively. An accurate determi 
nation of the pressure at which firing begins is difficul 
because the slope of the curve at that point is steey 
enough to cause considerable reading error. 

The oscillograms in figure 10 (A, B, and C) wer 
made at progressively increased compression ratios 
with a mixture of 0.5 stoichiometric. The superposed 
pressure curves of the secondand third cycles show 
earlier firing than in the first cycle at all three compres 
sion ratios. This fact seems to indicate that residual 
exhaust gases from one cycle had an accelerating 
effect on combustion in the next cycle. At a richer 
mixture the first cyele in figure 10, D, showed no 
perceptible reaction, but the second and third cycles 
fired. At higher compression ratios (fig. 10, E, and F 
all three cycles fired. 


INCH ABSOLUTE 


Figure 11, A, shows a slight first-stage reaction 
beginning at about 10° after top center on the third 
cycle with a stoichiometric mixture. Several cycles 
later (not shown in the figure) the engine began to 
fire. When the compression ratio was raised (fig 
11, B), no firing occurred on the first cycle, but the 
second and third cycles fired progressively earlier 
It appears that the residual gases from a motored 
cycle may have an accelerating effect on the reactions 
in the next cycle even though no perceptible reaction 
took place in the first 
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With a very rich mixture, 1.5 times stoichiometric 
figure 11, C, shows that firing began very late on the 
second cycle and near top center on the third. Only 
a slight evidence of a first-stage reaction is apparent; 
in the third cycle it appears as a fillet (shown by 
arrow) between the compression and firing parts of 
the pressure trace. 

The pressure rise due to first-stage reactions 
decreases as the mixture becomes richer; with th 
leaner mixtures it is apparently necessary for the 
first stage to progress farther before the second stag« 

a can begin. At the higher fuel concentration it was 

Toc 30° necessary to use much higher compression ratios to 

cause firing when the engine is being motored than 

CRANK ANGLE those shown in figure 9 for steady firing conditions 

This phenomenon may be attributed to the larg: 

difference between the temperatures in a motoring 
and a heavily knocking engine 
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Ficurse 11. x re ( angle oscillograms at richer 
miziure. 
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V. Conclusions 


n-heptane and n-hexane both evidence a pro- 
nounced two-stage autoignition reaction over a 
wide range of operating conditions. 

2. When all other conditions are held constant, 
the pressure and crank angle at which the first-stage 
reaction begins in a steadily firing engine are nearly 
constant over a range of mixture compositions from 
critically lean to 1.2 times stoichiometric. 

The duration of the first stage of the reaction 
decreases as the mixture is made richer up to about 
1.2 stoichiometric, then remains nearly constant. 

4. Considerable ionization occurs in the cylinder 
during firing, but none was observed during the first- 
stage reaction. The ion concentration is affected 
greatly by fuel-air ratio and is relatively unaffected 
by changing other operating variables. 

When residual gases from either a first- or 
second-stage combustion are in the cylinder, the 
charge will burn earlier than when uncontaminated. 


The authors are grateful to E. F. Fiock, Chief of 
the Combustion Section, for his advice and assistance 
in preparing the manuscript of this report. 
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Creep of High-Purity Aluminum 
William D. Jenkins 


A study was made of cold-dra 
tions of creep rate 


-n aluminum at 105° F, 
on the mechanism of deformation. 
crease monotonically with increase in strain during first stage of creep 


to determine the effects of varia- 
Strain-hardening was found to de- 
The extension on 


loading and the extension at the beginning of the third stage both increased with increase in 


stress. 


Cyclic temperature changes appeared to increase the ductility. 


Discontinuous flow 


was evidenced in constant-load creep tests at 105° F and in tensile tests conducted at 80° F. 


Conformance to 
rates 
creep rate 


I. Introduction 


Many empirical formulas have been developed and 
used with rather limited success to evaluate the re- 
lationships between temperature, stress, and creep 
rate of metals and alloys. Interpretations of these 
relations are complicated by the various types of 
flow occurring during the first, second, and third 
stages of creep. Some investigators, however, agree 
that movement can occur both at the grain bounda- 
ries and within the grains in creep, and that the 
parent grains break up into crystallites whose dimen- 
sions are dependent on such factors as chemical 
composition, temperature, strain rate, and the stage 
of the creep test [1, 2].' 

Discontinuous flow has been observed in 
tests made at this Bureau on annealed high-purity 
copper [3] and on cold-drawn ingot iron [4]. Further 
evidence of this phenomenon, observed in the tests 
with high-purity aluminum, will be discussed in this 
paper. 


creep 


II. Material and Testing Procedure 


All specimens were cut from a i-in. round bar of 
high-purity aluminum that had been cold drawn to 
10-percent reduction of area. This aluminum con- 
tained the following percentages of other elements as 
determined spectrochemically 


0.012 My 0. 004 
Ol! ri 003 
007 Mi 002 


Some properties, at room temperature F), of 


the cold-drawn aluminum are as follows: 





000 
500 


Yield strength, 0 2-perce nt offs« :. psi l 
Tensile strength, ps! l 
Elongation in 2 in., 


1, 
7, 
percentat maximum load 

Elongation in 2in., percent at fracture 

Reduction of area, percent at maximum load 

Brinell oumber, 500-kg@ load; 10-mm-di- 
ameter ball 


Andrade’s law of transient flow was found within a limited range of strain 
Sigmoidal curves were produced that depicted the relationship between stress and 
\ mechanism of flow based on observed structural changes is proposed. 


The creep tests were made on round specimet 
having a gage length of 4 in. and a gage diameter o! 
0.505 in. These tests were carried out in air at 
temperature of 105° F, with the exception of one test 
for which the temperature was changed from 105 
to 95°F and back to 105° F. The specimens wer 
held at the test temperature for 48 hr before loading 
and the rate of loading was controlled so that a: 
average rate of extension of 30 to 50 percent per 
1,000 hr was maintained until the desired creep load 
was obtained. Thereafter the load was not changed 
until the test was discontinued. The creep-testing 
equipment has been described in a previous paper [5 

Tension tests also were made at room temperaturs 
on specimens similar to those used for the creep tests 
and true stress-strain values were determined from 
simultaneous load and diameter measurements mac: 
during these tests. Usual procedures were used in 
preparing the specimens and in carrying out the 
metallographic examinations. 


III. Results and Discussion 


l. Influence of Stress on Plastic Extension During 
Loading 


A composite stress-extension curve for different 
specimens loaded at rates of extension of 30 to 50 
percent per 1,000 hr is shown in figure 1. The 
homogeneity of the material is indicated by the lac! 
of appreciable scatter of the observed values With 
this rate of loading the change of strain hardening 
with strain decreases continuously with increase 1 
These conditions are evidence of the fact 
that the course of the reaction to a more stable stat: 
is speeded as the stress is increased. Similarity to 
chemical-rate reactions should be noted as the ability 
of metals to flow is increased by increase in strain 
energy. The subsequent work-hardening characte! 
istics are described in succeeding sections. 


stress. 


2. Extension-Time Characteristics as Influenced by 
Stress 


The general relation between extension and tim 
during creep of metals and alloys tested in tension a 
a constant load and a constant temperature can b 
simply represented by an extension-time curve. Fo 
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TEMPERATURE 105° F 

RATE OF EXTENSION DURING 
LOADING-30 TO 50 PERCENT 
PER 1000 HR 


NOMINAL STRESS,1I000 PSI 











0.05 0. O15 0.20 
PLASTIC EXTENSION, PERCENT 
} IGt " Efe 4 of nominal stress on plastic extension du? ing 
loading. 


Each symbol represents values obtained in individual tests 


ideal conditions, the course of such a curve is divided 
into parts as follows: (1) An initial extension (elastic 
and plastic), which occurs on the application of load, 
2) astage of decreasing rate of extension (designa- 
ted as primary or first stage of creep), (3) a stage of 
nearly constant rate that attains a minimum (desig- 
nated as secondary, minimum rate, or second stage 
of creep), and (4) a stage of increasing rate even- 
tually leading to complete fracture (designated as 
tertiary or third stage of creep). 


The family of extension-time curves for various 


applied stresses are shown in figure 2. A general ex- 
amination of these curves shows that, as the applied 
stress is increased, the time for initiation of each of 
the stages is decreased, and the initial extension on 
loading and the extension at the beginning of the 
third stage are increased. The variation of stress 
with time to attain some equal value of strain appar- 
ently conform to a power relationship; however, no 
simple qualitative relationship could be established 
between stress and extension at the beginning of the 
second stage nor the intercept extension produced by 
extrapolating the second stage extension-time curve 
to zero time due to the existing work-hardening and 
recovery characteristics of the material various 
stress levels. 

A suggested partial explanation for these phenom- 
ena is that the deformation occurs in part by the 
breaking up of some of the parent crystals into smaller 
units designated as subcrystals [3], which for a fixed 
temperature are smaller, the greater the initially ap- 
plied stress. The factors that produce the smaller 


units, tend to cause greater hardness and thus in- 
crease the rate of strain-hardening; however, due to 
the randomness of the adjacent sube rystals (figs., 10, 
A and B) the tendency to recover is more apparent 
than in larger crystals, and thus greater duc tility at 
the beginning of the third stage is ‘obtained. At this 
point the atoms are in a more stable position than at 
any previous strain state, and the increase in free 
energy of the system with time is ata minimum. The 
shape of the extension-time curves is further in- 
fluenced by local imperfections due to plastic de- 
formation in addition to the inherent structural 
stability of the material and the distribution of the 
parent grains that did not fragment during deforma- 
tion. These factors account in part for the transient 
effects observed when some of the test conditions were 
altered and then returned to their original state. 


3. Effect of Cyclic Temperature Changes on Creep 
Characteristics 


The extension-time and creep-rate time diagrams 
for a specimen tested in creep with a stress of 13,200 
psi (consecutively at temperatures of 105°, 95° and 
105° F) are shown in figure 3. The average creep 
rate before lowering the temperature was 0.37 percent 
per 1,000 hr. The change to the lower temperature 
did not cause an immediate steady-state rate condi- 
tion, but rather a gradual reduction in the slope of 
the extension-time curve. This indicates the presence 
of a first stage under these conditions. After the 
attainment of approximately constant rate at 95° F, 
the subsequent creep rate on raising the temperature 
to 105° F appeared to be unaffected by previous test 
temperature. These findings confirm the results 
presented in a previous investigation with copper [3], 
which showed that other disturbances could produce 
entirely different results. The conclusion in such 
cases was that the subsequent creep rate after tem- 
perature change depends essentially on the prior his- 
tory of the material, despite the fact that changes in 
test conditions may be expected to alter the flow 
mechanism. As would be predicted, the creep rate- 
time cycles, illustrated in the upper portion of figure 














Extension-time curves for specimens teste d in cree P 
eith diffe rent init ally applied loads. 
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3, are more prominent in the range of lower tempera- 
ture. The deformation at lower temperature and the 
thermal gradients set up as a result of the tempera- 
ture change altered not only the deformation process 
as reflected in the strain rate-time diagram, but also 
caused an increase in extension at the beginning of the 
third stage as was shown in figure 2. Thus the range 
at which the free energy of the process is at a mini- 
mum is shifted to a higher extension value. Con- 
cepts of this type must be accounted for in physical 
interpretations of the deformation process, as it has 
been shown previously [3, 4] that neither recrystalli- 
zation on a macro scale nor increase in average true 
stress is necessarily a prerequisite for initiation of the 
third stage of creep 
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Fiaure 3. Extension-time and creep rate-time relations for a 


specimen le sted with constant nominal stress af 13,200 psi, and 


and 95 F. 


tem perature cycled between 105 


4. Variation of Creep Rate With Plastic Extension 


The creep rate-plastic extension curves for two 
specimens tested with different initial stresses are 
illustrated in figure 4. The values are plotted on 
a semilog basis. It is evident that the magnitude 
of the creep rate-strain cycles, when the data are 
plotted in this manner, is greater in the specimen 
strained at the lower (extension readings 
were made at equivalent time intervals). Further- 
more, as is shown in the upper curve, after an initial 
adjustment period, the amplitude of the cycles 
becomes until third stage conditions are ap- 
proached. The strain-hardening process tends to 


stress 


less 
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Ficure 4 Variation in creep rate with plastic 

cause an increase and the recovery a decrease in 
the activation energy necessary for flow. Thus th 
time for attainment of each stage of creep is shorte: 
as the applied stress is increased, and the deforma 
tion produced as a result of the metastable stat 
adjusting itself toward a condition of lower fre 
energy is greater. (For example, see fig. 2). In 
order for a decrease in free energy to occur, this 
system must undergo considerable atomic rearrange- 
ments. The mechanics of this process are illustrated 
by the serrations in these curves (fig. 4) and con- 
firmed by subsequent metallographic examinations 
The activation energy of the process is supplied in 
part by heat gain, the externally applied stress and 
localized internal thermal fluctuations. The cycling 
is also due in part to nucleation and growth of 
fracture sites. 

Log-log creep rate-time curves for the intial stage 
of creep of two specimens tested with nominal stresses 
of 12,700 and 13,600 psi are shown in figure 5. It is 
evident that both curves have approximately the 
same slope although at different levels. Such an 
equivalence tends to show that when the applied 
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Ficure 5. Variation in creep rate with time during first stage 


of creep, 
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stresses are of the same order of magnitude, the 
mechanism of deformation is similar during the 
initial stage. Andrade [6], in his analysis of first 
stage conditions, concluded that with this type of 
plot, the curves should have a slope of minus % when 
transient flow was the predominating characteristic. 
These two curves have a slope of about minus %, al- 
though such good agreement could not be obtained 
for the specimens at other stresses. Until the many 
types of flow have been clearly delineated, the appli- 
cation of these physical theories should be used 
with caution. Several competing and supporting 
mechanisms, dependent on the material, grain and 
ervstallite size, total deformation and prior strain 
history, may be operating simultaneously. 
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Ficure 6. Relationship between nominal stress and average 


creep rate during second stage of creep. 


5. Effect of Stress on Creep Rate in the Second Stage 


The relationship between the applied nominal 
stress and the average creep rate during the second 
stage is illustrated in figure 6. Sigmoidal curves are 
produced when these values are plotted on either a 
log-log or semilog basis. It seems probable that 
either type of representation reflects the flow charac- 
teristics of the material. Only in those regions of 
stress in which the mechanism of deformation and 
the subsequent modes of fracture, transcrystalline or 
intercrystalline, are nearly equivalent, would linear 
relationships between the variables of stress and 
creep rate exist. Previously, it has been shown [3] 
that the size of the subcrystals produced as a result 
of plastic deformation vary inversely with the creep 
rate, and the tendency toward transcrystalline 
separation of the parent grains is greater the faster 
the speed of the process. Transcrystalline fracture 
predominates above the point _of reversal in the 
stress-strain rate curve. 


6. Stress-Strain Relationships at Room Temperature 


McReynolds [7] in an analysis of stress-strain 
curves of aluminum tested in tension noted that 
macroscopic measurements almost invariably show 
a smooth curve of stress versus strain. Such 
homogeneity would indicate that discontinuities in 
strain on slip planes occur in increments too small 
to be recorded individually and in some manner 
independent of one another. However, as illustrated 
by a portion of the autographic stress-strain curve 
(fig. 7), inhomogeneity of deformation is reflected by 
the serrations (S) in this diagram for a tensile test 
conducted at 80° F with a rate of extension of 60.000 
percent per 1,000 hr. This test was interrupted at 
various extensions and the load dropped to zero in 
order to evaluate the effects of the healing of the 
internal surfaces previously ruptured due to strain, 
on the subsequent work-hardening characteristics of 
the material. The test specimen was held at zero 
load during each interval for 10 mins before reloading 
at the same rate at which the specimen was originally 
loaded. The consequence of the ability of the atoms 
to rearrange themselves into a more stable system 
during the rest period was reflected after reappli- 
cation of load in the production of a “kink” (K, 
fig. 7) in the stress-strain curve. This phenomenon 
was observed and described previously by Farren 
and Taylor [8]. After the appearance of the kink, 
there appears to bea lag in stress before a continuous 
increase in rate of change of load with extension is 
noted (fig. 7). This lag extends over greater exten- 
sion values the less the strain to which the specimen 
was subjected before the rest period. The rate of 
change of load with extension in the last part of the 
test appears to be independent of the load; however, 
multiaxial stresses due to the formation of a “neck” 
make further analysis difficult. 

Whereas the above observations were based on 
nominal stress values, the true stress-strain curve 
reveals other interesting results, as is illustrated by 
figure 8. The values are shown only to a true strain 
of about 0.9 due to the fact that beyond this strain, 
the specimen deformed with a quadrilateral contour 
at the base of the neck and microscopic cracks oc- 
curred at the axis of the specimen. These factors 
prevented quantitative determinations of the re- 
maining area of the material. In figure 8, A, the 
points at which the load was dropped to zero are 
indicated. The extension at maximum load was 1.4 
percent—a value that was slightly less than the ex- 
tensions obtained during the creep tests at the begin- 
ning of the third stage. It should be noted that the 
curve at 80° F becomes approximately horizontal 
for a limited range of true strain values immediately 
after maximum load. This indicates that the strain 
hardening in this strain range is zero. However, 
the contour of the curve between true strain values 
of 0.2 and 0.9 indicates that the strain-hardening 
increases slightly or becomes approximately constant. 
Some of these findings appear to confirm the conclu- 
sion of Farren and Taylor [8] that the ratio of the 
increase in internal energy to the work done is a 
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Enlarged scales were emploved for plot- 
ting the data in the regions in which the specimen 
was unloaded and subsequently reloaded. These 
curves are reproduced in figure 8, B, C, D, E, and F 
It is obvious that In spite of the occurrence of the 
kinks in the nominal stress-strain curve of figure 7, 
the yield stress on reloading was lowered In an 
analogous investigation, Cherian, Pietrokowsky, and 
Dorn [9] concluded that during the rest periods two 
In one type, flow 


constant 


types of recovery were evident 
values on reloading are below those of the 
virgin material for all values of total strain. Fast 
reattainment of flow stress was evident in the test 
conducted with the creep specimen at two different 
temperatures (fig. 3 Slow reattainment of the 
original flow stress Is apparent for the present tensile 
test (fig. 8 The change in extension at which the 
maximum flow stress is reattained after reloading is 
ereater the higher the strain values at which the 
specimen was unloaded Self recovery of the parent 
grains after deformation more evident as 
the strain is increased, 


stress 


be omes 


am tor ante iptled fens 


7. Microstructure and Fracture Characteristics 
The appearance of a fracture and some of the 
structural features accompanying deformation that 
led eventually to fracture of an aluminum specimen 
deformed at 80° F with a rate of extension of 60,000 
percent per 1,000 hr are shown in figures 9, 10, and 
11. Necking of the specimen commenced immedi- 
ately after maximum was reached, and the 
cross-sectional area at the base of the neck continu 
ously decreased to the end of the test. The appeal 
ance of the outer surfaces at fracture is shown in 
ficure 9, A. The general direction of granular elonga 
tion is outlined by the strain lines appearing at the 
surface. During the process of deformation the 
stress distribution was such that cracks (fig. 9, B 
were nucleated at the axis of the specimen and prop- 
agated to the surface. As this process of deforma- 
tion continued toward total fracture, the specimen 
finally exhibited a rim effect characterized by the 
extreme outer surface being the last portion to sep 
arate after considerable widening of the macro crack 
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Fiaure 8 True stress-strain relationships for a@ specimen 


stressed interruptedly at 80° F. 
/ / 


formed by the linking up of a number of microcracks 
in the interior. Figure 10, A, illustrates the struc- 
ture representative of the original cold-drawn alum- 
inum; the direction of the drawing is evident from 
the strain pattern. Some of the elongated grains 
appear to have broken up into small fragments due 
to the severity of the deformation, others are in- 
tact. Asa result of testing in tension at a fast rate, 
the mode of fracture is essentially transcrystalline, 
as illustrated by figure 10, B. Transcrystalline 
cracks that did not link up to form the main fracture 
course are shown at various points in this region. 
Figure 11 illustrates some of the etched patterns 
developed in areas near the axis of the specimen. 
The reduction of area in the region in which these 
photomicrographs were made was approximately 90 
percent. Figure 11, A, shows some of the grains 
elongated in the direction of the applied stress and 
containing suberystals. These substructures, ac- 
cording to Heidenreich {10}, are fully developed re 
orientation domains corresponding to self-recovery 
after cold work. The absence of fracture through 
the grains in which these substructures were ob- 
served in the present study suggest that these grains 
are in a rather stable state; however, conditions in 
all grains were not favorable for the development of 
this substructure. Figure 11, B, illustrates the etch 
patterns developed in some of the grains that were 
not resolved into suberystals. The strain markings 
as shown, indicate that during the course of deforma- 
tion a number of slip systems were active as would 
be predicted from the serrated stress-strain curves 


previously shown. The geometry and quantitative 


FIcuRE 9 Photograph of the necked portion of a specimen afte 
testing in tension at SU F 


tud view: BB. & t 


studies of the nucleation and growth of strain mark- 
ings of this nature have been made by Barrett [11], 
Banerjee [12], Rosi and Mathewson [13] and Yen 
and Hibbard [14]. 


IV. Summary 


Tension tests were made on cold-drawn aluminum 
at 105° F at various creep rates 

Plastic extension at the beginning of the third 
stage was greater the higher the rate of extension in 
the second stage 

Cyelic temperature changes appeared to affect the 
mechanism of deformation and ductility at the 
beginning of the third stage to a greater degree than 
the second-stage creep rate 

Discontinuous flow was indicated by graphs 
obtained by plotting the logarithm of the creep rate 
against extension occurring during the course of each 


creep test Discontinuous flow also was illustrated 
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by serrated stress-strain curves for a specimen tested 
at 80° F 

Within a limited range of strain rates a slope of 
about minus % obtained in first-stage creep 
curves when transient flow was the predominant 
characteristic that is essentially in agreement with 
Andrade’s conclusions. 

Sigmoidal stress-creep rate curves were obtained 
when the data were plotted on either a log-log or 
semilog basis. 

Yield point lowering as evidence of recovery at 
room temperature in an interrupted tensile test 
served to confirm the results of previous investiga- 
tions in this respect 

Self-recovery processes during deformation were 
illustrated by the breaking up of some of the parent 
grains into suberystals of microscopic dimensions 


was 


-? 


$ PEA nD Ses 
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F with a strain rate of 60,000 percent per 
fracture 


“IGURE 11. 
tension at SU 
hr; longitudinal sections, 0.1 in. from 
etched in 10 parts nitric acid (concentrated), 10 parts hydr« 
chloric acid (concentrated), 3 parts hydrofluoric acid (47% 
and 77? parts water. 
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and the occurrence of slip indicated by strain mark- 
ings on other grains. 
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I. Introduction 


The National Bureau of Standards is undertaking 
a systematic review of published X-ray powder dif- 
fraction patterns and an experimental investigation 
planned to revise and supplement them. The work 
is being performed in cooperation with the Joint 
Committee on Chemical Analysis by X-ray Diffrae- 
tion Methods, composed of members from the Ameri- 
ean society for Testing Materials, the American 
Crystallographic Association, and the British Insti- 
tute of Physics The basis for the work is the card 
file of X-ray diffraction data published by the ASTM 
\j.". The information in the file was assembled 
from all available sources by a subcommittee, under 
the chairmanship of W. P. Davey, representing the 
organizations noted above Data for inorganic and 
chemicals and minerals, given in the 
cards, include interplanar spacings, relative inten- 
sities of the diffracted beams, and, where available, 
planar indices, unit-cell dimensions, specific gravity, 
and optical properties 


as 


organic 


The objectiv es of the project are to record standard 
X-ray diffraction patterns for compounds of high 
purity comparable to the best patterns produced by 
other types of apparatus for the ASTM file. These 
patterns are used as a basis for comparing and 
evaluating duplicate file patterns, and as new infor- 
mation to be added to the file. For this purpose 
X-ray diffraction equipment utilizing a Geiger- 
counter spectrometer with a 160° effective are is 
used 

Considerable progress has been made since 1920 
in the production of X-ray diffraction patterns. In 
most older patterns the interplanar spacings are 
fewer in number and the precision of the measure- 
ments less than that obtainable with present equip- 
In some instances the first low-angle line or 
and in other cases they do not fit 


ment 


lines are missing 


constants 


included to compare the patterns with 
follow for X-ray diffraction data based on four compounds not 
\merican Society for Testing Materials file and on four represented there 


and magnesium tungstate, MgWOQ,, 
- barium nitrate, 


HgCly, 


X-ray diffraction pattern file of the American Society for Testing 
at 
t Committee on Chemical Analysis by X-ray Diffraction Methods of 
for Testing Materials, the American Crystallographic Association, and 
The equipmen 
with a 


the National Bureau of Standards with the 


t used in obtaining the data presented here 
180° are calibrated with the calculated 


are known with high precision. Separate 


to obtain data for interplanar spacing and intensity measurements so 
al attention can be paid to mounting the sample for the best results in each case. 


those recorded in the literature. 


The 
forsterite), 
The 


mat- 


oxide, magnesium silicate 
are recommended as additions to the file 
Ba NO : lead fluochloride 
recommended to replace the present 


SseQ,; 


are 


the unit-cell dimensions determined from the mor 
accurate high-angle lines or from data obtained by 
other techniques. In addition, many of the patterns 
appearing on the ASTM diffraction pattern ‘4k 
cards are combinations of several sets of published 
and previously unpublished data, and the resultias 
patterns are consequently difficult to evaluate. Thi 


present program to improve and clarify the ASTM 


file should add to its usefulness 


II. Equipment 


Thespectrometerassembly is based on X-ray equip- 
ment of the type designed for film cameras, modified 
for Geiger-counter recording. A Machlett X-ray 
tube is mounted horizontally, so that the focal spot 
has an effective focusing cross section of approxi- 
mately 1.0 by 0.2 mm. For most work a coppe! 
target and nickel filter used. Two rectifying 
valves contribute to linear Geiger-counter response 
An arc of 180° was made by joining two 90° Geiger 
counter spectrometer ares [2]. The specimen holde1 
is designed to allow oscillation or rotation of the speci- 
men in the plane of its surface to increase the numbe 
of contributing crystal particles. The diffracted 
X-ray beam is received in the Geiger counter as it Is 
moved along the arc by a synchronous motor. The 
Geiger-counter the diffractior 
maxima as peaks on a potentiometer chart 

Auxiliary equipment includes a small blower di- 
rected on the specimen holder, and, in summer, an 
air-conditioning unit, which aid in maintaining the 
temperature of the specimen close to 26° C. A dry- 
box permits specimens of low stability to be mounted 
in dry air, or in an atmosphere of nitrogen or helium 
Many unstable samples can be mounted satisfactorily 
as flat specimens if the powder is mixed with enough 
petrolatum to coat each particle and the mixtur 
Samples that are affected 


ure 


response represents 


smoothed into the holder 
by petrolatum may be sealed into a cell with a thin 
cellophane window. 
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The Geiger-counter diffraction spectrometer is 
capable of greater precision in determining inter- 
planar spacings if patterns of a few well-known 
materials are used to construct calibration curves for 
the full extent of the Geiger-counter are travel. The 
unit-cell measurements given by Jette and Foote [3] 
are particularly useful as a basis for the calculation of 
spacings suitable for calibration. Aluminum, silver, 
vold, silicon, and tungsten dimensions are given at 
25°C. Sodium chloride [4] is also a good calibration 
standard, but the coefficient of expansion is relatively 
large, so that temperature is a critical factor. Dia- 
mond [5] dust is useful because its coefficient of 
expansion is low enough to be disregarded. 


III. Procedures 


Interplanar spacings and intensity measurements 
are obtained from separately prepared charts, be- 
cause the type of sample mount that gives the most 
precise interplanar spacing data is the one least 
suitable for satisfactory intensity data. 

For the spacing measurements the powder sample 
is packed flat on a glass slide, where it is confined 
within a thin plastic ring cemented to the glass. 
Special care is given to obtaining a smooth flat sur- 
face, which facilitates close alinement with the 
spectrometer center. An error of 0.05 mm in the 
position of the specimen surface relative to the spec- 
trometer arc center causes a change in the apparent 
diffraction angle 26 of 0.032° at 0° of are, which 
decreases to zero at 180° of are. Errors due to such 
surface alinement difficulty may be compensated by 
corrections based on a few lines from a calibration 
substance used as an internal standard. Errors due 
to the penetration of the X-rays into samples con- 
taining elements of low atomic number may be 
corrected in the same manner. A small amount of 
the standard suffices, as few lines are needed for 
comparison. ‘Tungsten contributes few lines to com- 
plicate the pattern, and its back reflection lines are 
unusually strong; less than 5 percent of the powder 
added to the sample is usually sufficient. In running 
the spacing charts a motor speed of \s rpm is used to 
move the Geiger-counter arm min, covering the 
160° of are in about 5 hr. The recorder charts are 
run so that 1° of are is equivalent to 1 in. on the 
chart paper. The interplanar spacings are obtained 
from the reflection angles directly by reference to 
spacing tables [6]. In general, the angles are meas- 
ured with an accuracy well within +0.04°. With 
copper radiation the largest measurable interplanar 
spacing is 44 A. The definition of the diffraction 
peaks is such that those due to the Ka; and Ka, 
wavelengths are usually resolved at angles greater 
than 70° with copper radiation. 

For intensity measurements an additional chart is 
run. Random orientation of the crystal particles is 
achieved by loosely drifting the powdered sample, 
the particles of which are usually less than 0.025-» 
diameter, into a cell fronted bv a window of 0.005- 
mm cellulose acetate. The cell is wide enough to 
intercept the full width of the X-ray beam at the 
lowest angles of the pattern and thick enough to 


accommodate the greatest penetration of X-rays at 
high angles. For most materials | mm is thick 
enough. The specimen preparation outlined avoids 
the packing associated with many types of flat 
specimen mountings, which usually modifies the peak 
intensities so that they do not compare well with 
those made from a randomly oriented cylindrical 
mounting. In preparing the charts from which the 
intensity measurements of the patterns given here 
were made, the Geiger-counter motor was run at | 
rpm. Recent work indicates that a motor speed of 
4% rpm reproduces peak heights somewhat better. 
As the amplitude record of the instrument has been 
found linear for the counting rates used, the heights 
of the peaks are regarded as close approximations to 
the actual intensities of the diffracted X-ray beams. 
Measurements at diffraction angles less than about 
70° include the intensities of both the Ka; and Ka, 
peaks. The peaks are measured vertically from a 
base line drawn at background height and are valued 
relative to the line of strongest intensity. Repeated 
intensity measurements for the stronger lines are 
found to vary less than 5 percent with separate 
samples of the same material. The resulting data 
are comparable to those from both the Debye- 
Scherrer type of camera and the usual Geiger- 
counter spectrometer procedures, if absorption errors 
introduced by the round specimen of the first are 
kept small and if random orientation is attained in 
the flat specimen of the latter. 

Table 1 includes a comparison of intensity meas- 
urements from film camera and Geiger-counter 
equipment. The data for MgO in column (a) are 
based on careful film measurements made by Frevel 
[7], using molybdenum radiation. Column (c 
comprises Frevel’s calculated intensity values. In 
columns (b) and (d) both sets were converted to 
measurements based on copper radiation by means 
of the conversion chart in the Index to the ASTM 
diffraction pattern file in order to compare them 
with the NBS measurements of columns (e) and (f 
A portion of a large single crystal of MgO was ground 
to pass a 325-mesh sieve. This powder, pressed 
lightly into the ring on the glass slide, formed a 
highly oriented specimen. The effect of the pressed 
sample on the intensity pattern, as shown in column 
(e), is typical of that of powder mountings used in 
many laboratories for routine identification with 
Geiger-counter diffraction equipment. The packing 
of the erystal particles is influenced by the prominent 
cleavage of magnesium oxide parallel to 100, and 
causes a reflection of considerably higher relative 
intensity for this plane than occurs with a ran- 
domly oriented sample. When the same sample 
powder was remounted by pouring the powder into 
the cell used in this laboratory for intensity meas- 
urements, the data shown in column (f) were ob- 
tained. The agreement of this pattern with Frevel’s 
intensity measurements is well within the limit of the 
experimental error. 

In addition to the tabulated interplanar spacings 
and intensity measurements, the standard patterns 
are accompanied by other pertinent data. The 
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Miller indices of the diffraction planes are calculated 
from published unit-cell dimensions. From the 
indexed lines the unit-cell measurements of the 
specimen used for the pattern are precisely de- 
termined and the density calculated. The phase 
purity of the specimens is checked by the petro- 
graphic microscope, and the indices of refraction of 
the material are determined when possible. Effort 
is made to obtain specimens of high purity upon 
which to base standard patterns, although, in general, 
small amounts of impurities do not significantly 
affect the pattern. Indeed, temperature variation 
may be just as important, for the presence of 0.01 
percent of a foreign material in solid solution may 
produce, at most, a barely perceptible change in the 
size of the unit cell, whereas the change of a few degrees 
produces in some materials a substantial difference 


TABLE 1 Comparison of intensity measurements for 


magnesium oxide 


same 
powder as 
(e poured 
into cell, 

NB 


Film camera, Calculated 
revel values, Fre 


Radi 
ition 


The patterns prepared are checked with all 
available patterns previously made by other workers. 
References to the published literature are obtained 
by systematic search through Wyckoff’s compilations 
[8], the Strukturbericht [9], Chemical Abstracts [10], 
and the ASTM cards [1]. The type of equipment 
used, radiation, unit-cell dimensions derived, space 
groups, density, and other data are noted for com- 
parison. Patterns usually are published in a variety 
of forms, requiring recalculation for comparison. 
The values of @, sin @, or sin? @ may be listed, rather 
than interplanar spacings. Intensity measurements 
may be given as absolute values, as relative values 
on the basis of 1, 10, or 100 for the strongest line, or 
as visual estimates. In the preparation of the 
following comparison tables, these are reduced to 
interplanar spacings, d, and relative intensities, J, to 
the base 100, except that visual estimates not 
expressed numerically are retained. 

Interplanar spacings and unit-cell dimensions for 
published patterns may appear in kX units or in 
angstrom units variously defined, depending on the 
accepted value of the radiation wavelength used at 
the time the pattern was made, which may or may 
not be given. In general, the difference is small; 
however for accurate comparison, the interplanar 
spacings and unit-cell dimensions are recalculated 
to the angstrom unit adopted in 1946 [11], which is 


used for the standard patterns. In tables 2 to 9 th 
recalculation of published patterns was usually mad 
by dividing the spacings by the wavelength of th 
radiation originally used and multiplying by th 
wavelength now in use. In many cases the wave 
length was not specified; then for data publishe 
between 1930 and 1947 it has been assumed that k.\ 
units were employed, and conversion to angstrom 
was made by use of the factor 1.00202; conversior 
of earlier data without specified radiation wave 
lengths was not possible, due to the variety o! 
wavelengths in use. The interplanar spacings of eac! 
pattern previously published are not rounded off 
although in many instances it is plain from the 
variation in the unit-cell calculations, particularly 
with cubic patterns, that one or more of the figures is 
not significant. For the sake of simplicity and 
uniformity of treatment, the practice was adopted 
of averaging the unit-cell dimensions from the last 
five lines, those of smallest interplanar spacing, of any 
This has generally given a serviceable figure 
For the tetragonal patterns @ was determined from 
as many A000 and hAkO lines as possible, and the 
average used to calculate c from the remaining lines 
These were averaged after eliminating the more 
divergent values. 


set. 


IV. Patterns 


The following X-ray diffraction patterns for 
exceptionally pure chemicals available to the Na- 
tional Bureau of Standards were measured. Inter- 
planar spacings and intensities are listed in the 
accompanying tables, in which they are compared 
with the results of other investigators. Cubic: Zine 
borate, ZnB,O,; strontium nitrate, Sr(NO;).; barium 
nitrate, Ba(NQO,),. Tetragonal: Lead fluochloride 
(matlockite), PbFCI; selenium dioxide, SeQ,. Ortho- 
rhombic: Mercuric chloride, HgCl,; magnesium 
silicate (forsterite), Mg,SiO,. Monoclinic: Magne- 
sium tungstate, MgWQO,. There are no patterns in 
the ASTM file for ZnB,O,, SeO,, Mg.SiO,, or MgWO,,. 
It is recommended that the following patterns in the 
ASTM powder diffraction pattern file be replaced by 
National Bureau of Standards standard patterns: 


Old file 
number 


New file 


Compound numbes 


| Hanawalt Rinn, and 
| Frevel 


| pritish Museum 


IN ieuwenkam p and 
| Bijvoet 


lH inawalt, Rina, and 
uate 


|r ikken and Harang 
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1. Zinc Borate, ZnB,O, 


No published pattern for zine borate was found. 
The sample used for the pattern of table 2 was one 
of the phosphor preparations of the Radio C orpora- 
tion of America [12], sample XII-17, of high purity. 
The unit-cell edge derived from an average of the 
pr obtained from the last five lines is 7.4726 A, 

26° C. The lattice derived from the powder pat- 
tern is body-centered cubic, “— two molecules in the 
unit cell. The density is 1.2017, from the X-ray 
data. The index of refraction for the sample was 
determined as 1.739. 


TABLE 2. Zinc borate ZnB, , 





Copper radiation, A = 1.5405. Temperature, 26° C 
Akl d I a hkl d I a 
110 5.29 f 7.48 444 1. O78 1 
20 4.74 ; 7.48 710 1. 0568 1 
211 +. O48 100 7. 406 40 1. 0365 1 
310 2. 364 23 7.476 721 1. 0169 3 
222 2. 158 1 7.476 42 0. 9001 1 
s21 1. 997 2 7.472 730 9812 2 
400 1. 860 13 7. 476 732 490 l 
411 1. 761 $8 7.471 S11 vIdS ; 
420 1. 672 2 7. 477 S20) O82 1 
$32 1. 504 3 7.477 653 SUS2 l 
22 1. 526 2 7.476 822 SSO7 1 
10 1. 466 7. 475 R31 SORT ; 
521 1.364 & 7. 471 “62 S573 1 
440 1. 321 4 7.473 752 M462 1 
530 1. 282 ; 7.475 v0 $252 1 
600 1. 246 l 7. 476 S42 S153 1 
ll 1. 213 2 7.477 921 ROSS 1 
420 1. 1817 l 7.474 4 7006 1 
‘41 1.1551 7.473 930 7877 1 
61 1. 1025 l 7. 478 
Average last 5 lines 7. 4726 


2. Strontium Nitrate, Sr(NO;), 


A pattern for strontium nitrate is compared in 
table 3 with two previously published patterns. 
The first, by Vegard [13] in 1922, was well indexed 
and misses few lines, although it is of less precision 
than the later patterns. The copper radiation used 
was assigned a wavelength of 1.54 A, close to the 
currently accepted value for the angstrom. The 
data were published as sin @ values and estimated 
intensities. The former were converted to inter- | 
planar spacings for table 3. The pattern by Hana- 
walt, Rinn, and Frevel [14], 1938, is included in the 
diffraction pattern file of the ASTM. It was ob- 
tained with molybdenum radiation with the wave- 
length given as 0.708 kX. The sample for the 
National Bureau of Standards pattern was a spe- 
cially purified material supplied by the Mallinckrodt 
Chemical Works. <A spectrographic analysis shows 
Ba < 0.01 percent and Na < 0.01 percent as the 
only impurities greater than traces. 

From the intensity measurements of the NBS 
pattern, the three strongest lines are the 111, 311, 
and 222, consistent with the index lines of the ASTM 
card for the Hanawalt, Rinn, and Frevel pattern. 
The unit-cell edges derived for the three patterns, 
7.81, 7.799, and 7.779 A, show considerable varia- 
tion. Other published unit-cell determinations in- 


TaBLe 3. Strontium Nitrate: Sr(NOs)> 





1922 1938 1951, 
Vegard Hanawalt, Rinn Swanson and Tatge 
Cn. LB4 A ind Frevel Cu, 1.5495 A 
4 Mo, 0.700 A 26° C 
Aki d I a d I a d I a 
111 154 m 7.86) 4.41 100 7.81 4.48 100 
20 «3.92 m 7.84 5 92 $3 7. 84 3. 88 13 
210 3. 53 “ 7.89 49 a3 7. 80 4.474 2! 
211 3. 22 “ 89 | 3.19 i 7.81 3.175 14 
220 | 2.7 “ 7. 87 2.749 19 
1 2.37 s 7.86 | 2.36 100) 7.3 2. 346 72 7.781 
222 | 2.27 s 7.85 | 2.24 100 7. 76 2. 246 4 | 7. 780 
mw 1% m 7.84 | 1.94 17 7. 76 1. 045 12 | 7. 780 
411 1, 836 2 7. 789 
31 180 m 7.86 | 1.78 17 7. 76 1. 785 16 | 7. 781 
420 | 1.75 m 7.84 | 1.75 17 7.83 1. 740 12 | 7. 782 
422 1.00 m 7.81 1. 58 17 7.74 1, 589 10) 7. 784 
33 1. 51 m 7.54 1.50 17 7.79 1, 498 12 | 7. 784 
21 1. 420 2 7. 778 
40) «(1.39 m 7.84 | 1.379 17 7. 801 | 1.376 1! | 7. 784 
531 1. 32 s 7.83 | 1.318 17 7.797 | 1.315 10 | 7. 780 
Hoo 1.30 w 7.81 1. 206 4 7.776 
620 | 1.24 “ 7. 82 1, 231 2 7. 786 
m3 1.20 m 7.34 1. 1867 4 | 7. 782 
622 1.176 m 7. 80 1. 1736 1 | 7. 785 
444 1. 128 w“ 7. 82 1. 1235 2 7.784 
711 1.0% m 7.81 1. ORGS $ | 7.779 
642 «1.045 m 7. 82 1. 0896 + 7.780 
731 1. 017 s 7.81 1. 0128 7. 780 
732 0. 9878 2) 7.778 | 
820 0.951 s 7.54 9435 4 7.780 
422 go? s 7. 86 vlOs 4 7.779 
751 wo2 s 7. 81 SUNS 5 | 7.780 
S40 873 m 7. 80 
911 857 8 7.81 
S42 SS! s 7.80 
tw AS “ 7.81 
| Average unit-cell 
| edge from last 
5 lines 7. 81 87.799 7.779 


* Average of last 2 lines only 


| clude 7.80 A, by Jaeger and Van Melle [15] in 1928, 
given without stating the radiation used; and 7.811 
A by Ringdal [16] in 1932, measured from a chant 
showing unit-cell variations in solid solutions, iron 
radiation, and with no wavelength given. Vegard 
and Roer [17] in 1942, studying the coefficient of 
expansion of nitrates, gave the unit-cell dimension 
of strontium nitrate at 10° C as 7.7629 A, the coeffi- 
cient of expansion between 10° and 70° as 2.58 * 107°. 
They state that the unit-cell edge of sodium chloride 
used for calibration was 5.628 A. Thus, assuming 
that the values of Ringdal and of Vegard and Roer 
are in kX units, the following table shows compara- 
tive values, after conversion to angstroms: 


Unit-cell dimensions, angstroms 


81 


1922 Vegard 

1W32 Ringdal | 7.827 
1938 Hanawalt, Rinn, and Frevel 7. 79 
1042 Vegard and Roer Cc 7. 7818 
1951 Swanson and Tatee (26° ¢ | 7.779 


The lattice of strontium nitrate is simple cubic, 
four molecules to the unit cell. The space group 
according to Jaeger and Van Melle is 7% (Pa3 
Vegard and Bilberg |16] confirm this, but indicate 
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the possibility of 7''(P2,3). The patterns of table 
3 show AkO only if / is even, adding confirmation of 
the 7{ (Pa3) group. The density from the NBS 
unit-cell edge value is 2.974. The index of refrac- 
tion is 1.587 


3. Barium Nitrate, Ba(NO;) 


that for 


pattern for barium nitrate closely parallels 
strontium nitrate. Vegard [13] and Hana- 
walt, Rinn, and Frevel [14] published patterns, of 
which the latter is included in the ASTM file. In 
Vegard’s paper the copper radiation was assigned, 
as for his strontium nitrate pattern, a wavelength 
of 1.54 A, and the data were published as sin @ 
values. For comparison with the data in table 4, 
they were converted to interplanar spacings. Hana- 
walt, Rinn, and Frevel used molybdenum radiation 
with a wavelength of 0.708 &.Y, converted in table 4 
to 0.709 A. The National Bureau of Standards 


sample was specially purified material supplied by 


raspie 4 





the Mallinckrodt Chemical Works. Spectrograp| 
analysis indicates the following  impuriti 
Al<0.01 percent, Na<0.01 percent, and Sr<0. 
percent. Copper radiation with a wavelength 
1.54050 A was used. 

The three strongest lines, used as index lines for t! 
ASTM cards, are the same for the NBS and t| 
Hanawalt, Rinn, and Frevel patterns: 311, 111, ar 
222 lines. The unit-cell edges derived for the thi 
patterns compared in table 4 are 8.14, 8.130, a 
8.119 A, respectively. Other published unit-c 
determinations include 8.10 A by Jaeger and \; 
Melle [15] in 1928, radiation not stated; 8.111 | 
Ringdal [16] in 1932, measured from a chart showir 
unit-cell variations in solid solutions, iron radiatio 
no wavelength given. Vegard and Roer [17] 
1942, studying the coefficient of expansion ( 
nitrates, gave the unit-cell dimension of bariu 
nitrate at 10° C as 8.0985 A, the coefficient of ex 
pansion between 10° and 70° as 1.7510~°. The 
state that the unit-cell edge of sodium chloride use: 
for calibration was 5.628 A. Assuming that th 
values of Ringdal and of Vegard and Roer are 
kX units, the following table compares the various 
lattice constants after converting them to angstrom 
units: 


‘nit-cell dimensions, angstroms 


Ve rd 

Ringdal 

Hanawalt, Rinn, and Fr 
Vegard and Roer (26 


Swanson and Tatge (26° ¢ 


The lattice of barium nitrate is, like that of 
strontium nitrate, simple cubic, four molecules to 
the unit cell. As found for strontium nitrate, the 
patterns of table 4 confirm the determination of the 
space group 7¥ (Pa3) by Jaeger and Van Melle 
[15] and by Vegard and Bilberg [18]. The density 
calculated from the NBS unit-cell edge dimension 
is 3.244. The index of refraction is 1.570. 


4. Lead Fluochloride, PbFCl 


The new ASTM X-ray diffraction file includes two 
patterns for lead fluochloride, one, also in the old 
file, of natural matlockite from the type locality 
of Matlock, Derbyshire, Eng., furnished by the 
British Museum (Natural History), London, the 
other, in the second supplement, from artificial 
material, published in 1932 by Nieuwenkamp and 
Bijvoet [19]. In 1933 Nieuwenkamp [20] compared 
the patterns of matlockite, whose formula was then 
given as Pb.OClh, and artificial PbFCI, showing their 
identity. In table 5 these two patterns are com- 
pared with one prepared at this Bureau from mate- 
rial of high purity obtained from the NBS chemical 
laboratory, which had been prepared as part of a 
project for the precise determination of fluorine 
compounds. 
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Taste 5. Lead fluochloride: PbFCI 


1932 No dat 1951 
N ve a, Swansor | 
British M 7 = 
iB , rats 
( 2 20) ( 1.5405, 26° ¢ 
u I 
“ 4 2 
mn 7 
lol ‘ is ‘ " 
” 2s “ 2.04 1 
Ps) « : 
we 9 24 2 2 409 
112 2 os st) 2. 03 is 
Ww j 2 2 70 J ; 14 
AM) 2 4 1.4 “i 2 24 
1 t 2 1. OS 40 74 
1 18 1&4 s 1 
(Hi SOS 
l 1. TS tit 177 si 780 
4 H54 4 t O i4 
212 2s 1.¢ ‘ l 7 
LD 4 { 1. 558 i) 1. 4 
i 1. 461 tu wil ‘ 
= 1. 452 
w) 4 4 1, 443 tit 1. 448 
221 1 
) 7] 
> ) 
ol |} 4 it 1. 343 ‘ 1. 34 
1 lad hel ~ 1.209 r 
15 »y 
214 ant su nit) 
Me | vi » 
1 ~ 1, 27¢ 21 
223 24 a) 244 
12 7H ” 22? tu 22 2 
T a1 
a) 191 2 
a 1.181 tit IND¢ 
lin a7 ir 1. 1 
I 42 ‘ 1.144 2 
2 1. 1386 
? a) ‘ 1. 1254 
4 1. 0922 
22 at) 4 1, ONE 1 
2 1, 0300 1 
H 02 ti 1. 2¢ ’ 
! O08 tu 1. OO78 
uw 0. ORT 
W753 
O664 9 
a4 O39 4 
412 GHOS 2 
"V7 av? l 
$ 120 1 
42 yINS 2 


Chromium radiation was used for the Nieuwen- 
kamp and Bijvoet pattern, but the wavelength is not 
recorded. The data published do not include inter- 
planar spacings; for table 5 they were calculated 
from the sin® @ values listed, using the wavelength 
2.2896 A for chromium radiation. The interplanar 
spacings of the British Museum pattern, presumably 
in kX units, were multiplied by the conversion fac- 
tor 1.002 for recording in angstroms. While the 
interplanar spacings of the patterns check closely, 
the intensity measurements vary. The NBS and 
British Museum patterns agree that 101 is the 
strongest line, but the Nieuwenkamp and Bijvoet 
pattern shows the last line (312) strongest, with the 
second and third strongest in close proximity. The 





002, the second strongest line of the NBS pattern, is 
unresolved in the others. The third and fourth 
strongest lines of the NBS pattern appear as second 
and third strongest in the British Museum pattern. 
The intensity measurements of the latter show effects 
due to absorption and to variable focusing common 
with the Debye-Scherrer procedure when using 
long wavelength radiation. 

Frevel, Rinn, and Anderson [21] in 1946 repeated 
the 1932 Nieuwenkamp and Bijvoet determination 
of a=4.09, e=7.21 kX, which, converted to ang- 
stroms compares thus with the NBS determi- 
nation 


Unit-ce dimensions, angstroms 


bOI Swanson and Tat 26” ¢ +. 1 


Bannister [22] in 1934 gave the structure as tetrag- 
onal, space group | (Pnmm), two molecules in the 
unit cell, indices of refraction #,=3.145, ¢p—2.006. 
The NBS pattern yields a density of 7.13. The 
sample was too finely powdered for us to determine 
the indices of refraction. 


5. Selenium Dioxide, SeO, 


\{ pattern was prepared from specially purified 
material supplied by the Mallinckrodt Chemical 
Works. Spectrographic analysis showed no impuri- 
ties greater than 0.001 percent. The lines are in- 
dexed in accordance with the structure and unit-cell 
dimensions determined by McCullough [23] in 1937. 
The material is hygroscopic; lines from the mono- 
hydrate are omitted from the pattern given in table 6. 
Although the crystals generally have been described 
as monoclinic [24], the NBS pattern confirms the 
work of McCullough, showing a tetragonal structure 
MeCullough gives the unit cell dimensions a4=8.353, 
e=5.051, the probable space group Di} (P4/mbc), or 
C8, (Cheb), with eight molecules in the unit cell, 
Frevel, Rinn, and Anderson [21] in 1946 used these 
same values. An average of 13 cell-edge values 
based on AkO lines \ ields a=8.35 A for the NBS pat- 
tern; an average of 11 values for c, after eliminating 
the more divergent ones, vields 5.08 A. Converting 
the McCullough values to angstrom units, they com- 
pare with the NBS determinations as follows: 


U'nit-cell dimensions, angstroms 


MeCullous 8.37 061 


51 | Swanson and Tat 26° 8 \s 


The material proved too unstable for a determination 
of the indices of refraction. The density is 4.16 
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TABLE 6 lenium dioxide: SeQO, 


1951 Swanson an 1951 Swanson and Tatge 
Cu, 1.5405, 26 Cu, 1 5, 26° C 


6. Mercuric Chloride, HgCl, 


One of the two patterns for mercuric chloride in 
the X-ray diffraction file of the American Society for 
Testing Materials was published in 1928 by Brikken 
and Harang [25], the other by Hanawalt, Rinn, and 
Frevel [14] in 1938. Unfortunately, the Hanawalt, 


Rinn, and Frevel pattern is omitted from the index. 
Instead, the index lines 4.35, 3.00, 2.70 belonging to 
mercuric chloride are assigned to mercuric chlorate. 
Index of 


(See “Alphabetical X-ray diffraction 
patterns, covering original set of cards’’ for June 
1945, p. 53.) In table 7 the two patterns are com- 
pared with one more recently prepared at this Bureau. 
The Brikken and Harang pattern was prepared with 
iron radiation for which a wavelength of 1.932 is 
given; for the Hanawalt, Rinn, and Frevel pattern 
molybdenum radiation with a wavelength 0.708 
kX was used. The NBS material was obtained 
from the J. T. Baker Chemical Co., No. 101742. 
Spectrographic analysis showed no impurity greater 
than 0.01 percent. 

The data of Brikken and Harang were published 
as a table of Aki indices, sin’ @ values, and intensity 
values visually estimated. For our comparison 
table the sin? @ values were converted to interplanar 
spacings, using the iron radiation wavelength 1.93597 
A. The spacings of Hanawalt, Rinn, and Frevel 
were converted from &X units to angstroms by apply- 
ing the factor 1.002. The 120 line is the strongest 
for all three patterns. The second strongest is the 
200, but this line is not resolved from the 031 line by 
Hanawalt, Rinn, and Frevel, so that their intensity 
measurement is a combination of the two intensities. 
The O11 and 111 are third and fourth strongest in 
the Hanawalt, Rinn, and Frevel pattern, reversed 
for the NBS pattern. The difference in intensity 
is probably too small to be significant; it is not due 
to the radiation used, as the conversion factor for 
molybdenum to copper radiation is close to 1 in this 
range 


TABLE 7. Mercuric chloride: HgCl, 


1928 1938 1951 
Briikken and Hanawalt, Rinn, Swanson and 
Harang ind Frevel 
Fe, 1.936 Mo, 0.709 


In indexing the pattern the unit-cell dimensions 
were taken in the Dana convention, e<a<b 
although the reverse order is given by Wvekoff ‘ 
The unit cell dimensions published in 1934 by 
Braekken and Scholten [26], 5.963, 12.735, 4.325 
converted from kX to angstroms compare thus with 
those derived from the NBS pattern: 


Unit-cell, angstroms 


ind Scholter 


1934 Braekken 
1951 Swanson and Tatge (26° C 


The presence of hkO lines only for k even, and AO 
lines only for A+/ even, agrees with the generally 
accepted orthorhombic space-group determinatio! 
D& (Pmnb) for the crystal orientation used her 
The density of the material calculated from the NBs 
X-ray data is D=5.49 The indices of refractor 
are very high 
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7. Magnesium Silicate, Mg,SiO, 


A sample of magnesium silicate (forsterite) of high 
purity, labeled X-9, was furnished by the Radio 
Corporation of America. It was prepared in connec- 
tion with their phosphor project [12], by heating at 
1,500° C for 1 hr. The indexed pattern is compared 
in table 8 with a pattern published by Clark [27] in 
1946. Clark’s pattern was made with cobalt radia- 
tion, for which he specified a wavelength of 1.786 A 
(kX), that of the NBS with copper, of wavelength 
1.54050 A. The large unit cell of magnesium 
silicate furnishes a large number of possible planar 


reflections. Thus indexing becomes increasingly 
difficult with increasing Bragg angle. As @ in- 


creases, Clark’s interplanar spacings diverge more 
and more widely from calculated values. The last 
17 lines of his pattern were omitted from the table 
because the divergence combined with the multi- 
plicity of possible lines makes indexing purely 
arbitrary. The Geiger-counter intensity measure- 
ments of the NBS pattern show 112 to be the 
strongest line, 131 second, and 222 third, rather than 
the order 222, 131, and 112 estimated by Clark. 
TABLE 8. 


Magnesium silicate: Mg,SiO, 


Magnesium silieath: Mg,SiO,—Con. 


1951, Swanson and 


1946, Clark Co, 
itge 


1946, Clark 


51, Swanson 


1.789 Cu, 1.5405, 26° C 
hkl d I d I 
lf i > 
Os] j 1. 246 . 
270 1. 242 2 
51 1. 208 1 
0) 1, 181 vvw 1. 189 l 
401 1.159 vvw 1. 16 2 
S52 1. 148 vvw 1. 155 1 
272 * 1. 1476 l 
263 \ 1. 1369 > 
421 j os - 
430 1. 1224 l 
412 1, O87 l 
34 \ 1. 0360 
055 { ‘ 
0.10.0 1.0201 
283 \ 0 a7o7 
is] j ei l 
4652 | 
414 1206 1 
2.10.1 J 
1 ge] l 
444 R748 1 
3.10.1 | ‘ 
0.12.0 { S49 I 
a3 } 
2t S331 7 
65 
} 111.3 S281 7) 
037 | 


one Natge 

Co, 1.789 Cu, 1.5402, 26° C 
hkl d I d I 
(20 5.1 vvw 5.11 yi) 
021 eal) ms 1 88 69 
101 . 71 vw 1.73 25 
Mm New ‘ 
120 3 49 vw 487 21 
121 2. 3 vw 000 17 
002 2. 87 vw 
130 2. 75 ms 2. 768 53 
131 2. 0 s 2. 513 73 
112 2. 45 s 2. 458 100 
041 2.34 vvVWw 2 348 9 
210 2.31 vvw 2. 316 Q 
122 2. 26 “ 2 268 50 
1” 2. 24 Ww 2. 2% 83 
211 2.15 vw 2. 161 15 
132 2.03 vvw 2. 034 5 
042 1. 034 vvVw 1. 045 4 
1%) 1. 864 vvWw 1. 878 F 
113 1. 798 vvw 1.811 2 
151 1, 776 vvw 1. 792 
222 1. 737 vs 1. 748 wo 
240 1.729 vvwWw 
241 1. 1 vw 1.470 13 
O61 1. 624 vw 1. 636 12 
133 1. 07 vw 1. 618 15 
152 1. 589 2 
x00) 1. 579 vvw 
043 1 vvw 1. 572 10 
wl 1. 523 vvVWw 1. 331 1 
S11 ] 
213 1.54 vvw 1.514 10 
20 
004 1. 487 “ 1. 497 27 
OH2 1.471 ms 7 0) 
30 1. 424 vvw 1. 448 4 
170 1. 385 vw 1. 306 12 
233 1, 304 4 
322 1. 341 NY 1. 351 17 
134 1. 305 vw 1,316 i) 
2 1. 285 vvw 1.29 2 
24 1. 2h 1 
2i4 | 1. 2h ' 


* The remaining 17 lines could not be indexed and are omitted 


Although several sets of unit-cell dimensions are 
available for the closely related mineral olivine 
(iron-bearing), only one was found for forsterite. 
Rinne [28] in 1923 examined a natural forsterite from 
Vesuvius, for which he found the dimensions a=4.74, 
b=10.26, c=5.99, using molybdenum radiation, 
wavelength 0.7076 A, which compares very closely 
with olivine measurements. Converting from /X to 
angstrom units, his values compare with those derived 
from the NBS pattern thus: 

Unit cell, 


angstroms 


1923 Rinne 4 75 10. 28 & 00 
1951 Swanson and Tatge (26° ¢ 476 10. 20 5 OO 


From the NBS data the cell dimensions were calcu- 
lated only from spacings of planes parallel to one or 
more axes. 

Forsterite is orthorhombic with a 
presumably the same as that specified by Bragg and 
Brown [29] for olivine, V'’, or D& (Pnma The 
density calculated from the cell dimensions of the 
NBS determination, allowing four molecules per unit 
cell, is 3.213. The material was too finely powdered 
for us to determine the indices of refraction 


space group 
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8. Magnesium Tungstate, MgWO, Paste 9. Magnesium tungstate: MgWO,—Continu 


Five patterns of magnesium tungstate are com- 28 194 
a , j M4 : 
pared in table 9 wo are by Broch [30], published ee Broct — ce 
in 1928 and 1930, one by Fonda [31], 1944, and one by i Fe, 1.9 M ” Cu, 1.54 ( ' 


Dunning and \Mlegaw |32 1946 To Broch we owe 
most of the indices He used « oppel radiation for the 
first pattern, iron for the second. In addition to the 
indices his data include diffraction, angles and, for the 
second pattern, estimated intensities The inter- : 
planar spacings listed in table 9 were computed from ‘ Ho 
his reflection angles with recent wavelength values, so 

that they appearin angstroms Fonda used moly bde- 0 

num radiation, other workers copper. The Fonda and 

the Dunning and Megaw interplanar spacings are 


converted from presumed &X units to angstroms z 
For the NBS pattern, material of exceptionally high > 1. 30 
purity was obtained from the Radio Corporation of 
America, marked No. 4, prepared at 1,000° C. The “ na 
NBS pattern includes the first accurately measured 2 A. 215 
intensities ( 202 n 0 
rasce 9. Magnesium tungstate: MgWO, 42 “ 
ia isi 
241 1776 77s 
v28 sane . Me 1951 D4 1 1764 
Br } ne : , “4 D g Swansor 10 ] 1 
Cu | ® ome WW Megaw und Tatg ‘ ) 
h 1. . , Cu 0 Cu, 1.5405, 26° 024 1273 
‘ ] ] ] ~ 
i] 1l4 
/ / I 
LM 1 19 4 
, ri " aN 2 i. 1. (On 1. OOS 
“ 4 ' ‘ s 5! ; 1% ORTS ! ~ 
70 ) { 1. O72 1. O768 
J 2s 4 
, 02 wu The unit-cell constants of the monoclinic mag- 
“4 | 260 2. 4 2. 34 20 nesium tungstate crystals were given by Broch from 
021 2 4 6) m |2% |1 2. 462 i his first pattern as a=4.67, 6=5.66, c=4.92, 8 
sa Vali 1 - (a) - 89°35’, from his second pattern as a=4.68, 6=5.66, 
-—| =o cA, w | se 2. 34 n 2. 344 10 c=4.93, B=89°40’. Converted from /¢X to ang 
8 | 219) : 220 | mw | 2194 m stroms, the later values compare with those derived 
. 217 | m 2.17 2s from the NBS pattern thus: 
“ 0 ve l'nat ce angstroms 
> 4 “ ® § “ 2 O47 
x 2 5 (2t 4 
5 "v ‘ 1. oo ] 
wy 1 ao s02 mw 1. 892 n0 Broch 4 49 ‘ 4.04 
-. ~ ° . ~~ -. 1951 Swanson and Tatge (26° ( 1.69 6S 1. 92 
4 14 _ = 
S| ow 1.7 2 Broch’s 1930 paper gives the space group as 
2 \\ 1. 724 . ; . y 
o| 1 (4,(P2/c) two molecules in the unit cell; the NBS 
= " ee ie pattern conforms to this determination. The den- 
sity calculated from the NBS data is 6.897. The 
eo |e on v2 1. 680 2. material was too finely powdered to work with satis- 
w | 1.652 factorily in determining the indices of refraction 
0 639 | vw 1.639 but it is known that they are very high. 
: , si | mw | 1.578 V. References 
2 | 1.50 05 | 0 1. 502 ( [1] American Society for Testing Materials, X-Ray Diffra 
19 104 100 tion Data Cards, Philadelphia, Pa. (19397); first sup 
il 1 420 ‘ n 49 | plement 1944 second edition, including second sup 
ao; Lan ie w 1.4 2 plement (1950 For a description of this file se 
‘54 - a : ~ 4 Bull. Am. Soc. Testing Materials No. 135, 64 (1945 
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The Half-Life of Carbon 14' 


George G. Manov? and Leon F. Curtiss 


Che half-life 


where C*Q, is used to designate inert carbon dioxide containing some C™O,) by 
# counters that are identical in construction except for the length of the cathodes 


of carbon 14 has been determined by gas counting of C*O CS 


mixtures 
usIng pairs 


The 


suitability of various equations proposed for calculating the true counting rate, Cr, from the 


resolving time 
urements of C, and ¢ 
half-life of C' 
workers, it is estimated 
5.400 — 200 vears 


that the most 


I. Introduction 


The half-life of carbon 14 has been the subject of 
many determinations, the values of which have 
ranged from a high of 7,200 to a low of 4,700 years 
Two general methods have been used: end-window 
counting of solid barium carbonate, and counting of 
gaseous carbon dioxide, in which the radioactive 
carbon is either a minor or a major portion of the 
filling gas of the counter. In general, larger and 
more numerous are necessary to the 
experimental measurements if solid counting is 
emploved 

Libby [1] * and Jones [2] have introduced the radio- 
active carbon dioxide as a minor constituent of the 
gas filling of Geiger-Miiller counters containing the 
customary alcohol-argon mixture. The technique 
employed by Brown and Miller [3] of counting 
Cc” ), Cc ), C'S, abbreviated as C*( ).4 Cs, mix- 
tures represents a notable advancement in the 
analysis of radiocarbon 

The National Bureau of Standards has prepared, 
and is in the process of certifying, standard sources 
of activity of C"™. The process of standardization 
consists of three parts: (1) A collaborative deter- 
mination by mass spectrograph of the isotopic ratio 
of C™ in carbon dioxide prepared at the Bureau from 
active barium carbonate, conversion at the 
Bureau of a weighed amount of barium carbonate 
into active sodium solution containing 
0.1-\f inactive sodium carbonate as carrier and 
approximately 0.0001-/ free sodium hydroxide, (3) 
collaborative analysis of the solution by absolute gas 
counting The standards after certification will be 
available * in the form of glass ampoules containing 
25 ml of the above solution and having an activity of 
the order of 1,000 disintegrations per second per 
milliliter. A full account of the collaborative 
analyses will be presented when all the results of the 
determinations have been received. The purpose of 
the present paper is to offer preliminary values 
obtained at the Bureau for this intercomparison. 


corrections 


carbonate 


Presented at the Chicago meeting of the American Chemical Society Petro- 
leum Division on Sept. 6, 1950 

? Present address, Isotopes Division, Atomic Energy Commission, Oak Ridge, 
Tenn 

* Figures in brackets indicate the literature references at the end of this paper 

* The ammpoules referred to above are now available for distribution by the 
National Bureau of Standards 


{.. and the observed counting rate. ¢ 
and is presented graphically. 
‘is obtained from these measurements 
probable 


, Is evaluated from experimental meas- 

A value of 5,370+ 200 vears for the 
By comparison with the results of other 
value for the half-life of carbon 14 is 


The value in this paper for the half-life of C™ wa 
calculated from the equation 


T,;=0.693 N/(dN/dt), (1 


where N, the number of atoms of C™, was deter 
mined by mass-spectographic analysis, and dN/dt 
(the rate of radioactive decay ) Was measured by vas 
counting of mixtures of C*O,+-CS),. 


II.tDetermination of the Number of Active 
Carbon Atoms 


Radioactive carbon dioxide was obtained by 
heating active barium carbonate * to approximatel) 
1,100° C in a quartz tube. Analysis of the isotopi 
ratio was made by Fred L. Mohler and Vernon H 
Dibeler, of the Bureau’s Mass Spectrometry Section 
The average value obtained from 10 determinations 
was 3.876 +0.002 percent of C', with an estimated 
accuracy of +0.04 percent (absolute), 
known limits to calibration and other procedural 
errors. The probable error was computed from 
the formula 


based on 


0.845335d if \ ‘i 


A quantity of active barium carbonate weighing 
0.1066 +0.0001 ¢ was transferred to a small round- 
bottomed flask containing 1.40 g of inert barium 
carbonate. (The exact quantity used is not critical 
The mixture was treated with an excess of perchloric 
acid in an evacuated system and the carbon dioxide 
absorbed in a tenfold excess of 2-\/ sodium hydrox- 
ide. The large amount of inactive CO, evolved 
aided in sweeping over the active carbon dioxide 
(It should be noted that neither the perchloric acid 
nor the sodium hydroxide were required to be 
carbonate-free.) Quantitative transfer was accom- 
plished by boiling the acidified solution several times 
Examination of the resulting sodium carbonat¢ 
solution by means of a Tyndall beam indicated that 
no barium ion was carried over. 

The sodium carbonate portion was diluted gravi- 
metrically to a volume of 39.00 liters at 25° C with 

* The BaC*O; was made available by the Atomic Energy Commission througt 


the Oak Ridge National Laboratory, and had been especially purified by W. J 
Skraba of the ORNL Chemistry Division 
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1e appropriate quantity of inert Na,CO, dissolved in 
ater to form a solution that contained all the active 
arbon as Na,C'0,;, 0.100-M inert Na.CO, as carrier, 
ind free sodium hydroxide to the extent of approxi- 
nately 0.001 M. The accuracies in weighing the 
itive barium carbonate and in preparing the large 
olume of solution were each of the order of 0.1 
percent. One milliliter of the above solution was 
therefore equivalent to 2.733 +0.005 mg of barium 
arbonate. Analysis by R. A. Paulson, of the 
Bureau’s Chemistry Division, of the original barium 
carbonate for carbon dioxide indicated a purity of 
97.1 percent (with an absolute error of + 0.4 percent), 
and 1 ml of the solution contained® (3.132 + 0.050) 

10'* atoms of C™. Bureau calibrated weights and 
volumetric glassware were used throughout. Experi- 
ments made with A. Schwebel, using an ionization 
chamber and a vibrating-reed electrometer, demon- 
strated that 0.1-\4 sodium carbonate carrier was 
sufficient to prevent the loss of any significant 
quantity of active carbonate by adsorption on the 
walls of the container. 


III. Apparatus and Experimental Technique 


1. Apparatus 


The counter, a sketch of which is shown in figure 
|, consists of an outer Pyrex jacket in which are 
assembled, in order, a l-in. glass spacer, a disk of 
aluminum approximately 1 mm in thickness, a 
copper cathode with a passivated surface, a 1-in. 
glass spacer, a second disk of aluminum, and finally 
a l-in. glass spacer. The outermost spacers are 
sealed to the glass jacket and prevent movement 
within the assembly. The aluminum disks define 








hu Cross section of carbon dioride counter showing 


effective length (L) and diameter (D). 


the active volume of the counter, and a small hole is 
drilled in each, through which passes a center wire 
of 3-mil tungsten. The ends of the center wire are 
sealed to the ends of the outer Pyrex envelope. 
The inner spacers are identical in length (+0.01 
in.); they serve to insulate the copper cathode from 
the aluminum disk and to make the ends of the 
counter identical. It has been shown [4] that the 
efficiency of these counters is 97 +1 percent. The 
number of observed counts per second is, therefore, 
very nearly equal to the disintegration rate. 

A number of experimental counters were built, 
and a total of six was used for the C' measurements. 
Their nominal dimensions (inside diameter and 
length) were 1 in. by 12 in. and 1 in. by 6 in.; % in. 
by 12 in. and % by 6 in.; and \ in. by 12 in. and 3s in. 
by 6 in. Four measurements of the diameter and 
two of the length were made of each cathode, and 
the averages are listed in table 1. 

§ The accuracy in the determination of N, of dN/dt, and of the half-life is given 
in section V 


TABLE Dimensions of copper cathodes used in counters 


Difference Length 


in cathode diam 
volume eter 


Cathode 
volume 


Inside 


Length 
. diamete 


+0. 31 
1. 92 t 
1.928 4 ‘ ; $85 +0. 27 
278 


1. 200 = 0. 10 


The preparation line (fig. 2) consisted of a 25-ml 
round-bottomed Pyrex bulb in which the carbon 
dioxide was generated ; a calibrated mixing bulb 
(146.8+0.1 ml) equipped with a freeze-out tube, a 
small flask with a stopcock through which pure CS, 
vapor could be admitted into the mixing bulb, a 
mercury manometer with a 's-in. bore, and lastly, 
a pair of counters that could be detached from the 
system by means of standard-taper glass joints. 








CALIBRATED BULB 





Glass train for preparation of carbon dioxide 
samples and filling counters 


Ficure 2. 


The total volume of the system (mixing bulb, 
counters, manometer, etc.) was readily determined 
by noting the change in pressure produced by ex- 
panding nitrogen gas from the known volume of the 
mixing bulb (146.8 ml) into the entire system. The 
total volume ranged from 390 to 650 ml, depending 
upon the pairs of counters used and could be de 
termined to +0.2 percent. 

The counting system used included a preamplifier 
and quench circuit of the type used by Mann and 
Parkinson [5] and a conventional scaler. The re 
solving time applied to the counter was measured by 
an oscilloscope and could be varied from 300 to 
1,400 microseconds. 


2. Experimental Technique 


A definite quantity of active sodium carbonate 
solution (0.5 to 1 ml, measured to +0.002 ml) was 
pipetted into the generating flask. The solution was 
frozen by means of liquid air and the flask evacuated. 
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The vacuum was shut off, and a large excess (25 to 
30 ml) of concentrated sulfuric acid was added 
and frozen in a laver above the solidified 
carbonate solution. The liquid air was then removed 
and the flask allowed to warm to room temperature. 
The sulfuric acid, because of its lower melting point, 
was the first to liquefy and, as the carbonate began 
to melt, the evolution of carbon dioxide proceeded 
very smoothly. This technique of adding the sul 
furic acid served not only to generate the carbon 
dioxide, but completely to remove the water from it 
as the reaction progressed This technique can also 
be used with barium carbonate; the barium sulfate 
is soluble in excess of concentrated sulfuric acid.) 
The necessity for the quantitative separation of the 
water from the carbon dioxide, as would have been 
required had the customary procedure of adding the 
acid dropwise to the carbonate solution been used, 
was therefore avoided, and a considerable simplifica- 
tion in the construction of the gas train was rendered 
possible. 

The carbon dioxide was admitted to the calibrated 
mixing bulb, which contained a pressure of carbon 
disulfide (approximately 2 cm) requisite for the oper- 
ation of the counters. The generating bulb was 
warmed and flushed with a quantity of inactive 
carbon dioxide that would be sufficient to yield a 
final pressure between 8 and 15 cm within the count- 
ers. The contents of the mixing bulb was frozen by 
means of liquid air, and any trace of permanent gas 
removed by pumping. The pump was then shut 
off, the liquid air removed, and the gases allowed to 
mix by thermal diffusion, after which the counter 
assembly and the manometer system were filled. 
The counters were then removed from the prepa- 
ration line and placed within a 2-in. lead shield. 


slow ly 


3. Correction for Resolving Time 


It is well known that the use of CO, as a filling gas 
may result in counters having resolving times as high 
as 1,000 microseconds and, further, that CO,.+CS, 


counters are not self-quenching. In this respect, 
such counter fillings are inferior to alcohol-argon, to 
which is added a small amount of carbon dioxide. 
The advantage of using the CO,+CS, mixture, on 
the other hand, is that it eliminates the necessity of 
working with very small volumes of CO, (0.01 to 0.1 
ml) and diminishes any concern regarding the pos- 
sibility of absorbing some of the CO, on the compara- 
tively large surfaces of the sample preparation line 
and on the walls of the counter. 

Various equations have been proposed to correct 
the observed counting rate, C,, to a “true” one, C; 


T- 


) 





Equations 2 [6, 7, 8] and 5 [10] have been deri 
from theoretical considerations. Equations 3 
and 4 represent expansions of eq 5 and 2, respectiv: 

So long as the resolving time, ¢,, is small (appro 
mately 300 microseconds), eq 2 to 5 all lead to ; 
proximately the same percentage corrections even 
fairly high rates of counting, but if the resolving ti 
is of the order of 1,000 microseconds, not all of the 
are equally suitable. ‘The differences in these eq 
tions will be accentuated for high values of ¢, and 
C 

An experimental test of the suitability of eq 
through 5 was made by observing the counting r: 
obtained for a mixture of C*O,+-CS, as the resolvi: 
time applied to the counter was increased systemat 
cally. All other variables, including the voltage a; 
plied to the counter, were kept constant, and t! 
measurements were made at a point approximat 
in the center of the plateau. Since C, 1s measured 
experimentally for each setting of t,, the equation that 
is the most nearly suitable for obtaining a true count- 














ing rate is that which yields calculated values fo: 
(, that show no trend with ¢t,. The results of such a 
test are presented in table 2 and in figure 3. Figur: 
3 shows that eq 2 overcorrects and that eq 3 and 4 
undercorrect. Equation 5 appears to be the most 
suitable, and it was therefore selected for correcting 
all observed counting rates, C,, to a true value, ¢ 
[he accuracy in the correction for resolving time j 
estimated at +1 percent. 
sed I , 
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TABLE 2 Suitahi ily of various equations proposed for 
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observed counting rate, C,, at a resolving time, 
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4. Effect of Resolving Time on Length of Plateau 


The effects of varying the total pressure within the 
vstem, the ratio of carbon disulfide to carbon diox- 
de, the resolving time applied to the counters, ete., 
vere determined. The effect on the plateau of vary- 
ng the applied resolving time is particularly worthy 
if mention. 

In figure 4 are plotted typical plateaus (for a given 
mixture of CS, and C*O,) corresponding to different 
values of the applied resolving time, ¢t,. It is evident 
that longer plateaus are associated with larger values 
of t,. As the correction for the effect of the resolving 
time can be applied by use of eq 5, and a compara- 
tively long plateau is desirable, most of the subse- 
quent experimental work was performed by using 

1,000 microseconds, with two measurements being 
made by using ¢ 1,400 microseconds. 

A typical plateau obtained by using a resolving 
time of 1,000 microseconds applied to a pair of 
counters (1 in. <12 in. and 1 in. 6 in.) filled with 
CS, and active CQ, is shown in figure 5. The lowest 
curve represents the differences in the two counting 
rates. The flatness of the plateaus indicates that the 
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distortion of the electric field of the ends of the count- 
ers is small. A counter 1 in. in diameter and 3 in. in 
length did not give a flat plateau, presumably because 
the distortion of the field at the ends of the counter 
amounted to an appreciable fraction of the length of 
the cathode. Counters whose ratio of length to 
diameter was 6 or more were found to be satisfactory. 


IV. Determination of the Rate of Decay, 
dN/dt 


The rate of decay of the solution contaiming 
NaoC"O; was measured as follows. The number of 
counts per second was determined as a function of the 
applied voltage for each of a pair of counters filled 
with a definite mixture of CS, and active CQ. 
Between 100,000 and 200,000 counts were totalized at 
each voltage setting to obtain a standard deviation 
(,total counts) between 0.3 and 0.2 percent. 
Corrections were then applied to each point for the 
resolving time (eq 5 and for the background to vield 
a net counting rate, (’,. The difference in C,, for the 
long and short counters were averaged. The activity 
of the solution in disintegrations per second per 
milliliter was calculated from 


ACLXV 


= = (6) 
AV X mlX 0.97 


(dis sec) ml 


where AC, is the average difference in net counting 
rate for the pair of counters used, AV is the difference 
in their cathode volumes (table 1), V is the total 
volume of the system, ml is the number of milli- 
liters of the active solution taken, and 0.97 is the 
efficiency [4] of the counting system. 

In table 3 is presented a condensation of the data 
that were obtained. In columns 1, 2, 3, 4, 5, and 6 
are shown, respectively, the nominal dimensions of 
the pairs of counters used, the differences in their 
cathode volumes, the volume of active Na.,CO; 
solution taken, the applied resolving time in micro- 
seconds, the total volume of the system, and the 
average net difference in counting rate. In column 7 
is computed the activity of the carbonate solution in 
disintegrations per second per milliliter. The aver- 
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age value is 1,280 (dis/sec)/ml with a probable error | carbonate. The accuracy ° in N is therefore 
percent 
O.8452d The accuracy in dN/dt involves the error in the ¢ 
of the mean (F Jo 8 (cdis/see ml counting 8/1280 or O06 percent; in the det 
mination of the total volume of the counting syste1 
0.2 percent; in the application of the correct 


iyi | 


solution of NasCO; containing C™ for the resolving time, 1.0 percent; and for t! 
efliciency of the counting system, 1.0 percet 
viffer- \ tay | AVE | aes The accuracy in dN/dt is therefore 2.8 perce 
Count ein| of |Applk ume “Yiter- | ity « The combination of +1.6 percent in N and of +2.8 
tak tem | ee IN| tion percent in dN/dt corresponds to an accuracy 
t.4 percent in the half-life of radiocarbon 
by , oo | es7.0| en7| 1384 VI. Calculation of the Half-Life of C 
by r ue) “2.4 4 ’ l A 
by 2 ~+ ' 13 oo | e37.8| tmns| 1308 The half-life of C' was computed according to 
OD is amd So bee oo ooo S54 we) 1258 | eg 1. N has been previously calculated to | 
2b 2 by . oe r fi) L221 | 3.13210" atoms of C'/ml of solution, and dN /d/ 
) sat sa a " from table 3 is 1,280 dis/sec/ml. The half-life of 
asia , | Cis therefore 5,370 200 vears 
( } itd ) 2 
; VII. Comparison With Other Determinations 
V. Sources of Error The half-life obtained here, 5,370 200 years, Is 
in fair agreement with that obtained by Jones [2 
The accuracy in the determination of the half- 5,089 75 vears) and by Libby [1] (5,580 15 
life of radiocarbon can be estimated as follows years), but differs from that reported by Mann and 
The accuracies in the determination of N involve | coworkers [12] (6,360 200 years and by Yafh 
0.04/3.876) o1 1.0 percent in the isotopic and Grunlund [13] 7,200 500 years A summary 
ratio, 0.1 percent each in the weighing of the of the values for the half-life is presented in table 4 
barium carbonate and in preparation of the solution, mig 
and 0.4 percent in the analysis of the barium PI ernnccanany AF me ne - innate ene 
TABLE 4 Comparison of values obtained for the half-life of C™ 
Method used for determination of 
Year Re Au l 
Isotoy Act 
Ye 
oe ‘) Reid, Dun Weinhouse, Grosse_| Mass spectrometer (1.71% of C" End-window 4, 700 +(235 to 470 
1s N al Inghran Mass spectrometer 23 and 3.35 lo 5, 100 +200 
af C 
14s Yaffe and Grunlund Mass spectrometer (6.2 of C™ lo 7. 200 +50 
Mo 2 Hawkins, Hunter, Mann, Stevens lo C*Oe+-Cs 6, 360 +20 
ee he) Engelkeme mill, Fnel I Ma 1M 2s 6.14 cro leot sO) j 
Libby rc 
1049 2 Jone Mass spectrometer (39.4% of C™ lo 5, 589 +7 
10m) rt Manov a Curt Mass spectrometer Sv of Cu C*Oe+CS: 70 +200 
pape 
The measurements of Libby and of Jones were The value found by Yaffe and Grunlund was 
performed by using single counters in which the | obtained by using an aliquot of the same Na,CO 
major constituent of the filling was alcohol-argon, | solution as used by Hawkings, Hunter, Mann, and 
with the carbon dioxide introduced in a small quan- | Stevens, except that the counting was performed by 
titv. End effects were determined for counters of | using barium carbonate and an end-window countet 
various sizes, and corrections were applied to the 
data, and the assumption was made that alcohol- VIII. Summary 
argon-C"O; counters were 100-percent efficient The isotopic ratio of C™ in a mixture of inert an 


The value for the half-life of radiocarbon would | radioactive BaCO, was determined by a mass-spectro 
be decreased if the barium carbonate were not | metric analysis of the carbon dioxide prepared ther 
entirely pure, or if the gas counters were not 100- | from. A known quantity of this barium carbonat: 
percent efficient. Adsorption of a small amount of | was used to prepare a solution containing activ: 
CO, on the walls of the glass system prior to entering | NasCQOs, the activity of which was determined by using 
the counter would tend to increase the value obtained | internal Geiger-Miller counters. The value for the 
for the half-life half-life of C™ was found to be 5,370+200 vears 
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Heat Capacity of Liquid Mercury Between O° and 450° C 
Calculation of Certain Thermodynamic Properties of the 
Saturated Liquid and Vapor 


Thomas B. Douglas, Anne F. Ball, and Defoe C. Ginnings 


rhe enthalpy 


of liquid mercury was measured from 0 


to 450° C by the drop” method 


lhese and other precise published data were used to calculate a number of thermodynamic 


properties of liquid and gaseous 

38.88", to 500° © 
compared with 
The calculated values 
to agree over a wide 
indepe ndently 
were applied 


temperature range 
derived gas-imperfection 


I. Introduction 


It is of considerable practical and theoretical im- 
portance that the physical properties of mercury be 
known accurately. The element has found an im- 
portant use as the fluid in certain heat engines 
operating at high temperatures. Furthermore, it 
can be highly purified, perhaps more easily than 
almost any other commonly available substance 
Because of their reproducibility, the physical prop- 
erties of mercury have often been used as standards 

The accuracy of the values of many thermo- 
dynamic properties, over a temperature range, often 
depends on how accurately the heat capacity is 
known. As various past observers have shown 
considerable disagreement above room temperature 
in their values for the heat capacity of mercury, the 
measurements reported in this paper were under- 
taken primarily to furnish accurate values of this 
property up to a vapor pressure of 4 atm. This 
investigation is the second in a current series of 
measurements at this Bureau of the heat capacities 
of liquid metals. 


II. Experimental Procedure 
1. Method and Apparatus 


The method and apparatus have been described 
previously [1, 2, 3].' 

In brief, the method consists in heating the sample 
in a furnace to a known temperature and dropping 
it into an ice calorimeter, thereby measuring the 
heat evolved in cooling the sample to 0° C. The 
calibration factor of the calorimeter was determined 
electrically to be 270.46 0.03 absolute joules per 
gram of mercury. The samples were sealed in 
cylindrical containers of stainless steel. The heat 
capacities of the empty containers were accounted 
for by “blank’’ experiments employing them, these 
experiments being carried out at the same tempera- 
tures as with the filled containers. The temperature 


I iW n bracket i te the literat erences at the end of this paper. 


mercury at 
‘The entropy calculated from data on the vapor and liquid was 
that derived from published low-temperature heat-capacity data for the solid 
of vapor pressure, also using data on the vapor and liquid, 


the vapor pressures from the triple point, 


were found 


with certain published experimental values when 
and 


published temperature-scale corrections 


of the sample in the furnace was measured by a 
platinum resistance thermometer that had 
calibrated at this Bureau. 


bee li 


2. Samples 


Two samples of mercury of about 130 ¢ each wer 
sealed in the containers made of stainless steel No 
347, each having about 10-cm* capacity. The cap 
sules had the same mass (17 g) and composition as 
those used in the sodium investigation [3]. The 
samples, purified and sealed, were furnished by the 
Knolls Atomic Power Laboratory, of Schenectady 
N.Y. Commercially pure mercury that was believed 
to have been triply distilled was redistilled four more 
times in vacuum. The samples were sealed in the 
containers under a pressure of helium of about 
atm. The sealing process [4] was completed by a 
pulse of high-frequency current induced locally at 
the top of the container. It was found possible to 
seal the containers in this manner without changing 
their weight by more than a milligram. The “empty 
containers were sealed in the same manner with the 
@me pressure of helium. The containers filled with 
mercury were tested for tightness at 450° C and were 
found to have a leakage rate of about 0.1 microgram 
of mercury per hour at this temperature, an amount 
that is without significant effect on the enthalpy 
measurements. 

The mercury sample actually used for most of th: 
thermal measurements was examined spectrochemi 
cally at this Bureau. Of 34 elements looked for as 
possible impurities, only copper and nickel were de 
tected by this means, and these were found to b« 
present only in traces amounting to less than 0.01 
percent of the mass of the sample. The stock supply 
of mercury from which the sample for the thermal! 
measurements had been taken was analyzed by th: 
Knolls Atomic Power Laboratory. The total non 
volatile impurity found, mostly silver, amounted t 
0.00001 percent. A mass-spectrographic examina 
tion by them for ‘‘volatile’’ impurities indicated th 
possible presence of traces of aluminum, manganes« 
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= iron, cobalt, nickel, zinc, and rhodium; the total TABLE 1. Corrected heat and enthalpy values of individual 
amount of these, however, was shown not to exceed experiments 


0.001 percent. 
C Measured heat . mere amy of 
° nercury, 
h III. Experimental Results Furnas 
t s 


temperature 
With Caleu 








A total of 111 measurements of enthalpy was made, Blank | mercury | ObServed) jateq | Difference | 
from 0° to 450° C (to about 4 atm pressure). Of 
these, 66 were made with the two empty capsules ( abs j ahe abs-j gt | abs-jg-? | abs-j gt 
and 45 with the two capsules containing mercury. { jos. 4 
The average deviation of a single measurement from a. | ae 
the mean at a given temperature was 0.05 percent. 0.00 1.1 1298. 7 6. 942 6. 946 0. 004 
One of the mercury samples, measured at only 250° C ae er 
to ensure that no systematic error in mass was 4 
present), gave a mean value for the enthalpy of pet 
mercury that differed from that obtained when using S140 
the other sample by only 0.01 percent. The detailed 9 : sank 6 
results of individual runs are given in table 1. 100.00 815.9 | 2003.8 | 13.835 | 13.831 | +0.004 
. . . . Ts S16. 2 204.8 
No corrections for impurities were made, as these 813.2 
were undoubtedly so small as to lie well within the art 
accidental error All weights were corrected for a cine 
Dy il buovancy. No corrections for temperature were ne a8 ee ° 4 20. 672 20. 669 +0. 003 
been made, as the thermometer reading was held to within ; oa p14. 8 
0.01 deg of the stated temperature in each case : 
The ice point of the resistance thermometer changed + 9 
so slightly during the course of the measurements | gy + 
as to indicate a negligible error in computing the 200.00 1674.9 0 (427.461 27.470 0.009 
wert temperatures Corrections were applied to account a? 2 4 
No for small differences in masses of capsule and exterior ans - 
cap- oxide coatings. In addition, corrections were made \°1675. 6 
on as at the higher temperatures for the small heats | 6548 
a evolved in condensing some mercury vapor inside aie 2121 og Be ae eee aoe 
the container, in order that the results would refe1 | 2120. 2 6484. 0 
aly to the liquid alone. |=. ..| 0eens 
eved This last-mentioned correction, and that needed 
more to evaluate the enthalpy change that would have 
L the resulted under maintenance of saturation from the 
ut heat measurements made on the system maintained — sotendil bsckaadll Wiainahoes 
by a at constant volume, were conveniently calculated 
ly at by an equation given by Osborne [5], 
le to ; 
Wing i lq PV +mHat+(\ me) L/(v’ vii, (1) aii ~~ oa 17.791 17.790 Sa oni 
itv w2S8. 2 9211.0 
: the where lV}; is the heat evolved in cooling a closed ( 3008.1 — 
with container in which there is a liquid in equilibrium ( Sans. 5 
wert with its vapor from temperature 2 to temperature 1; | 3494.3 © 10548. 9 
vrai qi; is the contribution to [Qi made by the empty | 34034 | 0850.8 | 4540800) BA 877 0.008 
ount container, shields, and suspension wire; P is the — ous | care 
ralpy vapor pressure of the liquid; ' is the internal volume 1e3489.5 | 105475 
of the container; m is the total mass of liquid and 01. 2 
f the vapor; /7,, is the enthalpy, per unit mass, of the ay eee 
lemil- ‘saturated”’ liquid (i. e., at pressure P); ¢ and ¢ | 3068.1 119020 
or as are the specific volumes of saturated vapor and 450.00 is | sen |éenae | aa £6,008 
e de liquid, respectively; and L is the enthalpy of vapori- oe | oo 
Lo bi zation per unit mass. 064 
0.0) The total of the various corrections did not exceed 
ipply 0.03 percent of the enthalpy, except in the case of Values apply to saturation (i. ¢., at the vapor pressur 
mal some of the first runs, where a correction of approxi- es one em ih meant 
y the mately 0.2 percent was necessitated by the use of a Theos a ee ee ~ 4 a> ant pengemenney to Ghe Hanae 
non shield system that was later broken and so had to from the acca cithe set. degre ee ee 
ed t be replaced by one of different heat capacity, BL abe wg Tigers A. a Cae Cer cages Bas oe 
nina 
d th 
nes 
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Within the accidental error the heats of the empty 
containers (“blank’” experiments) were found to 
vary perfectly smoothly with temperature up to 
150° C. By plotting against #(°C) the function 


H,—H 


0.0025 


L.SO6f 
which varied only 1 percent between 50° and 450 
these blank values were thereby smoothed The net 
enthalpies of mercury given im table 1 have been 
calculated by the smoothed values for the 
empty contamer, although the heat values listed for 
the empty contaimers in 
obtained 
The experimental values of 
liquid mercury (at 
enthalpy at 0° C 
equation, 


usInY 
the table are those actually 


the 


pressure 


enthalpy of 
the 
are represented by the following 
fitted by least 


Saturation 


less 


Ww hose constants were 


squares: 
H,—H,=0.139612t—1 


1673110 


a3 
150°C ) 


1L.6874(10 0° to 
where /7 is in absolute joules per gram and ¢ in deg C 
International Temperature The values in 
the fifth column of table 1 were calculated from this 
equation 


Scale 


IV. Calculation of Thermodynamic 
Properties 


1. General Procedure 


The more common thermodynamic properties of 
liquid and gaseous mercury, at the existing vapor 
pressures, were accurately calculated as a function 
of temperature over the range from the triple point 

38.88°, to 500° C. Except for what may be con- 
but unavoidable terms, 
values of of properties 
rest on three sets of precise experimental data, 
namely, (1) the enthalpy measurements of the 
liquid reported in this paper, (2) a previously pub- 
lished series of vapor pressure measurements cover- 
ing a 13-deg temperature range in the vicinity of 
the normal boiling point, and (3) the experimental 
values of the fundamental physical constants that 
made the evaluation of the 
entropy of the ideal The additional exper! 


sidered as minor 
the calculated 


correction 


most these 


possible statistical 


Vapol 


mental data employ ed, W hose aceuracy is of secondary 
effect on the ac uracy of the calculated properties, 


are (1) those giving certain fundamental constants 
of the Hg, molecule (used to arrive at data of state 
for mercury vapor), (2) accurate gas-thermometer 
measurements of laboratory (used to make 
from the International to the thermo- 
dynamic temperature scale), and (3) PVT data for 
the liquid (necessary in an accurate evaluation of 
the heat capacities from the measured enthalpy 

The procedure of calculation followed will now be 
outlined. The heat liquid 


corrections 


capacities of mercury 


were calculated from the enthalpy data by 
usual thermodynamic relations. The changes w 
temperature of the entropy of the saturated liq 
were next calculated from the values of heat capaci 
and were combined with the statistically calcula‘ 
value of the entropy of the vapor (at the vapor pl 
sure to give the absolute entropy of the liqu 
The heat capacity (C,) of the vapor was calculat 
by first assuming it to be an ideal monatomic 
and then making small corrections for imperfectio 
of the By integration of the resulting he 
capacity equation and evaluation of the integrati 
constant by using the value of the heat of vaporiz 
tion at the normal boiling point calculated from t 
vapor-pressure data and the Clapeyron equatio 
values of enthalpy of the vapor (relative to the liqu 
fixed temperature) were obtained. The fi 
energies of the liquid and vapor followed from t! 
calculated entropy and enthalpy values. By equa 
Ing the expressions for these free energies of “sat 
rated” liquid and vapor, a vapor-pressure relatio 
was obtained that is applicable over a much larg 
temperature range than the supporting 
pressure data. 

At the present time, the uncertainty in the co 
rections for gas imperfection and temperature scal 
referred to above, limit slightly the accuracy of thi 
calculated properties Because of unce! 
tainties it was considered advisable merely to ind 
cate these two corrections in the equations derive: 
The values of the tabulated (in tables 
3, 4, and 5) were arrived at, however, by assignin 
to the corresponding eorrection terms, on the basis 
of what are judged to be the best existing data 
specific values that are separately listed (in tables | 
and 7 and in the text preceding them This pro 
cedure facilitates an estimation of the uncertainties 
introduced into the present values and also shoul 
simplify any future desired revision of the values 
should 


gas 


ata 


vapo 


these 


propert ies 


more accurate corrections become availabl 


2. Nomenclature and General Assumptions 


In W hat follows, energy W ill be expressed in calor les 
1 cal=4.1840 abs j); all extensive properties, per 
gram-atom of mercury (atomic weight = 200.61); and 
temperature (7), in degrees absolute (0° C= 273.16 
K), unless otherwise stated Each value calculated 
from the equations applies at the existing vapor pres 
sure, except in those cases where pressure is an ex 
plicit variable in the equations given. In addition 
the subscripts in the heat capacities C,, C,, and ¢ 
signify that the respective heat changes are thos 
occurring under the maintenance of constant pres 
sure, constant volume, a. liquid-vapor equilibrium 
respectively. The equa’ ois for relative enthalp) 
and relative free energy couiain the term //,, the 
enthalpy of the liquid at the triple point, —38.88° ¢ 
Liquid and vapor are distinguished by the subscripts 
/ and g, respectively. In the numerical equations 
i. e., those in which some or all of the numerica 
constants have dimensions), ?, the pressure, is to b 
expressed in millimeters of mercury (760 mm 
standard atm); and V, the volume, in cm* g-atom 
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Log and /n indicate logarithms to the bases 10 and ¢, 
respectively. 

Saturated mercury vapor at pressures not exceed- 
ing 2 atm has been found experimentally to have 
densities that differ from those calculated for an 
ideal gas by less than 2 percent (the uncertainty of 
the measurements). It is therefore permissible to 
assume that practically all of the gas imperfection 
arises from binary collisions only. Statistical me- 
chanics shows that under this circumstance the cor- 
rect form for the equation of state of the vapor is 


PV, rr +i} (4) 


where FB, the so-called ‘“‘second’’ virial coefficient, is 
a function of temperature only. However, for con- 
venience it will be used here in the approximately 
equivalent form 
PV @=RT+PB. (5) 
In the numerical equations of this paper, P is to be 
expressed in cm® g-atom™'. The values of dB/dT 
and @RB/dT* at the temperature at which the equa- 
tion is applied will be designated by RB’ and BY”, 
respectively, and the corresponding values at the 
normal boiling point, 629.74° K, will be designated 
by By, Bs, and Bo. 
6and 7 will be used to designate the values of a 
given temperature on the International and thermo- 
dynamic absolute temperature scales, respectively. 


3. Enthalpy and Heat Capacity of the Liquid 


By extrapolating the experimental enthalpy values 
(represented by eq 3) to thetriple point, —38.88° C 
(234.28° K), and changing to the units as indicated 
above (to cal, g-atom, deg K), the enthalpy of the 
liquid /7,, at any temperature @, relative to the 
enthalpy Hyp at the triple point, becomes 

H Hyp 7.259390— 1.36651 (107-°)@ 
+- 8.0906 (1077) @— 1636.13. 6) 


The heat capacity equations that are derived from 
eq 6 must be multiplied by d@/dT to account for the 
difference between thermodynamic (T) and Interna- 
tional (@) temperature scales. However, in the 
temperature range in which the equations are appli- 
cable, the heat capacity varies by only 3 percent. 
Since dé/dT is, at all these temperatures, very close 
to unity, an error negligible in comparison with the 
experimental errors of the supporting data is intro- 
duced by the more convenient procedure of multiply- 
ing by dé/dT only the average value of (O//,,/08)s, 
taken to be 6.55. This is equivalent to adding as a 
correction term 6.55[(dé/dT) — 1}. 

The “saturation” heat capacity (C, 
from the thermodynamic equation 


) ean be found 


(Ss ) dé . =) dé 


Cra 00 /,dT 00/, dT 


For mercury below 500°C thelast termis small and 
therefore may be approximated with sufficient accu- 
racy by replacing d@/dT by land by substituting for 
Vin the constant value 15.9 cm* g-atom~', equal to 
the volume of the saturated liquid at 400° C [6]. 
A simple vapor pressure equation, sufficiently accu- 
rate for the present purpose, is 


P=6.345(10%)e ~7196-5/, g 


(This was derived from the values of P and dP/dé@ 
given at the normal boiling point by eq 19, and pro- 
vides on differentiation an equivalent of (0P/0dé),). 
With these substitutions, together with the equiva- 
lent of (077/00), obtained by differentiating eq 6, 
eq 7 becomes 


7.25939 — 2.73302(10~*)6+-2.42718(10 baa 


7136.5 >» (9) 


2.294(10*) ail a 
¢  +6.55(d6/d7 1). \ 


A 


The heat capacity at constant pressure may be 
found from the thermodynamic relation 


C,=C,+T(OV/dT),(0P/d7),. (10) 


The last term is small in this case and need not be 
highly accurate. In fact, it will be used with T re- 
placed by @. Deriving (0P/0@), from eq 8 as before, 
and substituting as a constant value for the slowly 
changing (OV,,,/0@), the value at 400°C [6], 3.09(10 


' 
em’ g-atom™'! deg, eq 10 becomes 


44585 
> é 


6 


(11 
The heat capacity at constant volume may in turn 
be calculated from the thermodynamic relation 


eVT 
B 


(12 


where a and are the coefficients of isobaric thermal 
expansion and isothermal compressibility, respec- 
lively. 

In the use of eq 12 to calculate C,, use was made 
of the following equation for the liquid (f being in 
deg. C) 


14.756 + 2.678(107°) ¢+1.36(10 
-9 §(107')8+-9.93 (107) #. 


This equation was derived from the equation of 
Sears [6,7] for the relative volume and from the mean 
experimental value of the density of the liquid at 
0° C. Sears’ equation seems to be based on data 
covering the range 0° to 300° C only. However, it 
fits well the available data for temperatures below 
0° C. Its use in this paper between 300° and 500 
does not introduce an appreciably increased uncer- 
tainty to any of the thermal properties calculated, as 
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the affected terms in the equations are relatively 
small 

For the isothermal compressibility, the values of 
Smith and Keyes |S], extrapolated to zero pressure 
the approximate vapor pressure up to 200° C), were 
used. The values of (, tabulated in this paper 
extend only up to 200° C, because reliable values of 
the compressibility above this temperature are not 
now available 


4. Entropies of Liquid and Vapor 
The entropy of the saturated liquid is 


Te 

a aT + Sow, 14 

n° 2 

where 7) is the thermodynamic temperature of the 
normal boiling point (@—629.74), So.» is the absolute 
entropy of the saturated liquid at that temperature, 
and (, is the heat capacity given by eq 9. Of the 
terms in the second member of eq 9, the sum of all 
except 6.55000 dé/dT is always relatively small, as 
discussed in the preceding section and may thus, 
without appreciable error, be conveniently divided 
by and integrated with to @ instead of 7. 
The remaining term must be divided by and inte- 
grated with respect to 7, 


respect 
fiViInYe 


"6 dé 


6.55000 
620.74 1 


and this will be replaced by 


; )de, 


¢ 
6.55000 In 
H2' 


in order that again the term including the correction 
will be small. Upon integrating, it is found that 


+ 


16.71536 log 6@— 2?.73302(10 “ 


1.21359(110 





$5.548507 S : 


‘Lhe value of the entropy of the liquid at the normal! 
boiling point, Sy), will now be evaluated as a fune- 
tion of the corrections for gas imperfection and 
temperature-scale divergence. The result is eq 24. 

The entropy of the ideal vapor (S°), with ideality 
assumed to imply monatomicity as well, is given by 
the well-known statistical mechanical 
Sackur-Tetrode equation The ground state of the 
mercury atom is singlet. As for all values of entropy 


SO-Ci lled 


328 


in this paper, nuclear-spin and isotopic contributi 
are omitted.) The equation is as follows 


Rin (=7"") 


J (ke T) 
P 


where FP? is the vas constant, m is the mass of t 
atom, A is Planck’s constant, & is Boltzman: 
constant, and P is the pressure. The entropy of 1 
real gas at the same pressure can be shown therm 
dynamically to be less by PAB’ (see eq 5). Subtra 
ing this and adding the temperature-scale correcti: 
lead to 
») . id 4 
S, Rin (=> ) a +2 r+? Rint PB 


1.98719 cal g-atom 
6.624107) erg 
and 
eal, eq 


Making the substitutions R 
m 3.3308(10 2 ¢-molecule”', h 
sec, k 1.38047(10 erg molecule 
cm'*-mm (of mercury) —=3.186484(10 
becomes, for any pressure P, 


deg 
1 


S 


» 
1S 
PR’ 


3.186010 


T 
11.44 log 
0 


Beattie, Blaisdell, and Kaminsky [9] have accu 
rately measured the vapor pressure of mercury from 
349° to 362°C. They have given the normal boiling 
point (P=760) as 356.58° C, and in this temperaturs 
range the vapor pressure equation (t=deg C Inter 
national), 

760) 


t= 356.580 +7.30951(10-°\(P 


3.986610 P—760) 


L.3.191(10~*)\(P—760) 
from which their individual determinations vary by 
only 0.001 deg on the average. Equation 19, which 
will be adopted in what follows, gives a value of 
13.6808 (mm deg for dP/dé at the normal boiling 
point. This value will be substituted into th 
Clapeyron equation to obtain a value for the entropy 
of vaporization at this temperature. Substituting 
for the vapor volume its equivalent from eq 5, th 
Clapeyron equation may be written 
sa 
So—S = ( 


oP , 
ar). +'B- Va) (Sp 


oP 
(Sr). 


If 7 be replaced by 6, only the first term of the secon 
member is large enough at the boiling point to requir 
the addition of a correction term, which is obviously 


R oP 4 dé 
P\ 0¢ ) ( dT 











ecu 


rom 


ling 
ture 


ite! 
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Equation 20 then becomes 


. Re , oP 
8 S ( pte Vu (55 ) 
R oP ple 
+5 (So ) (Top 0): (21) 
Letting @= 629.74, at which temperature V 15.748 


6], applying the subscript ‘‘0” to designate this 
temperature, and substituting the above value 
13.6808 mm deg™'! for (OP/08),, eq 21 gives 


, 22.519887 +4.359(10-B, } 


0.03577(T—8), ( 7) 


.., (d0 
22.53 (oa l)- 


At @=629.74, eq 18 gives 


r (22) 





T 
Sow =45.508660+ 11.44 log ( ) 0.02422 Bi. 
; G /o 
23) 
Substituting for S),,.. from eq 23 into 22 gives 
Sou) = 22.988773 —0.02422 Bi—4.359(10~) B, ) 


0.03577(T Bo ( =) 


11.44 log ( 4) 





dT), > 
oe da 
22.53 (a—1 )- ; 
24 
And finally, the value of Sy. may be substituted 


into eq 15 to give an equation for the absolute entropy 
of the “‘saturated’’ liquid at any temperature in the 
range covered 


S 16.71536 log @—2.73302(10~)@ 
1.2135910-%6 1.511( 9 1 )e 
> T) > » + 
sie ; if] , 


22.559734—0.02422B)—4.359(10-9B, 


6.55 [ . ; - ) de + 11.44 log ( 2 ) 


dé _, (dé 
Ao ( a ) 22.53 (op ] )- 


O0.03577(7 





4 


a 


Ul 


Enthalpy of the Vapor; Free Energy; General 
Vapor Pressure Equation 

The relative Gibbs free energy of the liquid may 

next be found from the entropy and the relative 

enthalpy, by using the thermodynamic definition 


F=H—TS. (26) 


For this purpose eq 26 may be written in the form 


Fy —Hrpa AS H A rpa))+(T—OS « 
From the thermodynamic relation 
(0S 
p=T | or). (28) 


and eq 18, there is obtained for the heat capacity 
of the gas at constant pressure 

Cy @ = 4.96797 —3.186(10-°) PeB” 29) 
Integration of eq 29 with respect to 7’ at constant 
pressure yields the enthalpy of the gas at pressure P, 


H ¢, =4.967976+-3.186(10~°)[B— eB |P 
4.97(T—0)+ A. 30) 


The integration constant A may be expressed in 
terms of Hyp, the constant used in the case of the 
liquid, by utilizing the value obtained from eq 22 
for the enthalpy of vaporization at the boiling point. 
Multiplying eq 22 by the temperature 7), (629.74 

T—@),) and neglecting the very small terms as 
A ool there results: 


Ao, —H. 14181.674+ 0.2745B, > 
—— dé 





1.419(108 ( - |) 
f 4 


Evaluating eq 6 at 6=629.74 gives 
Hy, H-,; 2595.531. 32) 


By evaluating eq 30 at 6=629.74 and substituting 
the resulting equivalent of H),,, into the sum of 
eq 31 and eq 32, an equivalent of the constant A in 
terms of Hypa, is provided that when substitut ed 
for A into eq 50 vives 





H «¢, —H yp, = 13648.676 + 4.967970 7 
20226 ~*-5/*| B— 0B’) 4-0.2503B 
15.25Bj+22.53(T— 6), ( +1) | 

dé 
1.419(10 (or 1) 
4.97((7 — 0)—(T—9®),] J 


In this equation, the vapor pressure has been evalu- 


ated (with sufficient accuracy, for the correction 
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term, from eq 8) and has been substituted for P 
Subtraction of eq 6 from eq 33 will give the heat of 
vaporization at any temperature in the range studied 

The relative Gibbs free energy of the gas, as does 
that for the liquid, now follows by the use of eq 26, 
which mav now be written in the form 


Fig — rea aS, H..,— Hr; T—@S 


By equating the free energies of saturated liquid and 
gas from eq 27 and 34; substituting for the enthalpy 
and entropy of vapor and liquid from eq 33, 18, 6, 
and 25; and solving for log P contained in S,, the 
following vapor pressure equation is obtained 


gs 


3340.449 


“ 


log P=11.259839 1.153092 log @ 


L? ORG47(10~> 5 A@—S8 S4090(107>5)8 
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2 500 log : | 


“ 
- 131 . (7 , a8 P35 
+7.817(10~%) ( 


[ 620.74 ( oS 


L dT 


l \ a x 3099.3 
+ L.Sf05 1.1551 log o-+ 
A / 7) 


5.97T3(10~ 48 





2? 652010 et / A § §38(7 —e@ } 
4 


(In the part of the last bracketed term which repre- 
sents the entropy of vaporization, all minor terms 
have been omitted, since they would be multiplied 
by the very small factor (7'—@ In this bracketed 
term, P (contained in the equivalent of S,,)) has been 
replaced by its equivalent in terms of @, by using 
eq 5 

If the specific corrections for gas imperfection and 
temperature-scale differences described in the next 
two sections be adopted, eq 35 may be approximated 
by the following more readily applicable vapor pres- 
sure equation 


1.153092 log @ 


L929 95697(10~*)@—7.45S88(107°5)80 


1.560510 & + 3.600e — 0? 
Between 100° and 500° C, this equation reproduc 
the vapor pressures of eq 35 to within 0.01 percer 
Below 100°, however, it gives lower values, the di 
crepancy amounting to 0.1 percent at 25° and 0 
percent at 39°, the triple point. 


6. Calculation of the Second Virial Coefficient; 
Corrections for Gas Imperfection 


Smith and Menzies [10] found the density of m« 
cury vapor to be ideal between 360° and 400° C t 
within their probable experimental error, about 
percent. Apparently no one has measured such 
densities with much greater accuracy. 

Approximate values for the degrees of gas im 
perfection at various temperatures have been caleu 
lated for many gases by the use of some equation ol 
Berthelot’s To employ the last 
equation, the critical constants must be known 
In the case of mercury, however, the values that 
have been reported for the critical temperature 
(around 1,500° C) and critical pressure (from 1,00 
to 3,500 atm) show such wide variations that this 
method is hardly applicable 

There is abundant spectroscopic evidence that 
mereury vapor contains appreciable amounts of the 
dimer Hg, which is known to possess a singlet ground 
state. Adopting suitable constants for this molecul: 
as described below, values of the second virial co 
efficient B of mercury vapor, treated as an imperfect 
monatomic gas in which only binary collisions are 
important below 500° C, were calculated from th: 
following statistical mechanical equation [11] 


state, such as 


BT QaNn e~* a Llyredr. 37 


where N is Avogadro’s number, and [(r) is the 
potential energy of a mole of Hg, molecules, all at 
an interatomic separation of r. U(r) was calculated 
from a Morse potential function in which it was 
assumed that the equilibrium separation is 3.2(10 
em, and that the first vibrational constant w, is 36 
em~'. The former figure, derived from data on 
liquid mercury, has been generally accepted as ap 
proximately correct. The latter Meure, considered 
by Kuhn [12] to be the most likely choice among 
alternative multiples for which there is spectroscopi 
evidence, is supported by the theoretical calculation 
of a value of 35 em™~'! by Heller [13], whose values fo: 
a number of other diatomic molecules are in good 
agreement with well-established experimental values 
The third parameter needed in the use of the Mors 
function is the molar dissociation energy of the Hg 
molecule in its ground state. Evidence for various 
values has been described in the literature. London 
[14] made a theoretical calculation of the polarization 
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force and arrived at 2.0 keal mole however, 
Kuhn [12] peints out that this is undpubtedly too 
high because of the neglect of repulsion forces. 
Kuhn measured the total absorption in the discrete 
band system (2,540 A) at 175°C and calculated a 
value for the concentration of Hg,. He claims that 
the extrapolation to 0° K (giving a dissociation 
energy of 2.1 keal) yields too high a value because 
of the neglect of vibrational anharmonicity. If 
this be true, then, as he claims, the values obtained 
by measuring the decrease in absorption with rising 
temperature are too low. By the latter method 
and with the assumption of vibrational harmonicity, 
Koernicke |15] found a dissociation energy of 1.4 
keal, Kuhn and Freudenberg [16] found 1.6 keal, 
and Winans and Heitz [17] found 1.38 +0.07 keal. 
Gaydon [18] recommends the latter value for the 
dissociation energy, whereas Herzberg [19] lists 
1.84 keal. 

Equation 37 was solved graphically by the authors 
for values of B at the temperatures T=430, 530, 630, 
and 730, by using for each temperature values of 1.4, 
1.5, 1.6, and 1.7 keal for the molar dissociation 
energy of Hg,. In these ranges the calculated value 
of B for a given temperature varies almost linearly 
with the value of the dissociation energy used. The 
thermodynamic properties of mercury listed in this 
paper were calculated assuming a value of 1.5 keal 
as a reasonable weighted mean of the above values 
of the dissociation energy. This particular value 
was adopted partly because of the support it re- 
ceives in figure 4, where calculated and experimental 
vapor pressures are compared. The following em- 
pirical equation represents the corresponding values 
of the second virial coefficient B, in em*® g-atom™', in 
this temperature range 

B=56.4—43.82¢", (38) 
which agrees with the directly calculated values to 
within 0.5 percent from @=430 to @=730. On 
differentiation with respect to temperature, eq. 38 
o1lVvVes: 

655 


Bb’ —(56.4—B). (39 
- 


Values calculated from eq 3S and 39 are listed in 
table 2? 


Taste 2 Second al coefficient B) and its te m perature 
derivative (B’) for mercury apor 
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7. Corrections to Basis of Thermodynamic Tempera- 
ture Scale 


The results of recent investigations [20] of the 
differences between the thermodynamic and Inter- 
national temperature scales have been formulated 
by an equation [21] equivalent to the following 


T— 6=0.6381 t.809(107°)@4-1.1096(10 A 
7.481(107°)@ (40 


(The value of (7—@) is so small and changes so slowly 
with temperature that when desirable, 6 may be 
replaced by 7 in the second member of eq 40 with- 
out appreciably changing the corrections calculated 

It is believed that the differences given by eq 40 
are the most probable values in the light of the 
evidence now available, and that their uncertainties 
may be considered to correspond to probable errors 
of about half their respective values. All thermo- 
dynamic properties given in this paper have been 
corrected by this equation to units involving deg Kk 
thermodynamic, though for convenience of usage the 
tabulated values of these properties have been 
calculated for rounded temperatures on the Inter- 
national scale of 1948. 

By use of eq 40 it may be shown that the integral 
appearing in eq 25 and 35 is approximately as follows 


0 l | Le O.635S81 
e ) de 0.01107 log @+ 
620.74 ‘a 


i 


1.1096(10~°86 41) 


-3.7405(10 7°) —0.0264958 


V. Tables 


1. Thermodynamic Properties 

In tables 3, 4, and 5, values calculated from 
preceding equations are given for saturated liquid 
and gaseous mercury at temperatures from the triple 
point, —38.88° C [22], to 500° C. The numbers in 
parentheses given under the symbols heading the 
various columns indicate the particular equations 
from which the values were calculated. 

Though the tabulated values of the heat capacities 
of the liquid were calculated to apply at the vapor 
pressures, the corresponding values of C,,,) at any 
small fixed pressure are practically identical. For 
it may be readily shown by using eq 13 that C, 
varies by not more than 0.001 percent in this temper- 
ature range when the pressure varies by | atm. 

It will be noted that the free-energy values caleu- 
lated from eq 27 and listed in table 5 apply to eithes 
the saturated liquid or vapor. 





Relative j ( 7 vapo ation; compressi TABLE 5 Entropy; re r apor pressure 


TABLI 
2. Effects of Gas Imperfection and Temperature Scale 


The values listed in tables 3, 4, and 5 contain con 
tributions from the small terms (in the equations 
that contain /, or its derivatives, or 7 and that thus 
correct for gas imperfection of mercury vapor and 
for the deviations between the thermodynamic an 
International temperature scales. When these con 
tributions for a given property vary appreciably wit! 
temperature, they are listed separately in tables | 
and 7 for a number of temperatures. (Except fo 
vapor pressure these are listed in the same units as 
the values in tables 3, 4, and 5.) The magnitudes of 
these contributions depend in most cases on the pat 
ticular way in which the properties have been cal 
culated in this paper, as well as on what basic data 
have been used. Consequently, they are not intrio 
sically characteristic of the properties themselves 

Table 3. There are no contributions to the relative 
enthalpy of the liquid, as the values used are thos 
measured essentially directly The contributions at 
hence identical for the relative enthalpy of the vapo 
and for the heat of vaporization. These are listed | 
tables 6 and 7 

Table 4 The three calculated heat capacities 0 
the liquid have no contribution from gas impertfe: 
tion. At a given temperature they have identica 
contributions from the temperature-scale correctio! 
as listed in table 7. The contribution to C, ,. fro 
gas imperfection is the deviation from 5/2 R (i. « 
from 4.968 cal g-atom”~'deg 
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Contributions to thermod jnamic prope rlies from qas 


TABLE ¢ 
imperfection of mercury vapor 


H Hrs ro S To To vapor 
a ’ F-Hy; pressure 


cal g-atom 
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40 13.8 oO000 7.7 +). OS 
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TABLE 7 Contributions to thermodynamic properties from 


correction to thermodynamic te m pe rature scale 
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Table 5. To the entropy of the liquid the contri- 
bution from gas imperfection is the same at all 
temperatures, +0.0245 cal g-atom™'deg™!. The con- 
tribution to the entropy from the temperature-scale 
correction is practically constant, being equal to 

0.002 cal g-atom~'deg™! at all temperatures in the 
ease of the liquid and not lying outside the range 
0.000 to +0.001 in the case of the vapor. The other 
contributions to the entropy of the vapor, the free 
energy, and the vapor pressure are listed in tables 6 
and 7 


VI. Discussion of Results 


1. Reliability 


An index to the reproducibility, ol “precision”, of 
the enthalpy measurements on liquid mercury is 
afforded by the deviations from the means, as shown 
by the results of the individual measurements, which 
are recorded in table 1. Another index is provided 
by the deviations (also listed in table 1) from the 
smoothed values as represented by the empirical 
equation adopted. All the results lead to an average 
probable error of the mean of about 0.02 percent on 
enthalpy, and a corresponding magnitude of about 


0.1 percent on the derived heat capacity values (C 
or C, «)). 

One check on the over-all accuracy of the appara- 
tus in measuring enthalpy was made by measuring 
the heat delivered to the ice calorimeter by a Monel 
capsule containing water and dropping from 250° C 
By thus determining in several measurements the 
difference in heats for two amounts of water differing 
by about 6 g, a mean value of 1042.05 abs j g™' for 
ale’ of water, an enthalpy function defined else- 
where [23], was obtained. This figure differs by only 
0.02 percent from the value of 1041.85 published in 
the latest report [24] on the thermal properties of 
water as accurately measured earlier in this labora- 
tory by an adiabatic calorimeter. 

The various sources of appreciable systematic 
error were examined in order to estimate their 
likely contributions. As measurements by the 
General Electric Co. [24a] of Nitralloy and of 
Swedish iron in mercury at various temperatures 
indicate a solubility of less than 1 part in 100,000,000 
at 500° C, no error in the enthalpy should be caused 
by the dissolving of the container. Considering the 
uncertainty in each possible error, the authors be- 
lieve that the values of the enthalpy of the liquid 
given in table 3 are accurate to Og percent, except 
below 100° C, where small errorgJin measurement 
become relatively more importan?® as 0° C is ap- 
proached. As a consequence, it is believed that the 
corresponding uncertainties in the heat capacity 
values (’ and (, (table 4) may be as large as 
0.3 percent between 25° and 425° C. Outside this 
temperature range the values given should be con- 
sidered much more uncertain, as they resulted from 
the extrapolation of an empirical function beyond 
the range of experimental measurements. 

The calculated properties other than the enthalpy 
and heat capacity of the liquid are, as pointed out 
earlier in this paper, subject to varying uncertainties 
caused by uncertainties in the corrections for gas 
imperfection and temperature scale. The values in 
tables 6 and 7 should aid in estimating such uncer- 
tainties in specific cases. In many cases the differ- 
ence between the values at two temperatures for a 
given property will be much more accurate than the 
listed absolute magnitude of the property itself. 

The value of the absolute entropy as calculated 
here may be compared with that arrived at through 
use of low-temperature heat-capacity data for solid 
mercury. For liquid mercury at the triple point, 
for example, table 5 gives a value of 16.50 cal 
g-atom™'deg™', based principally on vaporization 
data. If eq 19, representing the vapor pressure 
data of Beattie, Blaisdell, and Kaminsky |9], can be 
accepted as having an accuracy comparable to that 
which they claim, the entropy value just quoted 
should not be uncertain beyond a very few hun- 
dredths of a cal g-atom@'deg™'. This estimate 
includes the consideration of the aforementioned 
uncertainty in the extrapolated heat capacity of 
the liquid between 0° C and the triple point 

Pickard and Simon [25] have recently measured 
the heat capacity of solid mereury down to 3° kK 
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By using their data and those of others, the authors 
have computed a value of 14.34 cal g-atom™'deg 
for the entropy of solid mercury at the triple point 
Addition of the entropy of fusion, 2.38 eal g-atom ‘- 
(26, 27], gives a value of 16.72 for the liquid at 
the temperature. This is 0.2 cal g-atom™'- 
deg™' higher than the value in table 5. The dis- 
crepancy seems to the authors as more likely due to 
errors in the low-temperature heat-capacity data 
than to anv other source 


deg 


Sane 


The effect of certain possible errors on the accuracy 
of vapor pressures as calculated by eq 35 will now be 
examined. As stated earlier, this equation has been 
adjusted to agree with the empirical eq 19 at, and in 
the immediate vicinity of, the normal boiling point, 
which both equations give as 356.58° C However, 
according to the experimental results on which eq 19 
is based, there may be an absolute error of 0.01 deg 
in this temperature, corresponding to an error of 0.02 
percent in the vapor any neighboring 
temperature Equation 19 values of vapor 
pressure at 350° and 362° C (temperatures near the 
extremes of the range of claimed validity of the 
equation) that are higher by 0.004 and 0.001 percent, 
respectively, than given by eq 35. These 
divergences are within the precision of the data on 
which eq 19 is b@ed 

The small uncértainty in the temperature assumed 
for the normal boiling point will affect all vapor 


pressure at 


vives 


those 


pressures calculated from eq 35 by the same percent- 


There four other factors contributing 
significant uncertainty to these calculated pressures 
and as the temperature becomes increasingly higher 
or lower than the boiling point, the effect of each of 
four factors becomes acceleratingly 
These factors are (1 
thermodvnamic 
vaporization at 


age 


are 


these greater 
degree of gas imperfection, (2 
temperature heat of 
the norma! boiling point, and (4 


scale, (3 


y 
° 


+ 


BSERVED MINUS CALCULATED VAPOR PRESSURE , 
+ 





average heat capacity 

boiling point and the 
Taking the uncertainty in the last factor to be | 
percent, as estimated above, it was computed tl 
the resulting uncertainties in the calculated vay 
pressures would be as follows: At 39° C, 0.7 pe 
cent; at 100°, 0.15 percent; and at 250° or 500°, 0.1 
percent.” The assumption of these and reasonab 
uncertainties in each of the other factors has led 

the assignment of the following uncertainties 

individual values of vapor pressure calculated fro: 
eq 35 and listed in table 5: At 39° C, 1.5 percen 
at 100°, 0.5 percent; at 250° or 500°, 0.2 percent; ar 
at 357° (the boiling point), 0.03 percent. The 
figures were arrived at without consideration « 
agreement of the calculated vapor pressures with an 
direct measurements, except those represented b 
eq 19. 


of the liquid between t 
temperature in questi 


A number of investigators [28 to 33] have direct] 
measured the vapor pressure of mercury in th 
temperature range. The deviations of most of thy 
more precise of these experimental values from thos 
given by eq 35 are shown in figure 1. In most of 
these measurements there was an average variatio! 
of several times the uncertainties just stated 


2. Comparison of Liquid Heat Capacity Values With 
Those of Other Experimental! Investigations 


A number of other investigators [34 to 47] hav 
measured the heat capacity of mercury above 0° ¢ 
Most of these results are shown in figure 2 for com 
parison with the results reported in this paper 
The experimental points labeled NBS were caleu 
lated by dividing by the temperature interval the 
differences in the mean unsmoothed experimental! 
heats for pairs of successive furnace temperatures 
A very small correction for curvature was applied 
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as found by various investigators 
/ ) 


Most observers are in substantial agreement below 
100° C, but those who extended their measurements 
to higher temperatures show considerable disagree- 
ment Gaede [42] emploved an essentially adiabatic 
calorimeter; Hirobe [46] used an isothermal one; and 
Dixon and Rodebush [41], claiming an accuracy of 
approximately 1 percent, resorted to adiabatic com- 
pression. The observations of Milthaler [39], Nac- 
eari [37], and Winkelmann [38], all of whom employed 
the method of mixtures, agree approximately, as to 
a value of the negative temperature coefficient, with 
the later measurements of the General Electric 
Co. [45]. The measurements of Barnes and of 
Barnes and Cooke [34] were made, apparently, with 
considerable care by using a continuous flow method 
previously employed in measuring the specific heat 
of water. They considered their measurements 
below 100° to be accurate to 0.1 percent. How- 
ever, above 150° their values differ markedly from 
those of all other observers, including the present 
authors. In spite of this difference, their values 
have often been considered by some in the past to 
have a reliability superior to the heat-capacity values 
of mercury of the other observers 

Recently Kleppa [48] has used an_ electronic 
pulse-circuit technique to measure ultrasonic veloc 
ities in mereury at 50° and 150° C, obtaining a 
precision of about 1 percent. When these velocities 
were combined with the coefficients of thermal 
expansion given by eq 13 and with the values of 
( interpolated from table 4, they gave values 
of ¢ of 5.72 0.01 and 5.49 0.01 cal g-atom™'!- 
'_ respectively, at these two temperatures. The 
corresponding values of C, interpolated from table 
1, based on the direetly measured compressibilities 
of Smith and Keyes [8], are 5.68 and 5.43, respec- 
tively 
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3. Isobaric and Isochoric Heat Capacities of the 
Liquid 


The values of the heat capacities of the liquid C, 
and (, that are given in table 4 are represented 
graphically in figure 3 as multiples of the gas con- 
stant AR. It is of interest that the curve for C, 
shows a minimum at a temperature somewhat below 
the boiling point. This is analogous to the results 
obtained by the authors for two other liquid metals, 
sodium [3] and potassium [49], whose values of C, 
were found to exhibit even deeper minima than 
mercury, and at temperatures somewhat below their 
respective normal boiling points. 

It is well known that the theoretical interpretation 
of C, is simpler than that of C,. The curve for 
C’, of liquid mercury in figure 3, in contrast to that 
for C',, shows no definite evidence of extrapolating 
to a minimum below the critical temperature, 
indicating that the distinet minimum in C, is caused 
by the increasing difference between C, and C, as 
the temperature rises. (Though Smith and Keyes 
[8] measured the compressibility of mercury up to 
300° C, their values above 200° show a very rapid 
and unexpected increase, which would lead to a 
minimum in the calculated curve for C, above this 
temperature. However, they discredited their re- 
sults above 200° because of invalidating experi- 
mental circumstances 


It is of interest that the value of C, calculated in 
this paper (table 4) for liquid mercury at its triple 
point is only 0.1 percent different from 3R, which 
is the value predicted at this temperature by a 
recent fairly successful theoretical treatment of 
liquid mereury [50]. According to some recently 
published theories [51], a liquid exhibits a fairly 
continuous transition from the crystalline to the 
gaseous state as the temperature rises. The heat 
capacity C, of a monatomic liquid, taken at small 
pressures, may be expected to vary from approxi- 
mately 3f at the triple point and approach R 
near the critical point. The decrease in C, as the 
temperature rises may be interpreted as being due 
te a decreasing contribution to disorderliness in the 
liquid. The curve for C, of mereury shown in 
figure 3 is consistent with this picture 


4. Vapor Pressure 


Near the normal boiling point, 357° C, the terms 
for gas imperfection have very little influence on the 
values calculated for the vapor pressure from eq 35 
However, as pointed out in an earlier section, the 
effect is many times as great at much higher or lower 
temperatures. If the calculated vapor pressures at 
such temperatures agree with corresponding meas- 
ured values, there is thus afforded an independent 
experimental check on the adopted value of the 
principal parameter determining the gas-imperfec- 
tion corrections in the various equations for thermo 
dynamic properties. This assumes, of course, that 
the experimental values of vapor pressure have 
sufficient accuracy and precision. 
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The vapor pressure of mercury has been measured 
workers [52] from the triple poimt to 


by various 
Though it very low and very high tem- 


400° C 


peratures in this range the percentage accuracy of 
the vapor pressure measurements need not be very 


great, the precision of most of these results Is so 
relatively poor as to exclude their use for this test 

In 1910 Smith and Menzies [53] carefully measured 
the vapor pressure of mercury from about 250° to 
(35° C 0.1 percent The 
temperatures were recalculated by Menzies [52] in 
1927 on the basis of the new value for the sulfur 
boiling pomt. Menzies added three more experi- 
mental pomts between 120° and 200° C and derived 


the following equation to represent the data: 


with a precision of about 


3285.92 


log r 0.665240 log A’. 


9.957094 ‘ 
ad 


Ww here 


t’ being deg C Int as measured by the platinum 
thermometer of Smith and Menzies 

The authors have recalculated Menzies’ 
grade temperatures, as given by eq 45, In accordance 
with the present knowledge of the platinum-ther- 
mometer scale. Smith and Menzies [54] had re- 
ported for their thermometer a 6 value of 1.6147, on 
the basis of their assumed sulfur boiling point of 
145° C Correcting this to the basis of the tempera- 
ture assumed by Menzies in 1927, 444.6°, yields a 
value of 6 of 1.5919. Equation 42 gives a mercury 
normal boiling point of 356.711° C, which ts 0.13 deg 
higher than reported by Beattie, Blaisdell, and 
Kaminsky [9]. Waidner and Burgess [55] and 
Beattie, Blaisdell, and Kaye [56] have found inde- 
pendent evidence that a platinum thermometer 
accurately calibrated at the ice, steam, and sulfur 
points indicates for the mercury boiling point a 
temperature depending fairly systematically on the 
5 value of the thermometer In fact, Beattie et al 


centi- 





have argued that by adding the mercury boili 
point as a fourth calibration point and using a cul 
temperature-resistance relation, the platinum-th 
mometer scale will be approximately independent 
the thermometer constants. 

This recommendation has been followed here 
cubic temperature-resistance equation has bi 
derived that gives the same temperatures for 
ice, steam, and sulfur points but a value of 356.58° ¢ 
for the mercury boiling point when there are subst 
tuted the resistances given for these four temper 
tures by a Callendar equation with 6=1.5919, a 
with the mercury boiling point taken as 356.711° 
Equating the resistances given by the two equatio: 
indicates that the Centigrade temperatures giv: 
by eq 42 should be changed by the amounts show 
in table 8 to accord with a four-point calibratior 
Though a fourth calibration point is obviously, 
itself, an asset toward greater accuracy, it is believ: 
that considerably more uncertainty should | 
attached to the interpolations of temperature pr 
vided by a thermometer with such a bigh 6 valu 
than to those by thermometers meeting present-da\ 
standards 

After correcting the Centigrade temperatures of 
eq 42 by the values of table 8, figure 4 was obtaine: 
for the differences between the experimental vapor 
pressures of Smith and Menzies and those calculated 
from eq 35. The four curves represent the diffe: 
obtained depending on what 
energy is made the basis of calculating the gas-im 
perfection terms in the latter equation. This equa 
tion has been derived to give the experimentally 
measured normal boiling point 
magnitudes of the gas imperfections assumed. Fo 
a range of 200 deg below the boiling point, the best 
agreement between experimental and = calculated 
vapor pressures Is seen to correspond to a gas im 
perfection equivalent to a dissociation energy of | 
keal mole~' of Hg. This value was claimed earlie: 
in this paper to be a reasonable weighted mean o! 
those indicated by independent spectroscopic ev! 
However, in of the smallness of thi 


ences dissociatior 


regardless of thi 


dence view 





Fieure 3 


He ul capactires of liquid mercury al constant pressure and constant 


v rl “nm 
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deviations of figure 4 and the uncertainties in some 
of the small corrections applied to both experimental 
and calculated vapor pressures, the agreement may 
be somewhat fortuitous. 


TABLE 5 Correction of temperature values of a platinum 
thermometer with 6 1.45919 to accord with a four-point 
calibration in which the merc ury howling point 18 changed 


from 356.711° to 356.48° ( 


Premperature Correction 
C Int C Int 
1 0. 000 
1m s 
200 O78 
250) 11 
we 1S 

) 140 
+). SS 131 
wn OsY 
450) O34 
44.6 Ooo 


It will be noted that above the boiling point all 
four curves indicate lower experimental than cal- 
culated vapor pressures. This may be attributed 
quantitatively, in large part, to the fact that as the 
temperature rises the saturated vapor 
denser, and as a result the empirical eq 42, derived 
to fit the data well near the boiling point and at 
lower temperatures, is too simple to represent the 
increasing importance of gas imperfection and 
liquid volume 

The ordinate differences in figure 4 are of the same 
order of magnitude as the absolute uncertainties 
assigned in an earlier section to the vapor pressures 
calculated from eq 35. However, the several factors 
named there as affecting the accuracy of this equation 
have effects that vary in roughly comparable ways 
with temperature. Therefore, errors in the adopted 
magnitudes of their effects on the calculated vapor 
pressures at various temperatures would be capable 
of being compensated considerably by the choice of 
a single somewhat erroneous value for the dissocia- 
tion energy of Hg, for calculating the effects of gas 


becomes 





imperfection. It is thus possible that the choice 
of 1.5 keal per mole (made partly on the basis of the 
apparent agreement in fig. 4 between observed and 
calculated vapor pressures when this value is selected 
is of this nature. 

Nevertheless, it is believed that the comparison 
afforded by figure 4 provides confirmatory evidence 
that the values of dissociation energy and second 
virial coefficient of mercury that were selected in 
this paper are not far from the correct ones. That 
this comparison between observed and calculated 
vapor pressures has such significance is due in no 
small degree to the accuracy of the experimental 
values of liquid heat capacity recently measured and 
reported in this paper. For a given temperature, 
an error in the heat capacity produces an approxi- 
mately proportional error in the vapor pressure 
calculated in this manner. The previously available 
heat capacity values, because of their disagreement 
in the region above and below the boiling point, 
would no doubt have been considered up to 10 times 
as uncertain Therefore had it been necessary to 
rely on these previous values, the relatively small 
differences of such a comparison graph as figure 4 
would have had very much less significance 


The heat-capacity 


measurements were greatly 
expedited by 


the commendable cooperation ol 
Leo F. Epstein and his associates, of the Knolls 
Atomic Power Laboratory, Schenectady, N. Y 
Under his direction the mercury samples were pains- 
takingly prepared and supplied in a state of very 
high purity and ready for the thermal measure- 
ments. Collaborating with him in this work were 
L. W. Hibbs, Jr., who purified the mercury, filled 
and sealed the containers, and tested them: George 
Strichman, who tested the containers for tightness; 
and R. E. Schofield, who performed the 
spectrometer analyses. 

The authors express their special indebtedness also 
to certain members of this Bureau. B. F 
and = his performed a 


Mmass- 


Seribnet 


associates spectrochemical 
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